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1 Introduction 



During the 1950s an intense search for new elementary particles was performed at various 
particle accelerators which led to the discovery of several hundreds of hadron states. Based 
on their spin hadrons were classified as mesons (spin = 0, 1) or baryons (spin = 1/2, 3/2). 
The richness and variety of the hadronic spectrum made it likely, that hadrons are not 
fundamental particles. First attempts to classify the hadronic spectrum by properties as 
mass, charge, isospin and flavor were undertaken by W. Heisenberg, M. Gell-Mann and 
different Japanese groups. In 1961 M. Gell-Mann and Y. Ne'eman independently proposed 
hadron representations based on the flavor group SU{3), [1-3]. The mesons and spin 1/2- 
baryon states can be grouped into octets and the spin 3/2-baryon states into dccuplcts of 
similar properties. The prediction of the unknown baryon fl~ with strangeness —3 and 
electric charge —1 was a triumphant success of this theory. A mathematical description 
of the hadron states was introduced in 1964, when M. Gell-Mann [4], and G. Zweig [5] 
proposed the static quark-model due to which hadrons are bound states out of three quarks 
in case of baryons and a quark-antiquark pair for mesons. Here the quarks are assumed to 
be spin 1/2-fermions. Quantum numbers of all hadron states known by that time could 
be described assuming three different quark flavors : up (m), down (rf), and strange(s). 
The quarks had to carry fractional elementary charges = +2/3, ea — eg — —1/3. The 
Q~-baryon, which is composed out of three s-quarks, was obtained as a fermionic state 
with a spin and flavor-symmetric wave function. Since quarks are fermions, this is not in 
accordance with the spin-statistics theorem [6-8]. This contradiction was resolved, after 
some intermediary steps, by introducing the three-valued quantum number color [9-12] 
for quarks. Since colored particles could not be observed experimentally, it was proposed, 
that all physical hadron states are color-singlets, [10, 11, 13]. 

Because of their anomalous magnetic moments, protons and neutrons were suspected 
to be composite states since the 1930s, cf. [14, 15]. In a series of lepton-nuclcon scatter- 
ing experiments Hofstadter and collaborators [16-19] have shown during the 1950s, that 
baryons possess extended charge distributions irrespective of their net-charge, which was 
a further clear sign of compositeness. 

In the late 1960s deeply inelastic lepton-nucleon scattering (DIS) experiments at the 
Stanford Linear Accelerator SLAC, [20,21], made it possible to study the substructure of 
nucleons at much higher spatial resolution. One may parametrize the differential scat- 
tering cross section in terms of different non-perturbative nucleonic structure functions 
Fj, [22,23]. They describe the structure of the corresponding hadron at a given energy- 
transfer p in its rest frame and the virtual 4-momentum-transfer from the lepton to the 
nucleon = —Q^. The spatial resolution reached is Ax ~ Based on current 

algebra techniques J. Bjorken predicted, that in the limit, Q^, u ^ oo, Q^/i^ = fixed, 
the structure functions do not depend on and i/ independently but on the Bjorken 
variable x — Q'^/2Mv, [24]. Here M is the mass of the hadron and < x < 1. This 
behaviour has been discovered in the SLAC-MIT experiments for nucleons, [25-28] cf. 
also [29-31]. R. Feynman provided a phenomenological explanation of this behaviour 
by creating the parton model, [32-34]. The strong correlation between and v in the 
Bjorken limit can only be understood if the nucleon is composed out of point-like con- 
stituents, Feynman's partons, at short distances. Partons appear as basically massless 
quantum fluctuations which live long cnoTigh to be resolved during the short interaction 
times Tint ~ l/?o, cf- [35], implying strict conditions for the validity of the parton model. 
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The hadronic structure functions can then be understood as incoherent weighted sums 
over individual parton distribution functions fi{x) at lowest order. The discovery, that 
the longitudinal structure function Fi vanishes in the Bjorken limit, known as Callan- 
Gross relation [36], showed that the quark-partons possess spin 1/2. Feynman's parton 
model and Gell-Mann's quark model were finally joint by Bjorken and Paschos [37], who 
identified the partons as the quarks of the group theoretic approach. 

In the 1950s C.N. Yang and R.L. Mills, [38], studied gauge field theories based on non- 
abelian gauge groups. The non-abelian nature implies three- and quadrilinear couplings 
between the gauge bosons contrary to the case of Abelian groups, like in Quantum Electro- 
dynamics. The renormalizibility of field theories of this type was proven by G. 't Hooft [39] 
in the massless case. In 1972/1973 M. Gell-Mann, H. Fritzsch and H. Leutwyler, [40], cf. 
also [11], proposed to describe quark interactions by a SU{3)c Yang-Mills theory. Quan- 
tum Chromodynamics (QCD) has been introduced as the dynamical theory of quarks and 
color-octet vector gluons as gauge bosons. In 1973, D. Gross and F. Wilczek, [41], and 
H. Politzer, [42], proved that Quantum Chromodynamics is asymptotically free, i.e. the 
strong coupling decreases toward shorter distances or larger momentum transfer. This 
allows to perform perturbative computations for scattering processes at large enough 
momentum scales. 

A useful tool to apply perturbative methods in deeply inelastic scattering processes 
is the light-cone expansion (LCE) [43-46]. Applying it to deeply inelastic scattering 
processes leads to a factorization theorem, which separates hadronic bound state effects 
from short distance effects allowing for a perturbative analysis of the operators at the 
respective twist-level. The long-range bound state effects are collected in the parton 
distribution functions, which have to be determined from experimental data or com- 
puted using non-perturbative methods. The short-distance effects are associated with the 
Wilson-coefficients and allow for a perturbative analysis. 

In 1975 logarithmic scaling violations in the deep inelastic cross section were discov- 
ered, [47,48], which are due to the fact that QCD is not a free field theory neither is it con- 
formally invariant, [49]. These scaling violations were predicted perturbatively in higher 
order calculations and constitute a great success for the Standard Model, [50,51]. Beyond 
the established light-quark picture a fourth quark, charm, has been proposed [52-54]. It 
was discovered in 1974 almost simultaneously at SLAG and at Brookhaven National Labo- 
ratory (BNL) through the J/tp meson. In 1977 the T resonance was observed, [55], which 
was interpreted as a meson consisting of a new quark, the bottom quark (b). The search 
for quarks continued and in 1995 the top quark was discovered at TEVATRON, [56,57]. 

Today QCD is established as the theory of the strong interaction and together with 
the electrowcak SUii'i) x f/y(l) sector it forms the Standard Model of elementary par- 
ticle physics. The electroweak sector of the Standard Model has been introduced by 
S. Glashow, [58] and S. Weinberg in 1967, [59], cf. also [60,61]. G. 't Hooft and M. Velt- 
man proved that this theory is renormalizable, [62] , see also [63-65] . 

Deep- inelastic scattering provides a clear method to probe the short distance substruc- 
ture of hadrons in the space-like domain. Many DIS experiments have been performed 
through the last forty years, [66-71]. The proton substructure has been examined in- 
tensely at HERA at DESY [72-76]. In case of unpolarized deep-inelastic scattering via 
single photon exchange the cross section is described in terms of the structure functions 
F2{x,Q'^) and Fl^XjQ"^). While F2 has been measured over a wide kinematic range, F^ 
has mainly been measured at fixed targets. Experimental data from HERA showed that 



6 



the structure functions receive a substantial part due to charm quark pair production, 
beyond the hght flavor contributions which amounts 25-35 % in the small x region, cf. 
e.g. [77-80]. 

The scaling violations of the DIS structure functions are determined by the anomalous 
dimensions and the Wilson coefficients both for the light partons and heavy flavors. The 
anomalous dimensions have been computed at leading order (LO), [50,51], and next-to- 
leading-order (NLO) [81-88]. At next-to-next-to-leading-order (NNLO) fixed moments 
were obtained in Refs. [89-92] before the general result in the Mellin variable N was 
determined in 2004 [93,94] by Vermaseren et al. A first independent check on the NNLO 
moments of the unpolarized anomalous dimensions was performed in Ref. [95]. The 
massless Wilson coefficients for the structure functions F2 and F^, resp. their moments, 
were computed at first order in [96-98], at second order in [99-107], and at third order 
in [89-92,108,109]. 

Due to the large size of the charm quark contribution to the DIS structure functions 
the precise computation of these terms is very important. The leading order massive 
Wilson coefficients have been computed completely in the late 1970s, [110-114]. The 
NLO-corrections are available in semi-analytic form, [115,116]. A precise numerical im- 
plementation was given in [117]. At three- loop order the computation of the heavy flavor 
Wilson coefficients over the whole kinematic range appears to be rather difficult at present. 
However, a very important region for deep-inelastic experiments at HERA can be covered 
by computing the massive Wilson coefficients for the structure function F2 in the limit 
of ^ m^. There the heavy flavor Wilson coefficients factorize, cf. Ref. [118], into 
the process dependent light flavor Wilson coefficients C(q^g)^(2,L){x, Q'^ / n'^) and the pro- 
cess independent operator matrix elements (OMEs) Aij{x, jj? /w?). The OMEs contain 
all mass dependence and are obtained as matrix elements of the leading twist local com- 
posite operators between partonic states | j) = q,g). For the structure function 
F2{x,Q'^) this representation becomes effective for Q^>10 m^, cf. [118]. In the case of 
the longitudinal structure function Fl this factorization theorem applies for much larger 
momentum transfers, (5^~800 m^, only, [118], which lays outside the kinematic region 
probed at HERA. The NNLO corrections for the longitudinal structure function have been 
computed in Ref. [119]. 

At two-loop order the quarkonic OMEs Agj have been calculated analytically in 
Ref. [118] and conffimed, applying rather different methods, in Ref. [120]. The contribu- 
tions linear in the dimensional parameter e were calculated in Ref. [121]. The remaining 
gluonic matrix elements Agj have been derived in Ref. [122]. They were confirmed and 
extended to 0{e) in [123]. The gluonic operator matrix elements are required to describe 
parton distribution functions in the variable flavor number scheme (VFNS). Furthermore 
they contribute to the NNLO quarkonic singlet OMEs through renormalization. The 
complete renormalization procedure for massive 3-loop OMEs was developed in Ref. [95] . 
There also the general structure of the NNLO OMEs was derived. In Ref. [95] a large 
amount of Mellin moments for all contributing massive OMEs were computed. All log- 
arithmic terms oc In'' {Q^ / m^) , k = 3,2, 1 arc known in full analytic form by now. The 
mathematical structure of these results is determined by nested harmonic sums [124, 125] . 
Previous analyses of known results of different single-scale hard scattering processes have 
shown, that at least up to massless three-loop order calculations, single scale results are 
most simply expressed through nested harmonic sums, [93,94,109,126-130]. They obey 
algebraic, [131], and structural relations, [132,133]. If one considers fixed moments only. 
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one obtains representations in terms of multiple zeta values, [134, 135]. 

In the present computation Feynman diagrams are evaluated by direct integration 
leading to a representation in terms of generalized hypergeometric functions [136,137], cf. 
also [132]. In the case of massless computations summation algorithms as in Refs. [125, 
138, 139] can be applied. In the massive case, various additional infinite and finite sums 
occur which possess a much more involved structure. These sums can be treated applying 
modern summation technologies encoded in the package Sigma, [140-143], written in 
MATH EM ATI CA [144]. 

In Rcfs. [95,118-123] all necessary calculations to describe the massive Wilson coef- 
ficients in the asymptotic region to 2-loop order and to determine all quantities which 
have to be known to renormalize these quantities at 3-loop order have been performed. 
Thereby all logarithmic contributions are known at general values of the Mellin variable 
N. For an essential piece ^ in the constant part of the renormalized OMEs the general N- 
dependence is not known yet. A series of Mellin-moments has been computed in Rcf. [95]. 
The general result is of numerical importance and has therefore to be computed exactly. 
This applies also to its small-x behaviour which may cause large effects. 

In this thesis we perform a first step within this larger programme and compute the 
0{o?gTpnf) contribution to the massive operator matrix elements Aqg, ^q^, A^q,Qi A^q,Q 
and ^qq^Q for general values of the Mellin variable A^. Due to the large fraction of the 
heavy flavor contributions to the decp-inclastic structure functions the precise knowledge 
of the respective Wilson coefficients is of essential importance to consistently derive the 
parton distribution functions at leading twist for the gluon, the valence quarks, and for 
the different sea-quark species, along with a precision measurement of the strong coupling 
constant q;s(M|), cf. e.g. [146]. as(M|) in itself is one of the fundamental quantities in 
nature and has to be known as precisely as possible. Note that, despite of a large number of 
precision analyses using different high-energy observables, a final agreement on the value 
of q;5(M|), being measurable with an accuracy of ~ 1% at present, could not be obtained 
yet, [147]. Needless to say that all measurements at the Large Hadron Collider LHC at 
CERN crucially rely on both the precision knowledge of the parton distribution functions 
and q;s(M|). Furthermore, as{M'^) is an essential input-parameter in scenarios of the 
potential unification of the fundamental forces of the strong-, weak-, and electromagnetic 
interactions, cf. [148,149]. Its value is decisive for the question whether, and in which 
theory, the fundamental forces of nature unify at high energy scales or not. This, in 
turn, touches the respective scenarios of the physics in the early universe, and is thus also 
connected to the major challenging problems in physics. 

The outline of this thesis is as follows. In Section 2 we describe the basic high-energy 
process, deeply-inelastic lepton-nucleon scattering, to which the QCD corrections, which 
are calculated, belong. This includes the QCD-improved parton model at short distances, 
which is established through the light-cone expansion. The heavy flavor contributions to 
the decp-inclastic structure functions can be viewed as a linear contribution in addition 
to the light parton contributions. The leading order corrections in the strong coupling 
constant are re-calculated in Section 3. We then discuss the general scenario which allows 
the analytic computation of the heavy flavor corrections to higher orders in the asymp- 
totic region 3> in terms of massive OMEs and the massless Wilson coefficients. 
The formalism is outlined to 3-loop orders in Section 4. In Section 5 the details of the 
renormalization of the massive OMEs are summarized to 3-loop orders. The leading or- 

^For numerical studies of this aspect see [145]. 
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der massive operator matrix element is then re-calculated in Section 6 as an introductory 
calculation, in which we discuss the asymptotic factorization theorem for the heavy flavor 
Wilson coefficients following the formalism of Section 4 and compare to the explicit cal- 
culation in Section 3. Wc then turn to the computation of the analytic 0{a^gTpN fC f,a) 
contributions to five heavy flavor OMEs at general values of the Mellin variable in 
Section 7. We describe the contributing Feynman diagrams, give details of their analytic 
evaluation, and illustrate the methods by different explicit examples. In Section 8 we 
present the analytic results for the OMEs Ag], AJJ^^^^^ and A^^^J^^'^^. The 

new results arc the constant parts of the respective unrenormalized 3-loop OMEs and a 
first independent recalculation of the corresponding contributions to the 3-loop anoma- 
lous dimensions. This calculation generahzes results obtained for fixed integer moments in 
Ref. [95] to general values of N, required by the experimental analyses, for the first time. 
The final results arc expressed in terms of nested harmonic sums, although intermediary 
results require the treatment of generalizations thereof, which finally cancel for the class of 
graphs computed. In Section 9 we summarize the main results. A series of technical details 
of the present computation is given in the Appendix. Basic conventions are summarized 
in Appendix A. A consistent set of Feynman rules for QCD, including those of twist-2 
composite operators, is given in Appendix B. Relations for D-dimensional momentum 
integrals are summarized in Appendix C. In Appendix D the results for the individual 
diagrams, which were computed in this thesis, are given for reasons of documentation, to 
allow other groups for comparison, and to discuss their mathematical structure. Useful 
variable transformations are summarized in Appendix E. The present computation relies 
on the Tise of special higher transcendental functions, which allow a particularly compact 
treatment. Main results and relations of these functions, as Euler-integrals, generalized 
hypergeometric functions, as well as harmonic sums, are given in Appendix F. The final 
results of the present work were obtained summing multiply nested sums of the hyper- 
geometric type and their extensions, which are of a sophisticated nature. They could be 
uniquely solved using general modern summation technologies encoded in C. Schneider's 
programme package SIGMA [140-143]. As an illustration for several thousands of sums 
which had to be computed we show a few examples in Appendix G. Partly they contain a 
large amount of also generalized harmonic sums. In Appendix H reference values for the 
moments of the 3-loop anomalous dimensions and the various constant parts af^ of the 
massive OMEs, being computed in [95], are summarized. These values were used to test 
the results of the present calculation. 
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2 Deeply inelastic scattering 



Deep-inelastic scattering (DIS) denotes the scattering process of highly energetic leptons 
off hadrons, and provides a very precise method to probe the substructure of hadrons at 
short space-like distances. The 4-momentum transfer = —Q^ is at least of the order 
> 4GeV^, such that space-like distances of approximately can be resolved. 

Different deeply inelastic scattering experiments exploring charged and neutral current 
reactions allow to probe the flavor and the gluonic structure of the hadron. By performing 
polarized scattering experiments, also the spin-structure of hadrons can be investigated. 




Figure 1: Schematic diagram of deeply inelastic scattering via single boson exchange 
2.1 Kinematics 

The schematic diagram of deep- inelastic scattering at tree level is shown in Figure 1. A 
lepton with momentum / is scattered of a nucleon with momentum P and mass M via a 
virtual vector boson. In this process the nucleon state disintegrates and F denotes a new 
linear combination of hadronic final states with allowed quantum numbers. The momenta 
of the outgoing lepton and hadronic states are denoted by /' and Pp, respectively. The 
4-momentum q of the virtual vector boson is space-like and the virtuality is defined 

by 

g2 ^ _^2^ q = l-l'. (1) 

Two more independent kinematic variables are sufficient to determine the scattering pro- 
cess: 

s ^ {P + l)\ (2) 
= (P + qY = Pi . (3) 

Here s denotes the total center of momentum energy squared and W is the invariant mass 
of the final hadron state {Pp |. The kinematic variables can be measured from the final 
lepton or hadronic states, depending on the specific experiment, cf. e.g. [150-152]. In 
the following analysis we will neglect the lepton mass and describe the process using the 
following Lorentz-invariant kinematic variables: 

^ Pg _ W^ + Q^- 
^ ~ M ~ 2M ' ^ ^ 

^ 2P.q 2Mv W^ + Q^-M^ ' ^ ^ 

P.q _ 2Mv _ W^ + Q^-M'^ 
^ " p7 ~ S-M2 ~ S-M2 ■ 
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Here v is the total energy transfer in the rest frame of the nucleon, x is a Bjorken vari- 
able and y is the inelasticity, cf. [153]. In general, the momentum q is transferred via 
the exchange of a 7, Z, VF^-boson. In this thesis we limit the investigation to photon 
exchange in unpolarized charged lepton-nucleon scattering. For not too large virtuali- 
ties, i.e. (5^<500 GeV^, single photon exchange dominates this reaction, cf. [154]. Thus 
from now on this region will be considered only and weak gauge boson effects caused by 
the exchange of a Z-Boson may be disregarded. The invariant hadronic mass obeys the 
condition 

W'^>M'^. (7) 

Prom (7) one obtains 

^ (P + qf = + Q2 - 1^ > . (8) 
Thereby x is limited to the region 

< X < 1 . (9) 

For X = 1 the process is elastic, whereas x < 1 describes the inelastic region [155]. By 
considering the nucleon's rest-frame and demanding a positive energy-transfer, we obtain 
further restrictions on the kinematic variables: 

i^>0, 0<y<l,s>M^. (10) 

The cross section for ep-scattering is determined by the transition matrix element for the 
electromagnetic current. In Born approximation it is given by 

Mf, = eMl',v')ru{l,v)^{PF I I P^^) ' (11) 

cf. e.g. [156-158]. Here i]{ri') and cr denote the spin components of the leptons and 
gluons, respectively. The initial and final hadron states are denoted by | P, a) and {Pp \ , 
respectively. 7^^ denotes the Dirac- matrices and u{u) are the bi-spinors of the electron and 
its conjugate, respectively, see Appendix A. Furthermore e denotes the electric charge 
and J^"^(0 is the quarkonic part of the electromagnetic current operator: 

JliO = MO ■ (12) 

For QCD, the electromagnetic current is given by 

^r(o = E^/(^)^^^/7'*/'(^)' (13) 
/,/' 

where ^E'/(0 denotes the quark field of flavor /. X'jJ, describes the electromagnetic charges 
of the different quark flavors. For three light quark flavors it is given by 

1 1 
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where the X^iavor Gell-Mann matrices of the flavor group SU (3) flavor, cf- [159, 160] . 

The unpolarized cross section is obtained by averaging over leptonic and hadronic spin 
degrees of freedom. The differential cross section, cf. [156-158,161], reads: 

Inserting the transition matrix element shows, that the cross section can be decomposed 
into a tensor Lf^i, depending only on the leptonic states and a purely hadronic tensor W/^i, 
with 

L^l^n = 5^ V)7^«(/,r/)] *[«(/', r^OT'^wl/, r;)] , (16) 

ri',r] 

W,.{q,P) = ^J2(2n)'6\PF-q-P){P,a\r;-{0)\PF){PF\Jr{0)\P,a). 

(17) 



47r 

cr,F 



In terms of these quantities the cross section reads 

= mp^'^^ = w^)i''"'^'- '''' 

Here a is the fine structure constant. The leptonic tensor can be computed easily by 
applying the conventions in Appendix A. One obtains 



(19) 



The hadronic tensor cannot be evaluated purely perturbatively due to the non-perturbative 
nature of the matrix elements^. Using the integral representation of the 5-distribution 
and applying elementary quantum mechanical identities, Eq. (17) can be rewritten as, 
cf. [158,168], 

w,.{q,p) = ^J^j d'^^wmip I [jr(o,jrm i n 

(7 

= i- y d\eM^<iO{p I [^r(o, jrm i p) ■ m 

Here the bracket [a, h] denotes the commutator of a and h. The hadronic tensor obeys 
various symmetry and conservation laws, cf. [169]. These impose conditions on the Lorentz 
structure of the hadronic tensor and allow to parametrize it by different scalar structure 
functions. They contain all information about the structure of the proton. In the general 
case 14 independent structure functions exist, [170, 171], but in the case of unpolarized 
DIS via single photon exchange only two structure functions contribute. Here the leptonic 
tensor (16) is symmetric. Since any tensor of rank-2 can be decomposed into a symmetric 
and an anti-symmetric part, only the symmetric part of the hadronic tensor contributes. 
Thus the hadronic tensor must be a linear combination of the following tensors 

(21) 



^Ab initio calculations would have to be based on lattice QCD methods. During the last years, an 
increased numerical precision has been achieved in this field, cf. e.g. [162-167] 
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Prom the conservation of the electromagnetic current, 



Lorentz- and time-reversal invariance it follows that 

q^W^" = . 



(22) 



(23) 



Furthermore strong interactions preserve CP-invariance, cf. [172]. Making a general ansatz 
in terms of (21) and imposing gaiigc invariance leads to the following representation of 
W^i,, containing the two structure functions and F2 



^^P^^^^P^.-^^g^}jFM^). (24) 



2x 



Due to the hermiticity of the hadronic tensor, the structure functions F2{x,Q'^) and 
Fl{x,Q'^) are real functions. Their arguments are Bjorken-x and Q^, whereas in the 
elastic case the cross section is determined by one kinematic variable, e.g. the total 
energy transfer, only. The differential cross section in terms of the structure function is 
obtained by inserting (24) into (18): 



da 



dxdy xyQ^ 



l+{l-yy F2{x,Q')-y'F,^{x,Q') 



(25) 



The two structure functions F2 and Fl can be extracted from (24) by applying the fol- 
lowing projectors in D dimensions: 



ix"- 



2x 



Arv.2 

{D - l)—P^^P^W,,{q, P) - g'^'-W.^iq, P) 



D-2 

Here and in the following we neglect target mass corrections and thereby set P'^ — 



(26) 



2.2 The Parton Model 

In general, the structure functions in (24) depend on two kinematic variables and x. 
However, in the Bjorken limit [24] ^ inf , x fixed, the structure function depends 

on X only. 



lim F(2,L)(2;,g') = F(2,L)(a;) , 

XQ , v) — >• 00, x=const. 



(27) 



which is called Bjorken scaling. Experimental observations from electron-proton collisions 
performed at SLAC in 1968, [25-28], confirmed the existence of an approximate scaling 
behaviour and thereby supported Bjorken's predictions. Furthermore it was shown ex- 
perimentally, that the cross section remained large at high momentum transfers Q^. This 
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behaviour indicates point-like particles in the target, as no further substructure could be 
probed with increasing momentum. The earlier picture according to which the size of the 
proton was about about 10~^^ cm with a smooth charge distribution, [16,18,19], which 
is valid at lower momentum transfer, was superseded. Feynman solved this evanescent 
contradiction by introducing the parton model, [32,33], cf. also [34,37,156,157,173,174]. 
At large enough scales the proton is a composite object, consisting of several point- like 
particles, the partons. During the interaction time with the virtual photon, the long-lived 
partons behave as quasi-free particles. The photon scatters elastically of a single parton, 
while the other partons act as "spectators" and do not interfere with the process. Thus 
the total cross section is given by the incoherent sum of the individual parton-photon 
cross sections, weighted by the probability to find a specific parton i with a momentum 
fraction z inside the proton. Feynman introduced this probability as parton distribution 
function (PDF), fi{z). From now on we will use the collincar parton model, according 
to which the momenta of the specific parton p is taken to be coUinear to the nucleon 
momentum P, 

p^zP . (28) 
Analogously to the scaling variable x one may define a partonic scaling variable r, 

T=^. (29) 

Combined with (28) one obtains 

Tz = x . (30) 

Feynman's original approach, the naive parton model, neglected the radiative corrections. 
Its main issue consists in the strict correlation 

^(^-^V (31) 

This condition implies z — x. Furthermore, according to this model protons are always 
composed of two u and one d valence quarks. With the advent of QCD this model was 
modified and also virtual quark states and gluons were incorporated as additional partons. 
This more advanced model is known as QCD-improved parton model and can be derived 
by applying the light-cone expansion. Here the assumption is made, that the hadronic 
tensor factorizes into the parton distribution functions (PDFs) and a partonic tensor 
yyi : 



(32) 



The partonic tensor is given by (20), where the hadronic states (P | are substituted with 
the corresponding partonic state (p \ of the struck parton. fj{z) denotes the PDF of the 
respective anti-parton to parton i. Assuming that the electromagnetic parton current 
takes the following form 

I I ^) = -teilfj.u' , (33) 
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where is the charge of the respective parton i, one obtains 



27re^ r i 

K^^ir, Q') = -^Sil - r) [2p>t + P;?. + pIq, - M-P'\ ■ (34) 



Combining the ^-distributions in (34) and (32) leads to Fcynman's assumption of the 
naive parton model: z = x. Applying the projectors (26) to the hadronic tensor (32) 
yields the following structure functions at lowest order : 

Fl{x,Q'') = 0, 

F,{x,Q') = xJ2enfi{^) + fM] ■ (35) 

i 

The parton distribution functions are determined from DIS- world data analyses by differ- 
ent groups. Currently they are known to NNLO, i.e. at 0(q;^), in the unpolarized case, 
cf. Refs. [146,175,176]. 
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3 Calculation of the Wilson coefficients H)Jr\ 

In the following we calculate the leading order massive Wilson coefficients contributing 
to the structure functions F2{x,Q'^) and Fl{x,Q^) for pure photon exchange. It is given 
by the Bethe-Heitler fusion process of virtual photon-gluon scattering: 

^*+9^Q + Q (36) 






Figure 2: Feynman diagrams contributing to the leading order heavy flavor Wilson coefficient 

The Feynman diagrams for the lowest order contributions to this process are shown in 
Figure 2. The calculation is performed referring to the collinear parton model, i.e. the 
partons are assumed to act like free particles during the interaction and their momentum 
is parallel to the proton momentum, 



k = zP . 



(37) 



Here P denotes the proton momentum and k is the gluon momentum. The cms velocity 
V of the outgoing heavy quarks of mass m is then given by 



1-4 



, ra T 



1/2 



(3^ 



where denotes the negative squared momentum of the space-like photon and r is 
defined as the ratio of the Bj or ken- variable x and the momentum fraction z. 



X 



T :- 



(39) 



Applying the Feynman rules, cf. Appendix B, to the diagrams in Figure 2 yields the 
matrix element 

ffi — d + m _ ffi — ft + m 

M^y = u{p2)^gsl^,i- ^2 ::i<^qluv{pi) + u{p2)igs-ff,i- -iegjuv{pi) . 

[Pi -q) [Pi -k) -m^ ^^Q^ 

The massive Wilson coefficients are obtained as projections of the squared matrix element 

Here the sum over the gluon polarization states is performed by contraction with —g^p- 
The average over the spin directions gives an additional factor of 1/(D — 2), with D the 
space-time dimension. The massive Wilson coefficients for D = A read 



16 



(1)/ rrl 



2 



(42) 
(43) 



Here Tp is a color factor which is given by 1/2 for SU{N), see Appendix A. The phase 
space integral J dR2 in this case is given by, see e.g. [156] , 



/ 



dR2 



dt , 



(44) 



where 



Pc 



= -i^AV2(,,o,-g^) = -^AV2(,,^^^2) 



(45) 



denotes the modulus of the center-of-mass 3-momentum of the initial (final) state particles 
and s, t the are Mandelstam variables of the 2^2 process [177], with 



\{x,y,z) = {x-y-zf-Ayz. 



(46) 



Performing the integration yields : 



,9 \ ' ^2 ' ^2 



8Tf<V 



+ 



-l + 4r-4r^ + 2^(r^-r) 



m 



- + t-t' + 2— (3t'-t)+4— ,t 



<(-'^2>„2 



rrf 

9 ' 



16 Tf 



.r(l-r)+2-rMn^^ 



In 



1-v 
1 + v 



(47) 
(48) 



The heavy quark coefficient functions (47) and (48) differ from the results of Ref. [114] 
by a factor of t/Q^, which is due to different definitions of the structure functions F2 and 
Fl. In the asymptotic limit the cms velocity v is written as an expansion in 

(mVQ2) 



2'w? T ^ / 



In this limit the logarithms in (47) and (48) become 
'l-v 



(49) 



In 



l + v 



o 



m 



(50) 
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The expansion of (47) and (48) yields 



2,0 ' ) 2 ' 2 



ATf ■ { 8r(l - r) - 1 + [r^ + (1 - rf] In I ^ 



+ [r2 + (l-T)2]ln 



1 -r 



+4rF^<{ -lOr^ - T + 4(3t2 - r) 



In 



r 



1 -r 



(51) 



H 



(1) 



Q2 ^ 



16Tf ■ r(l - r) 



-16Tft' 



- 1 + 



in 



1 -T 



+ ln| 



m 



(52) 



In the limit m^/ — > the Wilson coefficients obtain the following structure: 



H. 



{1)1 m 



H 



2,9 
(1) 



r,^,^l oc 16T^-r(l-r). 



(53) 
(54) 



Often one considers Wilson coefficients in Mellin-space by performing the integral trans- 
formation 



dr T^-' H{t) . 



(55) 



The factor in front of the logarithm denotes the leading order splitting function [50,51]. 

Pli.r) = ST^r^ + (1 - r)^] , (56) 
which is a process-independent quantity. In Mellin-space it is given by 

+ N + 2 



''^ N{N + 1){N + 2) 



(57) 



The logarithm in (53) indicates the presence of a coUinear singularity, if — > 0. Both 
(53) and (54) were calculated in the so-called on mass-shell scheme for the outgoing quarks. 
While (54) is a scheme invariant quantity, redefinitions of the logarithmic contribution in 
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(53) for — > would yield different expressions, absorbing the divergent term and part 
of the constant contribution into the gluon distribution. The structure functions F^2L) 
are obtained by a convolution with the gluon density G{z, /i^): 



— , —5- ) G{z, ix^) , 



9 ' 9 



(58) 



where a = 1 + Aw? /Q"^. Inserting (58) into (25) yields the cross section 



dxdy xyQ^ 



(59) 
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4 The Heavy Quark Coefficient Functions in the Hmit 

> 

We compute the inclusive DIS heavy flavor production cross section in the asymptotic 
region 3> . There the heavy flavor Wilson coefficients Hi^j and Ljj factorize into 
massive operator matrix elements Aij and the massless coefficient functions Ci^fc, as has 
been shown in Ref. [118]. Here i^j^- are Wilson coefficients with a photon coupling to the 
heavy quark line, and Lj^ are those with the coupling to a light quark line, see Eqs. (80- 
84). All process dependent quantities enter only into the light flavor Wilson coefficients, 
whereas the complete mass dependence is contained in the massive operator matrix ele- 
ments Ajj, which are process independent. The power corrections which are proportional 
to {rn? IQ"^)^ , k > 1, can be disregarded in this limit. A quantitative comparison with the 
exact LO and NLO result in Refs. [110-114] and [115, 116] shows, that in the case of F^^ 
these power corrections can be neglected for Q'^/m'^ > 10, cf. [118]. 

Applying the light cone expansion to the partonic tensor corresponding to the inclusive 
Wilson coefficient Cf 'J'^'^^ yields the asymptotic factorization formula, [95] : 

E4--(/V,„,+ l.^)c|i2-(^,„,+ 1.^) +0(9!) , (60) 



where the quantifier (n/ -|- 1) denotes one heavy and n/ light flavors and /i is the factor- 
ization scale between the heavy and light contributions in Cj^a- The light flavor Wilson 
coefficients are denoted by Cjj and taken at + 1) flavors. In order to obtain the correct 
behaviour, it is necessary to include all radiative corrections containing heavy quark 
loops into the heavy quark coefficient functions. 

Let us consider ffist the unrenormalized OMEs Aij. They can be computed as projec- 
tions of truncated Green's functions. During the computation of the Green's functions 
trace terms emerge. However they do not contribute, since the local operators are trace- 
less. These terms are projected out from the beginning by contracting with 

Jn = A^,...A^^, (61) 

where is a light-like vector. The Green's functions used for the computation of the 
OMEs with external gluons are then given by, cf. [118], 

^'ip)Gt,u^''ip) = e^{p)MAlip) I Oq,,,...,, I Al{p))e%p) , (62) 
^'ip)GtQ,,.e%p) = e^{p)MAl{p) \ O,,,,...,, \ 4(p))Qe^(p) ■ (63) 

Here the the external gluon fields are denoted by with color index a and Lorcntz index 
^. The polarization vectors of the external gluons with momentum p are denoted by 
e^{p). The indices q, Q of the local operators O label the operator coupled to a light or a 
heavy quark. In the flavor non-singlet case the following Green's function contributes 

u{p,s)Gl'^'Xru{p,s) = Jiv(^.(p) I 0;^^^^...^J *^(p))q . (64) 

Here u{p, s) and u{p, s) denote the bi-spinors of the external massless quarks and anti- 
quarks, respectively, and the corresponding fields are ^ and ^. 
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Further OMEs are obtained from the following Green's functions with external quarks 
in the flavor singlet case 



u{p,s)G'^u{p,s) = Jm{%{p)\Oq,,,...,,\^^{p)) , (65) 

u{p,s)GIqu{p,s) = M%{p)\0,,^,...,,\^^{p))q. (66) 

S PS NS f \ 

The OMEs A^j ' iN,nf + l\ are obtained as expectation values of the following 
twist-2 operators, cf. [178,179], 

CgliWv./^iv = ^^"^S[V^7^iL'/,2 . . . ^MJvyV'] - trace terms , (67) 

OISmi,...,.. = l^''-'Sbf<T,,,D,, . . . D,,^^] - trace terms , (68) 

Oq-„,,...,„j, = i^-^S[V'7^,L)^2 . . . Df.^ip] - trace terms , (69) 

OLi,...,^. = 2^^-^SSp[F;^<,D,,...D,,_,F;;;^] -trace terms, (70) 



between corresponding on-shell partonic states {q | or (gr | and a^'^ — (i/2) [7^7*^ — 7''7'*]. 
Here S denotes the symmetrization operator of the Lorentz indices /ii, ...,/xjv, -D^ is the 
covariant derivative, ^ the quark field and Fj^^^ is the gluonic field-strength tensor. The 
flavor matrix of SU{nf) is denoted and Sp is the color trace. The operators are classi- 
fied as fiavor singlet (5") and non-singlet (A^^*) with respect to their symmetry properties 
under the fiavor group SU{nf). 

ForAg^ one obtains 

Aq, {myf,', e) = j^^i-9,.)Sa,{^.pr''G^i,. ■ (71) 
The OMEs obey the perturbative expansion 



^ i=i 

The singlet contribution has the following representation 

■^qq — ^qq + ^gg ('''3) 

in terms of the fiavor non-singlet (NS) and the pure-singlet (PS) contribution. 

Since any integral without scale vanishes in dimensional regularization due to the on- 
shell condition, all corresponding graphs but the 0{a^) term do not contribute. Due to 
this we have to consider only those matrix elements with at least one heavy mass line. 
For the singlet terms one has to distinguish the cases in which the operator is inserted 
on a light or a heavy quark line, respectively a vertex with a number of additional gluon 
lines, cf. Appendix B. These contributions are referred to as A^qQ and in the 
pure singlet case and by Aqg^g and Aqg in the gluonic case. In this thesis we compute 
contributions to A^^^q, A^^'q^, Aq^^, A^^q, and Agg, while the correction to Agg^g will be 
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given elsewhere [180]. Using the relation 



^S,PS,NS [ . ^S,PS,NS [ 



rrS PS / . 



Lr. 



S,PS,NS 

j,(2,L) 



Q2 rn" 



(74) 



one may split Eq. (60) into the following contributions, cf. [95], 



C;|2,L)(^/)+i^S2,L)(r^/ + l) = 



PS 



C5,(2,L)(n/) + Lg,(2,L)(n/ + l) 



Klqinf + 1) + <,%(n/ + 1) + >lg(n^ + 1) 
xnfCl%,L){nf + 1) + <,%(n/ + l)C,^f2,L)(% + 1) 

+ l)^/C'g,(2,L)(n/ + 1) , 

(75) 

^flfl,Q(^/ + l)"/C's,(2,L)(n/ + 1) 

+A.,Q(% + l)^i!(2,L)(% + l) 



+ 



^</p,q(% + 1) + + 1) nfCl%^L){nf + 1) . 



(76) 



<?2,L)(n/ + l) = Ag(n/ + 1) C,72,L)(n/ + l) + C,"?2,L)(ri/ + l) 



IPS 



yNS 



VPS 



+ K'q(^/ + 1) + + 1)JC'J?2,L)(^/ + 1) 

+ l)C'5,(2,L)(n/ + 1) , (77) 
^99,Q{r^f + l)C'g,(2,L)(n/ + 1) + A,g,Q(n/ + 1)C'J;^2,l)(^/ + 1) 
+AQ,{nf + 1) k^f2,L)(^/ + 1) + C',72,L)(^/ + 1)1 ■ (78) 



Here we used the notation 



/(%) 



f{nf) 



(79) 



The heavy flavor Wilson coefficients up to 0{al) can now be obtained by expanding in 
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Us, [95], 



LTi2,L){nf + I) = a 
+ a 



' ^?:f ("/ + 1) + 42f (»/ + i)c^;(iS('v + 1) 

i'ftDl"/ + 1) = 4[ + 1) + 4Sg(n/) n/C'<;|'2,i,(»V + 1) 

i(3),PS 1 
+%<^g,(2,L)(%)J , 

+ «s [ + 1) ^2 + <!Q(n/ + 1) ^^^(^(^^^^(n/ + 1) 

l(3) 



(80) 



(81) 



(82) 



+<,q(% + 1) C^lL)inf + 1) + >lgr'(% + 1) + 1)1 , (83) 



<Z,(2,L)^ 



(1) 



+ a. 



+ ai 



A^glinf + 1)6,+ Ag)(n, + 1) Cj^'^^^^inf + 1) 

+ + 1) ^i!i,L)(% + 1) + C^S2,L)inf + 1)] 

+ 1) ^2 + 41(n, + 1) ^(^(^"^(n, + 1) 

+ 4'!q(^/ + i)c'1;(U)(^/+i) 

+ + + 1) + + 1)} 

+ 4U^f + 1) ^fi,L)(^/ + 1) + + 1)] • (84) 

Here is given by ^2 = 1 for and by ^2 = for Fl . The relations (80-84) provide the 

basic scenario for the present work. To 3-loop order they connect the various massive 
OMEs, Aij, and the known massless Wilson coefficients, C^^^l)' f^'-'^]' ^'^ massive 
Wilson coefficients in the asymptotic region Q'^ ^ m^. Henceforth, we will calculate the 
massive operator matrix elements. 
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5 Renormalization 



The renormalization of the massive OMEs proceeds in four steps: i) mass renormahzation, 
ii) renormahzation of the strong couphng constant, in) renormahzation of the ultraviolet 
singularities of the composite local operators, and iv) subtraction of the colhnear singu- 
larities. To 3- loop order this formalism has been developed in Ref. [95]. In the following 
we summarize the main steps. 



5.1 Renormalization of the Mass 

In order to rcnormalize the heavy quark mass, wc will apply the on-shell renormalization 
scheme and define it as the pole mass. Hence the bare mass m is replaced by 



m 



1 + as i^-j^j 6mi + a; i^—j ^mgj + 0(af) . 



(85) 



The constants Snii and 5m2 are given by 



5mi 



e 

I ^ (18C,. - 22Ca + STpiuf + Nh)) + ^ (- yC^F + y 

\ /199 51 \ 

-UTj^iuf + TV,) j +Cf(^—- -C2 + 481n(2)C2 - I2C3) 

/ 605 ^ \ 
+Ca[-— + - 241n(2)C2 + 6C3) 



^/(f + 10C2)+VV,(^-14C2) 



/69 



+ 



(86) 
(87) 



(89) 



Here Cfe is the Riemann ^-function at integer values, cf. (C.IO), nf denotes the number of 
hght flavors and the number of heavy flavors, which we will set equal to N^ — 1 from 
now on. The pole contributions were given in Refs. [181,182], and the constant term was 
derived in Refs. [183,184], cf. also [185]. 

After carrying out mass renormalization up to O (a^) the OMEs have the following 
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structure 



+dl 



■ 5(2) /rn^ \ n 



/ dm ^ 



! 5(1) / m^ 



md 2''^^ /m? \ 5m\ (rn?\^rri^d'^ 
Here and in the following we unify all renormalization and factorization scales to one scale 



(90) 



5.2 Renormalization of the Coupling 

Considering only n/ light flavors and no heavy flavors yields the following relation between 
the bare coupling constant and the renormalized coupling a^^ 

^2 , 



The coefficients in Eq. (91) are given by [41,42,186,187] and [188,189], 

2 



-Mnf) 



with 
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^CA + CpjTFnf 



(91) 

(92) 
(93) 

(94) 
(95) 



If one considers also heavy flavor contributions it is important to take into account that the 
factorization condition (60) strictly requires the massless external particles to be on shell. 
This condition is violated by massive loop insertions to the gluon and ghost-propagators. 
Fortunately these corrections can be uniquely absorbed into the strong coupling constant, 
by applying the background held method, cf . [95] . The most direct way to do this, is to 
perform the renormalization first in the MOM-scheme and to transform afterwards into 
the MS-scheme. The following relations are obtained: 



= a 



+ 



MS 



5o,Qln(^)< 



(96) 
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or, 



-2Sa^,nnf + 1) [Sa^r ' (^/ + 1)J j + Oia^-^ , (97) 

vice versa, where a^^ = a^^(n/ + 1). These identities are vahd to all orders in e. In (96) 
the value for /3['^q, [95], are given by 

Pi,Q = Mnf)^-^(^lcA + Cp^TF , (98) 
/^S = -^TfCa + 15TfCf , (99) 
Pl% = -^TpCa-^TfCf-C2(^ItfCa + TfCf^ . (100) 

5.3 Operator Renormalization and Mass Factorization 

The ultraviolet singularities of the composite operators introduced in (67-70) is performed 
by introducing the corresponding Zjj-factors: 

O^...,,. = Zf.{f,^)dl^_^^,i^q,g. (102) 

Due to their identical quantum numbers mixing occurs among the different singlet oper- 
ators, (102). The anomalous dimensions of the operators are defined by 



dfjL 
d_ 

dfj, 

The NS and PS contributions are split in the following way 



tS = f,Z-'''{f,')^Zf^{t,'') . (104) 



A,, = A'^^ + A^^. (106) 
The anomalous dimensions can be expanded into a perturbative series in a^^: 

oo 

7r'''(«r,-/,A^) = E«?'7r''''''(-/'^)- (107) 

The renormalization is performed in two steps, see [95]. First only Uf light flavors are 
considered and then the renormalization scheme is extended to n/ light and one heavy 
quark flavor. 
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In the first case, one finds up to O(a^) up to 0{a^^ ) 



(0) 



1 /I 

72 



2e 



+a 



MS- 



1 ^1 (0) (0) (0) o (0) (0) 4 2 (0) 



+- 



1 /I 



(57. 



(i).,(o) 



7^ + 37l?^7j;^ + l^olir + ^ 



,Xi) 



At!"') 



+ 



(2) 

7tj 

3£ 



(108) 



As a second step the additional heavy quark is implemented. In order to limit the inves- 
tigation to the ultraviolet singularities for now, the external momentum is temporarily 
kept artificially off-shell. The corresponding Z-factors for the massive OMEs are then 
obtained by taking Eq. (108) at {uf + 1) flavors and applying the scheme transformation 

(97). Up to O (a^°^^^ one has 



7-1/ MOM ^ I 1 ,2 



Sij - a 



MOM I^J 



7 



(0) 



s 



+ a 



mOm2 



(1) 



Sa] 



'IOM^,(0)\ 

,1 7ij J 



<^<2°"7l°^) + ^(^/3o7g^ + 2(5a-r/3o7i?^ + ^/3i7g 



1 



, r MOM (0) (0) i (1) (0) 1 (0) (1) 

111 lij + q7ii 7ij + g7i; 7ij 



+ 



3' 

^(-3/^07, 



2.,{0) 



(109) 



The contributions oc Sa^^^ in (109) stem from finite mass effects and cancel singularities 
due to virtual processes at — > in real radiation. The OMEs are split into a purely 

light part^ij and a heavy flavor part^lj^-, which denotes any massive OME, we consider: 

n2 



4•(p^m^/x^ar",r^/ + l) = 4(-|-,ar^n/) 



+ A,(p^m^/.^ar",r^/ + l) 



(110) 



Here the light flavor part74jj depends on since the renormalization prescription for 
the strong coupling constant applies to the massive part only. The UV-renormalized 
expression is obtained by subtracting all terms that apply to the light flavors only : 

4(p^ m^ /i^ ar^, uf+i) = z-\a:-^, uf + 1, //^)ig(p^ m^ ar^, uf + 1) 

+Z^l(a^^ Uf + 1, /.2)4- (^, a?^, n/) 



*n ^2 '-^ 



;iii) 
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Here Zij can be expressed as a series in by 



^^' = '^.- + E«'^^'^'^- (112) 

k=\ 

Since in the dimensional regularization scheme all integrals without scale vanish, for the 
light flavor OMEs only the constant term 5jj remains in the limit 0. Expanding in 

yields the following UV-finite OME: 




In the on-shell limit — )■ collinear singularities emerge. These are absorbed into the 
parton distribution functions and only occur in massless parts of the OMEs. Thus the 
renormalized OMEs are obtained by 

The generic renormalization formula is given by 

A, = Zl^hu^l]. (115) 

In case of massless quarks only the F-factors would be given by, cf. e.g. [118], 

r,,- = Zr^ . (116) 

Due to the fact that collinear singularities emerge in massless subgraphs, only the F- 
factors have to be computed newly, cf. Ref. [95]. Finally the renormalized operator 
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matrix element reads : 



,MOM (^ig)'«(^) + ZT-^^in, + 1) - Z-'(^)(n,) 



+ 



p-l,(l) 



r-i,(2) 



(117) 



One notices that the order 0(£^)-terms of the leading order OMEs and the 0(£)-terms of 
the NLO OMEs are required to renormalize the 3-loop OMEs, cf. [118,120-123]. Trans- 
forming the coupling constant back into the MS-scheme and performing a series expansion 
in a^^ yields the final expression for the renormalized OME ^q^, cf. [95], 



A 



(3), MS 



-(0) 

Tqg 
48 



Iqg \ Iqq Igg 



-6^0 - 8/3o,q) + 8/3o' + 28/3o,q/3o + 24/3Jq| In^(^) + ^ 
-4/3o - 6/3o,g) + tJ? (7^;^ - 7« + (1 - nj^f'^ + 7«'^' + 7^^'^' " 2/? 



-2/3 



-(2) 
799 



£(2) J2) 



-/^ + ^(7j?-7S^-4/5o-4/5o,., 



- (0) ^ (0) y. 
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27j?-7S^-6/3o 
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-6/3o,q] -4/3o[2/3o + 3/3o,q]+7J?[-7J? + 6/9o + 4/3o,q]) lln(^) +agj( 



^(0) 
'99 



(1) 



+ ^0) 



7j? - 6/3o 



+ 8/30^ - 4/3o/3o,, 



-24/3o,Q + 



+ 



16 



'99 



(1),NS _ (1),NS _ ^(1),PS , 2/3^ 
I qq Iqq iqq ' ■^PJ- 



+2/3i,q) 



+ 



(-1) 



+ 7gg 



5m\ 



(0) 



- Sm^^^'j) + 5m? (7j? - 7^? - 2/3o - 4/3o,q) 



iX^ (0)^(0) ,„{3) 



(118) 



Similar expressions are obtained for the other OMEs, cf. [95]. Both from the unrenor- 
mahzed and the renormahzed OMEs one may extract the information which is new at 
0{al): 



the constant parts a]^ {N) of the unrenormahzed OMEs Aij 



the corresponding contribution to the 3-loop anomalous dimension ^^^^ [N] . 



Both quantities are computed in the following. 
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a(l) 

6 Calculation of the Operator Matrix Element 

As an example to illustrate the method used in principle we recalculate the LO operator 

matrix element Aq^ in the following. To the the lowest order in the coupling constant 
Qs = 5'g/(47r)^ it receives contributions from the two Feynman diagrams given in Figure 3. 
The evaluation of the Feynman amplitudes has been performed using the computer algebra 
system FORM [190]. One may split off an overall factor 



I [1 + (-1)"^] 



:ii9) 



from the OMEs, which we will do in the following. 




p, a,ii 




q-p 



q-p 



(1) 



Figure 3: Feynman diagrams contributing to Aq'^ 



p, h,p 



In the unpolarized case and for pure photon exchange only even moments contribute. 
This is a consequence of the current crossing relations, cf. e.g. [170,171,191] and the local 
light cone expansion [43-46]. One obtains 



A 



Qg,{a) 



A 



QgXb) 



2{N^ + 3N + 2)+ e{N^ + N + 2) 
N{N + 1){N + 2) 

32asTFSs { — I TTTT — ^ exp ( 2J 



/iV {2 + e)e ^\^^n2J J {N + l){N + 2) 
Here denotes the bare coupling constant and is given by 



(120) 
■ (121) 



Se = exp (7e - ln(47r)) 



2i 



(122) 



with 



r(i-e/2) 

(4^)(^/2) 



\l=2 



(123) 



31 



cf. Appendix F. One obtains 



1 

as 



Qg xQfl 



— —SfrTp I — 

/i2 



2\ e/2i 

\ ' 1 



- exp 

s 



m 



8{N^ + N + 2) 
N{N + 1){N + 2) 

^y/V 1 C2 CS 2\ 8(iV2 + iV + 2) 



1=2 



I \2 



e 8 24: J N{N + 1){N + 2) 



(124) 
(125) 

(126) 



The A'"-dependent term is identified as the MeUin transform of the LO sphtting function 
Pqg\N), (56). In 2;-space the unrenormahzed leading order OME reads 



^r)^ — Sp { — 



m 



2\ e/2 



with 



.(1) 



cf. Ref. [118]. Expanding up to 0{e'^) yields 



Qg 24" 



(127) 

(128) 
(129) 

(130) 



1(1) 
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1 1 , fm' 

m — r- 

e 2 Va* 



1 , f m 
— In — r 



, , 1 



16 



24 



+ 1 



;i3i) 



^Qs possesses a single pole, which is removed by operator renormalization later. The 
other contributions up to O(e^) all contribute to the renormalized OME at 3-loop order. 

Let us now discuss the results of Section 3, Eqs. (53,54) in the context of the asymptotic 
heavy flavor Wilson coefficient Eq. (118) to 0{as)- In case of the structure function 
F^'^ {x,Q'^), ifq^ is given by the massless result since 5^ — 0. This applies to all 



schemes, cf. Ref. [98]. Eq. (54) is actually the same in the heavy and hght quark case 
after performing the limit m^/Q^ — > 0. Comparing H^^, Eq. (53), calculated using the 
on-mass-shcll scheme in Section 3 with the massless result obtained in the MS-scheme, 
cf. [98], one obtains the same result. This is in accordance with 



^Qg 







(132) 



Eq. (132) determining the finite part of in Eq. (118). Corrections only occur in O(a^) 
and higher. This gives a first illustration of the formalism. The corresponding results at 
NLO are given in the literature, cf. Refs. [118,120-123]. 
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7 Calculation of the 0{alTpNfCjr(^A)) Contributions 
to the Massive Operator Matrix Elements 



In the following we describe the computation of the contributions to the massive 3-loop 
OMEs of 0{TpnfCF,A) being performed in this thesis. They concern Aqg, Aq^, ^qq,Qi 
A^^Q in the unpolarized case and A^^q^ for transversity. The results for Aqg^q will be 
presented in [180]. The unrenormalized OMEs are obtained by applying the following 
projectors to the corresponding truncated Green's functions. For external gluons the 
projector Pg, cf. [118], 



Gf, 



ab 



^ab 



-.fJ,L> 



m-iD-2 



(133) 



defines the corresponding scalar contributions. In Eq. (133) the summation over /x, u in- 
cludes unphysical gluon-states, which have to be compensated by adding ghost-diagrams. 
For the external quark contributions the Green's function is projected by 



Gf 



Another projector is needed for the non-singlet transversity OME. It reads 

^ _,^^^i;L_^(A |Tr[^j/p^ G^j,TR,NS, 



(134) 



- A.pTrb'^ GfJ^^'""'] + iA.pTr[ a^'^ jf GfJ^^'""'] ] , (135) 



cf. Ref. [179]. 



7.1 Contributing diagrams 

The Feynman diagrams contributing to the respective Green's functions have been gen- 
erated by a code [95, 192] allowing for local operator insertions in QCD-diagrams based 
on QGRAF, [193]. The diagrams Ai to i?^^, see Figure 4, can be generated by including 
one massless fermion loop into the gluon propagators of the corresponding massive 2- 
loop diagrams [118,120]. The color factors have been evaluated using the FORM-package 
color, [194]. For the diagrams that could be obtained from the 2-loop diagrams by in- 
serting a light quark loop, the color factor differs by TpUf. This massless insertion can 
easily be integrated out and is given by the following expression for the extended gluon 
propagator including the quark loop 



- 4^2+e/2 B (,2 + 2' 2 + 2 j r [-^) (_^2)2-./2 ' (136) 

with gs defined in D = 4 -|- £ dimensions. 

The main differences to the Feynman rule for the gluon propagator, see Appendix B, 
are slightly more complicated numerator striictures, and the e dependence in the power of 
the denominator. While the new numerator structure in most cases just increased the size 
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R R' S 

Figure 4: 2-loop diagrams, that were used to generate 3-loop diagrams by the insertion 
Eq. (136) 
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K2 L2 



Figure 6: 3-loop diagrams, that are not generated from massive 2-loop graphs 




A B 



Figure 7: Diagrams contributing to Aq^ 




A B 
Figure 8: Diagrams contributing to A^^q 




ABC 



Figure 9: Diagrams contributing to A^^^^ 
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of the computation, the occurrence of real exponents in many cases made it necessary to 
recompute the decorated 2-loop diagrams completely, since new structures occurred. In 
Figure 4 we show the 2-loop diagrams of Refs. [118, 120] into which the massless fermion- 
bubble is inserted. This usually concerns more than one line. In some cases the insertions 
on different gluonic lines lead to different integrals, that cannot be mapped onto each 
other via a symmetry relation. We labeled this accordingly, e.g. Jib, Jib, etc. Figures 4 
and 5 include the corresponding ghost diagrams which contribute, since the calculation is 
performed using a R^-gaugc. In Figure 6 wc show the other topologies which contribute to 
Aqg at 0{TpnfCF,A), but cannot be generated by loop-insertions into known topologies. 
In this thesis the diagrams A2 — H2 are computed. In the final result we will show the 
contributions to I2 — L2 as well, which are calculated in Ref. [180].^ Finally Figures 7- 
9 show the corresponding diagrams contributing to the OMEs A^^, ■^qq,Q^ ^gg^Q' 
A^^'q^. In Tables 1-4 we hst the combinatorial multiphcities through which the different 
diagrams contribute. These were determined from the foregoing computation of the fixed 
moments for the respective OMEs in [95, 192]. 



Table 1: Multiplicities of the individual diagrams contributing to A^q 



Diagram 


Multiplicity 


Diagram 


Multiplicity 


Diagram 


Multiphcity 


Ai 


2 




4 


A2 


4 


Bi 


4 


Lib 


4 


B2 


4 


Ci 


2 


Ml 


4 


C2 


2 


^1 


4 




4 


D2 


2 


E, 


4 


Nib 


4 


E2 


4 


Fi 


2 


Oi 


4 


F2 


8 


Gi 


4 


Pla 


2 


G2 


4 


Hi 


2 


Plb 


2 


H2 


4 


h 


4 


Sla 


2 


h 


4 


Jla 


2 


Sib 


2 


J2 


8 


Jib 


4 


Tla 


1 


K2 


4 




2 


Tu 


1 


L2 


8 


Ki, 


4 


Tib 


2 







Table 2: Multiplicities of the individual diagrams contributing to A, 



Diagram 


Multiplicity 


Diagram 


Multiplicity 


A 


4 


B 


4 



^Results on the scalar integrals for diagrams I2 — L2 are given in Ref. [195] also. 
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Table 3: Multiplicities of the individual diagrams contributing to 



Diagram 


Multiplicity 


Diagram 


Multiplicity 


A 


4 


B 


4 



Table 4: Multiplicities of the individual diagrams contributing to Aqq'^ ' and A, 



Diagram 


Multiplicity 


Diagram 


Multiplicity 


Diagram 


Multiplicity 


A 


2 


B 


4 


C 


2 



7.2 Evaluation of the Feynman diagrams 

The two-loop massive operator matrix elements were calculated in [120,121,196]. Proto- 
types of graphs were computed in [197]. Many diagrams contributing to the T|.n/-term of 
the three loop operator matrix elements can be evaluated by considering the correspond- 
ing 2-loop diagrams, see Figure 4, and replacing one gluon propagator by the extended 
gluon propagator containing the one loop self energy, cf. Figure 10 for an example. 

Figure 10: Replacement of the gluon propagator by an extension including 1-loop massless 
fermion contributions 

The first part of the evaluation of the individual Feynman diagrams has been per- 
formed in a procedural way using the algebraic manipulation program FORM, [190]. The 
occurring fermionic traces of 7-matrices have been calculated using the built-in functions, 
before the momenta were integrated. To integrate a respective momentum ki all de- 
nominators containing this variable were combined using Feynman parametrization, see 
Appendix (C.12). Here the order in which the momenta are integrated is not arbitrary. 
In many cases specific choices allow to avoid infinite sums emerging from generalized 
hypergeometric functions. As a next step the ^-distributions are integrated out and the 
momenta are shifted in order to symmetrize the D-dimensional momentum integral. Due 
to this factors of the form (A./cj -|- A.l)^ occur, where I denotes contributions from other 
momenta in the form 

I = Pl({Xk})kn + --- + Pn({Xn})K, (137) 

with polynomials in Feynman parameters Pi. These terms have to be expanded using the 
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binomial theorem prior momentum integration 



{A.k, + A.lf = J2(^■)(^■^f~'(^■^^y■ (138) 

All but the first terms of this series can be dropped, as the symmetric momentum integrals 
over odd powers of ki vanish and the rules for £)-dimensional integration, cf. Appendix 
C, lead to terms containing the contraction A. A = for higher powers of A.ki. In this 
computation a maximum of the first three terms had to be considered. The symmetric D- 
dimensional integral was evaluated by applying the rules in Appendix C. These steps were 
repeated for all internal momenta, which yields Feynman parameter integrals of various 
complexity. The application of simple algebraic transformations leads to representations 
of the following form : 

XP3(xi,-- - ,Xn,j) , (139) 



or 



Here the Pj are respective polynomials in the Feynman parameters, rii G Z, and P{{xk}) 



denotes a product in the Feynman parameters {x^}. The physical sum X^jL'o^ emerges if 
the operator insertion is located at a vertex, cf. Appendix B. In general also diagrams 
with an operator insertion being represented by a physical double sum Yl!j=Q Yld=j+i 
have to be computed. For the contributions considered in this paper, these operator 
insertions where always located on an external vertex. Due to this one of the sums could 
be performed at the momentum level: 



Y,{A.k 



- A.p)^-'-2(-A.p + A.q)'-^-\A.q) 



0<j<l 



J2{A.k 



0<j<l 



^ N-2 

— J2 {(A. A: - A.pf-^-' [{A.qy - {A.q - A.pY] } , (141) 
■P j=0 

A.pf-'-\A.ky-^-\A.qy 

^ N-2 

= E { [(^-^y - (^-^ - ^-py] i^-^f-"-'} ■ (142) 

■P j=0 

In some cases a simpler structure is obtained, since variable transformations that map 
that map P{{xk}) to zero are applicable, cf. Appendix E. If possible the physical sums 
"12^=0 were evaluated at the Feynman parameter level. All Feynman parameters that 
could be integrated through substitutions or in terms of Beta-functions were integrated 
at this point. In case the relations given in Appendix E were not applicable, and no 
further Feynman parameter could be integrated out directly, the binomial theorem was 
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used thereby introducing additional finite sums. These steps have been repeated until 
all Feynman parameters, that do not contribute to the denominator are integrated out. 
The denominator in (139-140) and remaining Beta-like factors in the Feynman parame- 
ters {xk} were collected in terms of integral representations of generalized hypergeometric 
functions pFg, [136,137]. Thus a representation in sums over generalized hypergeometric- 
and F-functions has been obtained. The hypergeometric structure of the Feynman dia- 
gram is determined through the mass distribution of the Feynman graph and is widely 
independent of the operator insertion. We have implemented the relations being listed in 
Appendix F.2 into a FORM-algorithm, which allowed to perform infinite sums stemming 
from generalized hypergeometric functions pFg on this level in many cases. 

The results are now given in terms of finite and infinite sums over F-functions depend- 
ing on the MeUin- variable N and the dimensional parameter e. The respective expressions 
were simplified and the integrals were expanded into a Laurent series in e using the com- 
puter algebra system MAPLE [198]. In some cases remaining sums could be evaluated by 
applying known results from the literature, cf. Refs. [120,121,199,200]. More complex 
sums were evaluated using the MATHEMATICA-based program Sigma, [142,143]. Some 
examples for typical sums that were obtained during this work are given in Appendix G. 

For the cases of the diagrams Ai-Di, Gi, Jia-Kn,, the insertion of the massless quark 
loop changed the structure of the occurring Feynman parameter integrals moderately. 
The evaluation of these diagrams could be performed in terms of a modification of the 
computation of the 2- loop results [120]. All other diagrams had to be computed newly. 

According to the complexity of the topology, between one and three sums had to be 
performed after the e-expansion. In the following we demonstrate details of the calculation 
considering some examples. 



A typical case for a representation without any remaining sums prior the Laurent expan- 
sion is given by diagram E2. After applying the projector (133) the color and 7-algebra 
is performed. The momentum assignment is as follows: pi — ki, p2 — ki — k^, ps — k2, 
P4 — k2 — p, P5 — ki — p, P6 — ks, p7 — ks — k2, see Figure 11. One obtains the following 
representation : 



7.3 Diagram E2 




Figure 11: The Momentum flow for diagram E2 




N-2 



/(A^,A:i^,A:2^,A;3^,p^,m,£)(A./c3)-^(A.A:3 + A.A:2)^-^'-^ 



(143) 
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The factors {A.ksY {A.k^ + ^.^2)^'^'"^ in (143) stem from the operator insertion on a 
vertex. Applying the Feynman-parametrization (C.12) to combine all the denominators 
which contain the momentum and integrating out the 5-distribution yields 



^ N-2 



[\ f f d^k2 [ dPkz^ 

io V (27r)^y (27r)^7 (27r)^ 



/(A^, /ci^, fc2^, fc3^,p^, m, £)(A./c3)^'(A.fc3 + A.fc2)^-^'-^ 
-D(A;i, A;2,P, m){m? — k^ — 2xoks.k2 — a^o^i)^ 

Here D{ki,k2,p,m) denotes all the factors in the denominator of (143) which do not 
depend on the momentum k^. In order to symmetrize the momentum integral we shift 
the momentum /ca — > Z3 — xok2- This yields 

1 ^ f\ f d^'ki f d^k2 f d^k 



^ J{^u.,kif,,k2f,,h^,p^,,m,XQ,e){/\.h - xoA.k2yiA.l3 + (1 -a;o)A.fc2)^ ^ ^ , . 
D{ki, k2,p, ra){m'^ — k^ — Xq{1 — Xo)kiy 

Now we perform a Wick rotation to obtain an Euclidean momentum integral, which can 
be evaluated as described in Appendix C. The results depend strongly on the momentum 
structure in the numerator, and become rather lengthy in many cases. During this com- 
putation many intermediary results amount to several Mbytes of data. For this reason 
we consider a numerator of the form 

f{^n,kin,l2n,Pix,m,N,e,xi,X2) = m^A.p^ (146) 

to demonstrate further steps of the computation. In this case only the first term in 
both powers of the respective binomial expansions of the factors (A./3 — a;oA.A;2)-' (A./3 + 
(1 — a;o) A.A;2)^^-'^^ contribute. Wc thereby obtain the following intermediate result after 
integrating out the momentum 

r _ T^n.C.a' ^ m\A.pmD/2,2-D/2) ^^ 

Ie, - T^UfCAQs (A.p)^(47r)^/2 (47r)^/2r(-D/2) j^} "^"^^ ""'^ 

^ f {^■k2r-' .,47. 

J {2n)Dj {2n)DDik,,k2,p,m)im^-Xoil-Xo)p2?-^/^ ' 
The physical sum can now be evaluated in terms of a geometric sum 

7V-2 

Y,{-iy4i^-^of-'-' = (-l)^x^^ + (l-a;o)^-^ (148) 

j=0 

Thus the following representation is obtained 

2^61 m''{A.pfB{D/2,2-D/2) 



(A-p)^ (47r)^/2r(-L>/2) 

X I' dxo x-J^-\l - x^r^-^ [(-!)-<- + (1 - xof-^] I^J^ 

.1 (^•^-)"'^ (149) 

D{k2){-ki){-{k,-py){-{k,-k2r){my{xo{i-xo))-kiy-^/^ • ^ ' 
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In the term, which contains the factor {1 — xq)^ ^, the variable shift xq — )■ {1 — xq) is 
performed. The same steps as above are apphed to integrate over ki. This yields 

, _ r ,6 1 m^BiD/2,2-D/2)B{D/2,3-D/2) 

X dxo x^^^'/'-'il - xo)^/'-' (1 + (-1)^) (150) 
Jo 

D{k^){-kl){-{k^ -pf){-{x2{l - X2)k2 - X2X3py){my{xo{l - Xo)) - klf-^/^ ■ 

The transformation x^ ^ x^{l — X2) maps the integral restricted by the ^-function to the 
domain [0, 1]. Finally, integrating over the momentum A;2, one obtains 



X 



X 



X 



da;oXo^+^/^-^(l - x,ri^-^ [1 + (-1)^] 

dx2 [' dxsx^/^-\l - X2TI''-'' [' dx^xril - x^y-""/^ 
Jo Jo Jo 

/ dx5 dxa 6{1 - X5 - xe)xl~^^'^{l - X5 - X6f~^^^{x6 + 
Jo Jo 

Xo(l- Xq) "I 



m^{l - Xr){l - X5- Xq) 



(151) 



with D = 4 + s. One notices, that after applying the transformation — )■ (1 — X5)xq the 
denominator factorizes into Beta-function like structures. No hypergeometric functions 
are required to represent this Feynman-integral. The integral over X3 can be directly 
performed, yielding 

X /' dx„i^°+''(l - Xo)'-" [1 + (-1)"] 
Jo 

X dX2 X^^^-\l - X2f^'-' t dX7 X,{1 - X7)-'+^ 

Jo Jo 

X [\x5 [' dxe xl-^'^^l - x,f-\l - xef-^ 
Jo Jo 

X I [(1 - 0:5)^-^0;^^ - (X5 + xe- x,xef^'] I . (152) 

In the term, which contains the factor (^5 + Xq — x^Xq)^"^, the transformations 

0:5 — >■ 2:5X6 (153) 
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and 



3:5(1 - Xq) 
1 - X5XQ 



(154) 



cf. Appendix E, are applied to obtain a representation which can be represented in terms 
of Beta-functions : 



Ie2 = TpUfCAg. 



Jo 

X rdx2 X?/'-'(l-X2)^/'-' rdXTX,{l-Xj)-'+'' 

Jo Jo 

X [ dx5 [ dxe (1 - x^y-'^/^x^y 
Jo Jo 



X- 



N-1 



D-\^D-b 



. (155) 



Performing the Feynman parameter integrals, and applying the conventions in Ap- 
pendix A, yields the following representation in terms of F-functions: 



Ie2 — TpUfCACi 
xF 



m^V^^' 1 1 



^2 J (47r)3/2e N - 1 

6/2, e/2, 1 + e/2, e~2,2- 3/2e, 3-e,2 + N -e 
l + e,5 + N -2e,N + e/2 



x< F 



N-s/2 
1 - e/2 



F (AT) ^ X 



N 



(156) 



In (156) no more sums remain to be evaluated. Thus, the final result is obtained just by 
expanding in e. 

7.4 Diagram Lia 



P3 




P5 ^ L9^ P4 

Figure 12: Momentum flow for diagram Lia 

For this diagram the following representation is obtained after performing the color and 
Dirac-algebra and integrating out the massless quark insertion : 
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-£/2,2 + £/2,2 + £/2 

4 + £ 



(4 



1 1 r d^ki r d^k2 

:7r)^/2 (A-p)^ J (2^7 (2^ 



A'-l 



-A;i)(m2 - A;2)2((-A;2 -p)2)2-e/2(rn2 _ (yti -p)2)(m2 - (A;i - A;2)2) ' 



(157) 



Here the momentum flow is Pi = ki, p2 = ki — k2, ps = k2, P4 = k2 — p, P5 = ki — p, 
cf. Figure 12. The factor {A.ki)^~^ in (157) stems from the operator insertion on a quark 
hne. A Feynman-parametrization is apphed to combine all the denominators, which 
contain the momentum k2- The emerging 5-distribution is integrated out. This yields 



TlufCAQl r 



-£/2, 2 + £/2, 2 + £/2, 3 - e/2 

4 + £,l -£/2 



(47r)^/2 (A-p) 



N 



xy dxi j dx20{l - xi - X2)x2^^^ J j^:^ J 



d^ko 



X 



, (2n)D J (27r)^ 

/(A^, fci^, /C2^, p^, m, AT, £)(A./ci)^-^ 

D{ki,p, m){—kl + xim? + 2xiki.k2 + 2x2k2.p — xikf)^'^/'^ 



:i58) 



Here D{ki,p,m) denotes all the factors in the denominator of (157) that do not depend 
on the momentum k2- In order to symmetrize the momentum integration we shift the 
momentum A;2 — >■ ^2 + Xiki + X2P- This yields 



-£/2, 2 + £/2, 2 + £/2, 3 - £/2 

4 + £,l -£/2 



X 



X 



(47r)^/2 (A-p)^ 



/(A^, /ci^, k2f„ p^, m, iV, £)(A.A;] 



\N-1 



D{ki,p, m)(— Z| + Xim? — Xi{l — Xi)kl + 2xiX2ki.pY~^l'^ 



(159) 



Considering only the contribution 



f{Au_,ki^,l2n,Pn,m,N,e,xi,X2) = (A.p)m^ 



(160) 
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for now and performing the I2 integral gives 



-6/2,2 + e/2,2 + e/2,l-e 

1 -£/2,4 + £ 



d^ki (A.y;) ///^ (A./,-i)'^'-^ 



X 



(27r)-^ D{ki,p, m)(xim^ — Xi{l — Xi)kl + 2xi{l — Xi)x2ki.py~ 



(161) 



-e/2,2 + £/2,2 + £/2,l-£ 

1 -£/2,4 + £ 



X 



X 



J dXi J dX2 Xl-\1 - X^y/^X^'^^ 

J {27r)^{m''-kl){m''-{ki-py){my{l-xi)-kl + 2x2ki.py-'^' ' 



Combining the remaining denominators inside the momentum integral and symmetrizing 
by performing the momentum shift ki ^ li + {x^ + X2Xq)p yields 



hia — 



X 



r 



-£/2,2 + £/2,2 + £/2,4-£ 
1 -£/2,4 + £ 



(A-p) 



N 



I dXi f dX2x\-\l - X{)'/^X~''^ 

Jo Jo 

X dx5 dxeXQ^6{l - X5 - xe){l - X5 - Xq) 
Jo Jo 

r d% {A.p) n,\A.I, + {.r-, + .r2.r,)A.pf-' 

""J (27r)^ (^2(1 1/(1 _^^))_^2)4-e • 



(163) 



When expanding the factor (A./i + (2:5 + X2Xq)A.p)^~^ as described above, only the term 
A.p^~^{xts + X2Xq)^~^ contributes. Applying the transformation x^ 0:5(1 — Xq) and 
performing the momentum integration yields 



-r 



£/2,2 + £/2,2 + e/2,2-3/2£ 

1 -£/2,4 + £ 



X / dXi / dX2Xl~^{l — X\fl'^X2^ 

Jo Jo 

X dx5 dxQ{l - X5)xq%1 - xef 
Jo Jo 

(A.p)m^ ((^5(1 - xe) + X2X6) A.p)^-^ 

^ (m2(l-X6(l-/(l-Xi)))2-3/2e 



(A-p) 



N 



(164) 



The factor {1 — X5) is expanded. For reasons of brevity we will only consider the first 
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term. The integral over can be performed directly. Shifting X2 ^ {I — X2) gives 

1 



-£/2, 2 + e/2, 2 + e/2, 2 - 3/2£ 

1 -£/2,4 + £ 



X 



X 



/ dxi / dx2x\ ^{1 — xiY ^^2^^^ 
Jo Jo 

/■I _ 1 

/ dxQXQ^{l- 
Jo 



{1 — Xi + XqXi^ 3/2£^^2~j3/2e 



x-[{l-xeX2f -{{1-X2)xef] 



(165) 



The factor (1 — xeX2)^ cannot simplified further and a binomial expansion is applied 



-£/2,2 + £/2,2 + £/2,2-3/2£ 

1 -£/2,4 + £ 



X /" dxx ! dx2x\-^{l-xxf-'x2'^'^ 
Jo Jo 

X y dxex^'il- xe)-^ 



Xi + XqXi)'^ 3/2£^j^2^3/2£ 



X- 



N 



N 



J2{-XeX2y^-{l-X2fx^ 

J 1=0 



(166) 



Finally the variable shift ^ 1 — a;6 is applied and the integration is performed. We 
obtain the following result: 



(47r)3/2s N ^ 



1 ^ 

-T 



-1) 



xT 



X3i^2 



-£/2, £, 2 + £/2, 2 + £/2, 2 - 3/2£, 3 - £, 1 + ji - £, 1 + ji - £/2 
1 - £/2, 4 + £, 2 + ii - £/2, 3 + ji - £ 

2-3/2£,2,£ 



3 + ji - £, 3 



(167) 



As the generalized hypergeometric function in (167) contains both one positive integer 
in the upper and lower indices, it can be summed by using the relations described in 
Appendix F.2. 

7.5 Diagram I\ 

As a last example we consider diagram 1\. It has the following integral representation: 



rf^fci /• d^k2 /■ d^k^ 



X 



(27r)^y (27r)^ 

/(A^, fci^, fc3^,p^, m, £) [Ci(A.A:i)'--^-i + ^^(A.fci + A.A:^ + A.p)'-^-^] 
(m2 - A;2)(m2 - (/ti - p)2)(-A;i)2(m2 - {k^ + A;2)2)(-A;i)(-(A;2 - k^f) 



(168) 



(169) 
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P4 

Figure 13: Momentum flow for diagram Ii 



The momentum flow is pi = ki, p2 = /c2, Ps = ki + /c2, Pa = ki — p, see Figure 13. One 
of the sums is performed according to Eqs. (141,142). The integration of the momenta is 
carried out analogously to earlier examples. Again we limit this demonstration to a small 
part of the numerator by considering only 



/(A/,, ki^, /C2^, k3^,p^, m, e) = {A.pf 



(170) 



Furthermore we hmit ourselves to terms which contain the color factor Ca- The following 
representation is obtained: 



I I, = T^UfCAa. 



m 



2\ 3/2£ 



S \ 9 



2 + 8/2,2 + 6/2,-3/26 

1 -£/2,4 + £ 

N-2 



N 



X 



(47r)3/2D 

/ dxo / dxi / dxi / dxa'^^x^l — xi)~^~^{l — xiX3)^~^~-'{l — xqXi) 

Jo Jo Jo Jo ,_n 



3/2£ 



j=0 



X [(1 - a;iX3(l - xo)y - {xiXail - xo)y] 



(171) 



Here the two terms in the last factor stem from the structure in Eqs. (141-142). Perform- 
ing the sum "^fS^ in the second term by a geometric sum yields 



nrifCAai ( — 
A* 



2\ 3/2£ 



2 + e/2,2 + e/2,-3/2e 

1 -£/2,4 + £ 

N-2 



N 



(47r)3/20" 

X / dxQ I dxi / dxi / dx3 Xq{1 - xi)~^~^{l - xpXi)^^'^^ 
Jo Jo Jo Jo 

x{-{l-X,X3f-^-^{xiX3{l-Xo)y 



+ 



X0X1X3 



[(1 - X,X3{1 - Xo)f-' - (1 - X.X^f-'] I . 



(172) 



Due to the stronger nesting of the Feynman parameters, no simplifications by using the 
variable transformations in Appendix E are possible. The terms {1 — x^x^ll — Xo))^~^ and 
(1 — xiXs)^"^ are expressed by binomial series. The remaining integral is then evaluated 
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in terms of the hyper geometric function 3F2 yielding the following representation: 



7/1 = T^UfCAttl 

N-2 N-j-2 

-E E 

j=0 ji=0 

X3F2 



2 \ 3/2e 



1 



N-j-2 
h 



(47r)3/2^ 

r 



-£/2, -3/2e, 2 + s/2, 2 + e/2, -£ 

4 + £,l -£/2 
l+£,l+j-£,A^-Jl-l 

2+J,A^-ii-£ 



+E 

ji=0 

N-l 

-E 

ji=0 



iV- 1 
N -I 



-3/2e,l + e,N -ji 
2 + j,N-ji-e 

l-e,e,N 
N-j,-e 



-3/2e,e,N-j, 
l,N-j,-e 



e,N-j^-l 



-3/2e,e 
N-ji-e 



;i73) 



The Gaufi-function 2-^1 can be mapped to a product of F-functions by applying Gaufi' 
theorem (F.14). Examples for typical sums, emerging after the series expansion in e are 
given in Appendix G. 
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8 The Massive 3— Loop Operator Matrix Elements 
and 3— Loop Anomalous Dimensions 



In the following wc summarize the main results of the calculation. We first obtain the 

unrenormalized massive operator matrix elements Aij for the various channels. Their 
analytic structure is known, cf. [95]. At 3 — loop order two new quantities can be obtained: 
i) the constant part of the OMEs A\f , a\f; ii) the three loop anomalous dimension 
corresponding to the process considered, 7^^^ The latter term is part of the 1/s pole 

term of A^. We consider the OMEs Aq„ A^q% AJS^ A^^q and A^J;^"" in the following 
and determine both quantities. They are represented in terms of the algebraic basis [131] 
of the harmonic sums, which leads to very essential structural simplifications and avoids 
redundancies present in other representations for e.g. the 3-loop anomalous dimensions 
[93, 94] and the massless Wilson coefficients if compared to the representations in [109, 
132, 133,201,202]. The contributions to the operator matrix elements calculated below 
constitute the first contributions for general values of N at 3-loop order and, due to 
Eqs. (80-84), to the heavy flavor Wilson coefficients in the asymptotic region ^ m^. 
The corresponding renormalized OMEs can be obtained following Ref. [95], see also (118), 
cf. Section 5. 

8.1 The gluonic contribution 

The unrenormahzed OME Aqg has the following structure, cf . [95] , 



;(3) 



-(0) 

6^3 



+28^o,qA + 24/3Jq + 7^°; [i'J + + 14^o,q] j + ^ (^^'J [H'J - 2^f^ 



-8^0 - 10/9o,Q 



^ 'qg 



4ff 



i(2) 
'^19 



-(0)/- 

16 



^(0) 
199 



{27. 



(0) 
99 



+4/3o,Q/3o + 24/31 



6m[ 



27« + 36m[-'m^'J + 25mf)7j? 



qg 

+Sm^^H'J [7^? - 7j? + 2^0 + 4^o,q] - Sm^^'J 



(3) 

Qg 



(174) 
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Here and in the following we always project onto the color factors TpUfCp and TpUfCA- 
The computation performed in this thesis yields : 



A 



(3) 
Q9 



N{N + 1){N + 2) 



'(N^ + N + 2)( 512(7V2 + 7V + l) 256 \ 

+^ + ^) [- 9(7V - 1)N{N + 1)(7V + 2) + T^') 



{N^ + N + 2) 



12^ 



2 , 256 128 
——0—2 -\ ;;~>-'2 



128(5Ar4 + 20N^ + 59N^ + 7QN + 20)Si 



32Qi(Ar) 



27{N +1){N + 2) 



27{N - l)7V2(iV + l)2(iV + 2)2 



1 

+ - 
e 



256 832 



64(A^2 + TV + 1) 



128 128 320 

^ "9~ ^ "9" ~ 3(A^ - 1)A^(A^ + 1)(A^ + 2) 



C2 



128(5A^^ + UN^ + 53N^ + 82N + 20) 

27(A^ + l)(A^ + 2) " 
32(10A^^ + A9N^ + 154A^2 ^ ^97^ ^ 53) 



32Q2(A^) 



27(A^- l)A^(A^ + l)(A^ + 2) 



27(A^ + l)(A^ + 2) 
16Q4(A^) 



32Q3(iV) 



81(A^ + l)2(Ar + 2)' 



81(A^ - 1)N^{N + l)3(Ar + 2)3 
1 



+- 



{N' + N + 2) 



64g5(iV) 



256 



9{N -1)N^(N + iy{N + 2) 



Si 



+ 



704 



128.,, 2 I .r I I 128(5iV3 + lliV^ + 28iV + 12)5i 

+ i (Ar^ + 7V + 2)(^-^S? 
64(5Ar3 + 20iV2 + 37A^ + 12) 



16{N-2)Qe{N) 



27{N - l)m{N + l)3(Ar + 2y 



32 



C2^1 



128 



27 3 ' 9 

8{N^ + N + 2)Q5{N) 



S2S1 + 



<2 + 



27N 

32Qr{N) 



3{N - 1)N^{N + l)2(iV + 2) ^ 9{N - 1)N^{N + 1)^{N + 2) 
64(10iV^ + 86A^3 + 483A^2 + 341iV + 114) 



^2 



81Ar(iV + 1) 



-Si 



81{N - 1)N^{N + l)4(iV + 2)3 



.(3) 



(175) 
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The contributions due to the individual diagrams are given in (D.2-D.200). The constant 

,(3) 



part QqI reads : 



16(iV2 + N + 2) 



^'q', = rifT^CA I 27iV(iV + l)(iV + 2) r^^-^'^'^ ' ^^^^'^'^ ' ^^^-^-^ + ^2^2,-2 - 6^3,1 
-1085_2,i5i + 42S2,i5i - 65_4 + 90S^sSi + llS^s^i + 120^4 + 185_2^2 + 54S_2Sl 
+33^2^1 + 15^2' + 2^1 + 185_2C2 + 9-52(2 + 9-51(2 - 42^1(3 

+ UN^ + 53N^ + 82iV + 20 r. 



+32 



27iV(iV + l)2(iV + 2)2 
64(57V^ + lliV^ + 50iV2 + 85iV + 20) 

27iV(iV + l)2(Ar + 2)2 
16(40iV^ + ISlA^s + 544A^2 ^ 779^ ^ 214) 

27A^(iV + l)2(iV + 2)2 
32(65A^6 + 429A^5 + 1155A^^ + 725A^3 + 370A^2 ^ 495^ ^ 543 



6S'_2,i — 5S-3 — 6S-2S1 
S2S1 

l_4QfiA^-KS4S'l 



81{N - l)m{N + l)2(iV + 2)2 
16(20A^^ + 107A^=^ + 344Ar2 ^ 439^ ^ -^^^^-^ g^^^^ 



S1N{N + iy{N + 2y ' S1{N -1)N^{N + lf{N + 2y 

32(47A^6 + 278A^^ + 1257A^^ + 2552N^ + UMN"^ + 2SAN + 448) 



-S. 



81N{N + 1)^{N + 2)3 
8(22Ar6 ^ 271A^^ + 2355A^^ + 6430iV3 + 6816iV2 + 3172Ar + 1256) 2 

81N{N + 1)3{N + 2f ' 
Qio(A^) ^ 448(iV2 + Ar + l)(iV2 + iV + 2) 



243(Ar- l)Ar2(Ar + l)4(iV + 2)4 9{N-l)m{N+iy{N + 2) 

16{5N* + 20iV3 + 59iV2 + 76N + 20) QuiN) 

-01Q2 



9Ar(iV + l)2(iV + 2)2 9(Ar- l)Ar3(iV+ l)3(iV + 2)= 

gi2(iv) 



5 



243(A^ - 1)N^{N + l)5(iV + 2) 

+n^r|C'^| ^^^^|j^'_^+^^+^^^^ [144^2,1,1 - 72^3,1 - 72^2,1^1 + 48^4 - 16^3^1 

-24.2^ - 12.2..^ - 2.^ - 9..^(2 + 42..C3] + 32 "^;;^7;^y^^y^ [3.2.1 + S!] 

128(iV2 - 3iV - 2) gi3(iV) 
-'-'2,1 — 



2 



3Ar2(7v+ l)(Ar + 2) ' 81(iV- l)iV3(iV+ l)3(iV + 2) 

gi4(iV) 32(10Ar4 + ISSAT^ + 789A^2 ^ 5217V + 141) 2 

.1 



27(iV- l)iV4(iV + l)4(iV + 2)3 81A^2(7v + l)2(Ar + 2) 

16(230Ar5 - 924Ar^ - bldSN^ - 7454A^2 _ io2i7Ar - 2670) 



243A^2(^ ^ l)3(7v + 2) 
16(5A^3^11Ar2 + 28Ar + l2) Q^^{N) 

"'^1^2 - 7^777 -.NArS/Ar , -,\S/Ar , ^\2''3 



9A^2(7v + l)(Ar + 2) 9{N -1)N3{N +lf{N + 2) 

C.^ .,, . , (176) 



9(A^ - l)m{N + 1)4(A^ + 2)3^ 243(A^ - 1)N^{N + lf{N + 2)5 
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Q9{N) 



with the polynomials 

Qi (N) = 15N^ + 90N^ + U6N^ - S2N^ - 501 A^^ - GlOA^^ - 2AAN^ - ASN^ 

+22AN + 96, (177) 
Q2{N) = lOA^^ + 105A^^ + 231Ar^ + 109A^^ + 143iV^ + 122A^ + 144 , (178) 
Qs{N) = llA^^ + 32A^^ + 363A^^ + 848A^^ + 186A^^ - 520A^ + 88 , (179) 
Qi{N) = 273A^^^ + 2565 A^^^ + 9838A^^° + 18902 A^'' + 14303A^^ - 8953A^^ 
-26402Ar^ - 19930A^^ - 1860A^^ + 4512A^^ - 2560A^^ - 2208A^ 
-576 , (180) 
Q^iN) = 3A^^ + 9A^^ - A^^ - 17A^=^ - 38A^^ - 28A^ - 24 , (181) 
Q6(A^) = 45Ari° + 405A^'' + 1606Ar^ + 3842A^^ + 6717Ar6 + 9325Ar5 + 10888Ar^ 

+9804A^^ + 6232A^2 + 3264A^ + 864 , (182) 
QjiN) = + mN^ + 150A^^ - 8A^^ - 357A^^ - 350A^^ - 300A^^ - 208A^ 

-144 , (183) 
Qs{N) = 873Ar^^ + 8298A^^^ + 29743Ar^2 ^ 38g927vii - 42545Ar^° - 211766A^'^ 
-227439A^^ + 95376A^^ + 573704Ar6 ^ 942576Ar5 + 995168A^^ 
+776576 Ar=^ + 527232 A^^ + 232704A^ + 48384 , (184) 
: 32A^'' - 936A^^ + 6448A^^ + 55208Ar^ + 126160Ar5 + 61760Ar4 

-53152Ar^ - 25024A^2 - 32256A^ - 13824 , (185) 
Q^q{N) = +7856A^^° + 84672Ar9 + 377648ArS + 985568Ar^ + 1395456A^*^ 
+470688Ar5 - 1183712Ar4 - 1180224Ar3 - 182528Ar2 - 42752Ar 
+13824 , (186) 
Qii (AT) = 60 AT*^ + SGOA^s + 584Ar^ - 12SN^ - 2mAN^ - 2440Ar^ - 976Ar3 

-192A^2 ^ ggg^ ^ ^ (^;L87) 

Qi2{N) = 28776iV^^ + 356112A^^^ + 1896088A^^''^ + 5538320A^^2 + 9112264Ari^ 
+6793968A^^° - 3019528A^^ - 11879520A^^ - 11673088A^^ 
-6450992A^6 - 3726976Ar5 - 2248128A^^ - 183296Ar3 + 268032A^2 
+ 147456Ar + 27648 , (188) 
gi3(Ar) = +464A^^ - 15616A^^ - 38112A^^ + 27776Ar5 + 146064iV^ + 119552Ar^ 

+ 109312A^2 ^ 86016A^ + 62208 , (189) 
Qia{N) = A56N^^ + 4376Ari° + 11328iV'^ - 3184ArS - 54552Ar^ - 111720Ar6 
-155376A^^ - 251072A^^ - 312192A^^ - 222mAN^ - 135936Ar 
-41472 , (190) 
Qi^{N) = 168iV^ + 672Ar^ + 784Ar6 - 3192Ar^ - 5600Ar3 - 7168A^^ - 4480A^ 

-2688 , (191) 
Qig(iV) = 90A^^^ + 630A^^° + 1592A^'^ + 1260A^^ - 1934A^^ - 8218A^6 

-15524A^^ - 23944A^^ - 26752A^^ - 18400A^^ - 11328A^ - 3456 , (192) 
Qn{N) = 15777A^^^ + 186525Ari^ + 879391A^^^ + 1874085A^^^ + 575913A^^^ 
-5568833A^^^ - 10465411A^^^ - 2970289A^^° + 11884298A^'' 
+12640320A^^ - 10343664A^^ - 40750480A^^ - 55711424A^^ 
-53947712A^^ - 42534912A^=^ - 23256576A^2 - 7865856A^ 
-1244160 . 5^ (193) 



We compared aJJ](iV), Eq. (176), to the fixed moments (H.28-H.32) of Ref. [95] and 
obtained agreement. The anomalous dimension appears in the 1/e term of (174). As all 
other contributions to this term are known, the anomalous dimension can be obtained by 
comparing with the 1/e term of the present computation. The following expression for 
m {^) is obtained: 



119 



{N + l){N + 2) 
Ca 



64 32 \ 128(57V^ + 8iV + 10) ^ 



S2 



64(5iV^ + 26iV3 + A7N^ + ASN + 20) 

9N{N + 1){N + 2) 
64(5Ar4 + 2QN^ + AIN^ + 497V + 20) ^ 



9N{N +1){N + 2) 

6AQ^s{N) 



160i9(iV) 



27iV(iV+ l)2(iV + 2)' 



27(iV - l)A^'i(iV + l)3(iV + 2)3 



QisiN) 

Ql9{N) 



Q20{N) 



32{N^ + AT + 2) 3 320(A^2 + tv + 2) 



9A^ 

32(A^2 + iv + 2) 
3iV 



9N 



32{5N^ + 3N + 2) 



32(10A^3 + i^]sj2 ^ 29iV + 6) 
97V2 



32(47A^^ + 145A^3 + 426A^2 ^ 4^i2N + 120) 
27A^(7VTl) 
4g2o(A^) 



Si 



27(N - 1)N^(N + l)4(Ar + 2)3 

19iV6 



165iV^2 + 1485A^^^ 



124iV^ + 492A^^ + 1153iV3 + 1362iV2 ^ ^UN + 152 , 
5194iV^° + 8534A^^ + 3557iV^ - 88997V^ 
-10364iV^ + 6800 A^^ + 25896A^^ + 30864Ar3 ^ i9904Ar2 
+7296Ar+ 1152 , 

99Ar^^ + 990A^^3 + 4925A^^2 + 17916A^^^ + 46649A^^° + 72446A^'^ 



+32283A^^ - 95592A^^ - 267524iV^ - 479472A^^ 
-455168A^3 - 269760A^2 - 122112A^ - 27648 . 



586928Ar^ 



(194) 
(195) 

(196) 

(197) 



It agrees with the moments (H.1-H.5), [89,90,92,94,95]. Due to the algebraic compacti- 
fication we obtain a lower number of harmonic sums Ss{N) if compared to Ref. [94], and 
agree with [201]. 
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lPS,(3) 



8.2 The pure-singlet contribution 

The general structure of the pure-singlet contribution A^^^^^ is, [95], 



^(3),PS 



■,(o).,(o) 



6£3 



7, 



1/ 47^''"= 



39 



/So + /3o, 



Q 



(0) (1) (0) 

7gg 7qg _|_ 7gg 



2^(1) _ ^(1) 



. /^(2),PS 



;;(2),PS 
Iqq 



-(0) (0)> 

, ^(0)^(2) ^(0)^(2) ./^ , o ^„(2),PS 799 799 (,2 
+^99 S9,Q ~ Tgq^Qg ~ ^^P^ + M^Qq " 



16 



7S-7i? + 6/3o 



(198) 



In this process the photon couples to the heavy quark line. Summing the individual 
results of the different diagrams (D.201-D.206) weighted by their respective multiplicities 
yields 



.PS,(3) 



1 256 



A^2(_i + iV)(2 + iV)(l + Ar)2 £3 9 



1 

+- 

e 

64 
+7^ 



128 ,,,2 , , o^2o 128 



9 



{N' + N + 2ySi + — 



{N'^ + N + 2f 

Q2l{N) 

27 N{l + N){2 + N) 



g2i(iv) 



Q22{N) 



+a 



The constant part is given by 



27 A^(l + A^)(2 + A^) ' %l N'^{l + Nf{2 + Nf 

PS,(3) 
Qq 



(199) 



PS,(3) 
''Qq 



T^UfCp 



N^{1 + Ny{2 + N){N -1) 



{N^ + N + 2f 



1760 



■ ^ ^2 5i - ^ - y ^1 C2 + ^ C3 



208 



Q2i(iV) 



^2 + 



16 



Q2i(iV) 



27 iV(l + iV)(2 + AT) ^ 27 iV(l + iV)(2 + iV) 



+■ 



16 



Q2l{N) 



:C2- 



32 



g22(iv) 



9 A^(l + A^)(2 + A^)^ 81 N'^{l + N)^{2 + N)^ 
32 Q23(iV) 



^1 



243 + Ny{2 + N)^ I ' 



(200) 



with 
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Q2i{N) = 8N'^ + 37 + 68N^ -IIN^ -86N^ -56N^ -lOAN -48 , (201) 
Q22{N) = +25 + 176 + A17 + 30 N'^ - 20 + 18A8 

+2244 + 1648 A^^ + 3040N^ + 2112iV + 576 , (202) 
Q23(iV) = 158 iV^3 + 1663 A^i2^ 7714 AT" + 23003 iVi° + 56 186 AT^ 
+89880 + 59452 A^^ - 8896 A^^ - 12856 - 24944 A^^ 
-84608 A^^ - 77952 A^^ _ 35712 N - 6912 , (203) 

agreeing with the moments (H.33-H.38) from [95]. 



8.3 The pure-singlet contribution g 

The analytic structure for this OME, [95], is given by 



X(3),PS 
A 

^qq,Q 



/ ^2 3^/2 



3£3 



^(0)^(1) 

ha 'm 



+ - 



i{2),PS 
Iqq 



-(0) (0)0 A \ ^V"^'' 

_L n,(0).,(2) _ Iqg Igq Po,q(,2 \ "gg,Q 
7oo "00,0 ^ 



(3),PS 



3 ■ 'qg "'gq,Q 



(204) 



For this OME we computed the complete result, since no other color factors than TpUfCp 
contribute in this case. The matrix element corresponds to the pure-singlet term in which 
the photon couples to a massless fermion line. One obtains 



A 



PS, (3) 
qq,Q 



T| UfCp J 1 256 

Ar2(l + Ar)2(2 + Ar)(iV- 1) I ~£3~9~ 



+- 



+- 



e' 

1 

e 

12 



{N^ + N + 2f 



(N^ + N + 2fSi , 

9 ^ ^ ^ 1 27 N{2 + N){l + N) 



{N^ + N + 2f^ 

B Q24{N) 



128 128 2 32 



-Si- 



64 



Q25{N) 



27 N{2 + N){l + N) 81 N'^ {2 + N)^{1 + Ny 



+a 



PS, (3) 
qq,Q ' 



(205) 
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with the constant part 



- 7V2(-l + 7V)(2 + 7V)(l + iV)2|^^ +^ + 2^ [-W^'^^^'^' 

128 224 \ 64 g24(iV) 

27 iV(2 + iV)(l + iV) ^ 9 iV(2 + iV)(l + Ar)^^ 

,64 g25(iV) 32 Q,,iN) 

81 Ar2(2 + Ar)2(i + Ar)2 1 243 Ar3(2 + Ar)3(i + Ar)3 r ' ^ 



and 



Q24{N) = 16 iV^ + 74 iV^ + 181 A^s + 266 AT^ + 269 A^^ + 230 TV^ 

+44 - 24 , (207) 

Q25{N) = 181A^^° + 1352A^9 + 4737Ar^ + lOlOlAT^ 

+14923 + 17085 A^^ + 14133 A^^ + 5944 + 568 AT^ - 48 A^ 
+144 , (208) 

Q26{N) = 2074 A^^3 + 21728 A^^2 + 105173 Ar^^ + 311482 A^^° + 636490 AT^ 
+966828 A^^ + 1126568 A^^ + 968818 A^^ + 550813 A^^ 
+169250 A^^ + 12104 A^^ - 3408 A^^ - 1008 A^ - 864 , (209) 

agreeing with the fixed moments (H.39-H.45) from [95]. From (198, 199, 204, 205) one 
may extract the anomalous dimension 



(N -1)N^{N +1)^{N + 2) 



3 ' ' J \ L ' ' g N{l + N)(2 + N) 

64 Q28{N) 
27 N'^{l + Ny{2 + Ny 



(210) 



Q27{N) = 68N^ + 37N^ + 8N^ -IIN^ -86N^ -5QN^ -lOAN -A8 , (211) 
Q2s{N) = +52 A^^° + 392 A^9 + 1200 A^^ + 1353 A^^- 317 AT^- 1689 AT^ 

-2103 A^^ - 2672 A^^ - 1496 -48N+ 144 . (212) 

Again we obtain agreement with the fixed moments from [95], Eqs. (H.6-H.12). 
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8.4 The non-singlet contribution A^^ q 

The following structure is known for the NS contribution: 

2 f /I (O)c. 1 /o (l).NSc. .-(1),NS 



X(3),NS 



/X^ / I 3£^ 



(/3o + 2/3o,q) + ^| 



+ 



1 / '^(2),NS 



, Im PoPo,qC,2 ^r. (0) r (-l)^(l),Ns\ , {3),NS 



3 h + 

(213) 



cf. [95]. We obtain the following result: 



.NS,(3) 

^qq,Q 



512 128 3iV + 3iV2 + 2 
(1 + Ar)Ar 



256 1280 32 20Ar + 47iV2 + 6A^3 + 3Ar4-12 
02 — — + — 



27 
128 



81 



81 



64 640 



(1 + Ar)2Ar2 
16 3 iV + 3 + 2 

¥ (l + Ar)Ar 



C2 



1280 



Si 



^_8_ g29(A^) 



81 (1 + NfN^ 



+ a 



NS,(3) 
99,Q ' 



(214) 



with 



NS,(3) 



_160 112 3A^ + 3iV2 + 2 640 

C2 -bi - — (1 + jY)jY C3 + ^ -^2 

4 20A^ + 47iV2 + 6A^3^3iV'^- 12 55552 
^27 (1 + Nfm ~ 729 



5i 



and 



QsoiN) 



2 Q3o(iV) 
729 (l + iV)4iV4 



321 A^^ + 963 A^^ + 1307 N"^ + 833 A^^ + 152 A^^ - 16 + 24 , 
+11751 A^^ + 47004 A^^ + 93754 A^^ + 104364 A^^ + 55287 A^^ 
+6256 - 2448 A^^ - 144 AT - 432 . 



(215) 



(216) 



(217) 



Again the anomalous dimension for the non-singlet case is determined by comparing the 
1/e term of our result to the corresponding term of the general structure. We obtain 



MS _ r T'n 1^2^ S 9 + ^ ^31 (iV) 



} , (218) 
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with 



Q3i{N) = 51N^ + 153N^ + 57N'^ + 35N^ + 96N^ + 16N-24: . (219) 

For fixed moments afgQ^ agrees with the values in [95] , cf. (H.46-H.60) and the anomalous 
dimension 7^^''^^ with (H.13-H.27), [90,93,95]. 



-NS,TR 



8.5 The non-singlet transversity contribution -^qq^Q 

The non-singlet transversity OME A^^q^, cf. Ref . [179] , obeys the same general structure 
as the non-singlet OME A^^q, (213). One obtains 



4(3),TR 



128 512 



9 



27 



Si 



1 



256 1280 32 
- -g^^i + 27 



1 

+ - 

e 



o 64 640 1280 16 



8(107Ar2 + 107A^ + 8) 
27iV(Ar+ 1) 



a 



{3),TR 
qq,Q ' 



(220) 



with 



,(3),TR 



qq,' 



TpUfCi.- 
112 



320 160 



640 4 55552 
9 + ^"^^ + 9^^ " ^29""^^ 
2(3917A^^ + 78347V3 + 4157A^2 _ 43^ _ ^44) 
243A^2(i + TV) 2 



(221) 



agreeing with the moments (H.74-H.86) for fixed values of N. Extracting the anomalous 
dimension from the 1/e term of the OME yields 



^i2),TR r T^r, J c 640 128 8 (17iV^ + 17iV - 8) I 

The results for the anomalous dimensions constitute a first independent check of the result 
obtained in [179,203]. It is interesting to note that for this color factor the vector- and 
tensor operators lead to the same structures in the harmonic sums for a^g^Q and 7^5^. 



8.6 Harmonic Sums 

The T|.n/CF^^-contributions at O(a^) to the massive operator matrix elements contain 
nested harmonic sums up to weight w = 4. This also applies to all individual Feynman 
diagrams, which wc calculated in the Feynman-gauge, cf. Appendix D. In intermediary 
results generalizations of harmonic sums occur, see Section. 8.7. As has been observed in 
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the computation of various other physical quantities before, such as anomalous dimensions 
and massless Wilson coefficients to 3-loop order [93,94,109,201], unpolarized and polarized 
massive OMEs to 2-loop order [120,121,123,204], the polarized and unpolarized Drell- 
Yan and Higgs-boson production cross section [128], time-like Wilson coefficients [129], 
and virtual- and soft corrections to Bhabha-scattering [130], the classes of contributing 
harmonic sums are always the same. For main properties of the nested harmonic sums see 
Appendix F.3. They depend on the loop-order and the topologies of Feynman diagrams 
involved. 

In the present case the following harmonic sums emerge : 

Si 

52, S-2 

53, S^3, 5'2,l, <S'_2,1 

Si, S-4, Ss^i, S-3^i, S-2,2, 'S'2,1,1, S-2,1,1 ■ (223) 

Note that this class, as for the other processes mentioned above, does not contain the index 
{— 1}. Moreover, we used the algebraic relations between the harmonic sums, cf. [131]. 
Furthermore, structural relations exist between harmonic sums, cf. [132,133,202], which 
reduce the set (223) further. Here the sums 

S-2,2, .53,1 (224) 

are connected by differential relations w.r.t. their argument N to other sums of (223). 
This is also the case for all single harmonic sums S±n, n E N, n > 1, using both the 
differentiation and argument-duplication relation, cf . [124] . Due to this represents the 
class of all single harmonic sums. I.e. only the six basic harmonic sums 

Si 

'S'2,1, 'S'-2,l 

S-3,1, 'S'2,1,1; 'S'-2,l,l (225) 

are needed to represent the 3-loop results for the Tj.n/ CF,A^contributions to the OMEs 
calculated in the present paper. In the final representation we refer to the algebraic basis 
(223) and consider the basis (225) for a later numerical implementation. We sorted the 
respective expressions keeping a rational function in N in front of the harmonic sums 
(223) and ^-values, like (2 and (3. 

The harmonic sums emerge from the series-expansion of hypergeometric structures 
like the Euler B and F-functions and the Pochhammer-symbols in the (generalized) 
hypergeometric functions pFQ{ai,bi;l) in the dimensional parameter e. This leads to 
single harmonic sums first, which, through summation, turn into (multiple) zeta values 
[135] and nested harmonic sums [124,125]. The principle steps on the way from single- 
scale Feynman diagrams to these structures have been described in Ref. [132]. 

8.7 Generalized Harmonic Sums 

At 3-loop order the expressions for individual Feynman diagrams are rather large. In 
Section 7.2 we described how the Feynman parameter integrals which emerge in the 
present calculation are transformed into nested infinite and finite sums. If these sums 
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could be computed analytically as a whole only nested harmonic sums would occur in 
the calculation. However, this is not always possible in practice. Usually the expressions 
obtained are split into different parts and the sums are then computed. In intermediary 
steps, depending on the summation methods used, more complicated sum-structures 
emerge. We will not deal with this aspect here, but only refer to typical structures, which 
occur in the final result separating the expressions to be summed over into tractable terms 
of an intermediate size. 

Here, so-called generalized harmonic sums occur [205,206]. They obey the following 
recursive definition [205] : 




+Sm^+rn2,m3,...{^l ' 3^2, X3, N) . (226) 

The sums S may be reduced to nested harmonic sums for e {— 1, 1}. In the present 
calculation the values of ,Xj extend to {—1/2, 1/2, —2, 2}. These sums occur in ladder like 
structures, cf. [109,180,195], but may also emerge, if contributions to 3-loop Feynman 
diagrams containing a 2-point insertion, are separated into various terms. They were 
even observed in case of massive 2-loop graphs of the type [121] if large expressions are 
arbitrarily separated, cf . [195] . The weight of these sums can reach w = 5 intermediary, 
depending on the e-structurc of the contribution, although only w = 4 sums will emerge 
in the final results. Examples for these sums are : 

Sj{l/2,N), S2(-2;N), ^2^(-l, 2; TV), Ss,,{-2, -1/2; N), 
5i,i,i,2(-l, 1/2, 2, -1; AT), 52,3(-2, -1/2; N), 

^2,2,i(-l,-l/2,2;iV), etc. (227) 

The algebraic and structural relations for these sums arc worked out in Ref. [206]. 
Similar to the case of harmonic sums, corresponding basis representations are obtained. 
These relations were provided in a code used in the present calculation. They finally 
lead to the reduction of the results for the individual diagrams to a representation just in 
terms of nested harmonic sums. In Appendix G one of the typical lengthy nested sums, 
which also contains the generalized harmonic sums, is shown. 
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9 Summary 



In the present work a first contribution to the computation of the massive Wilson coeffi- 
cients of the deep-inelastic structure function F2{x,Q'^) for unpolarized charged lepton- 
nucleon scattering at 3-loop order for general values of the Mellin variable N in the 
region Q'^>10 ■ has been made. Since the corresponding massless Wilson coefficients 
are known, cf . [109] , only the massive operator matrix elements remain to be calculated. 
A series of fixed Mellin moments for all contributing OMEs has been computed before, 
Refs. [95, 179]. Here, we extend this work to general values of N, calculating a very large 
class of terms for A^qI{N) and aU contributions to ^^/^^(Ar) A^^^;^\n) A^g '^^\N), and 
^^^qQ^'^^^i-^) for the color factors TpUfCp and TpUfCA- The complete results on Ag^^lN) 
were obtained in a larger team and will be given, including A^^^q{N), in [180] in detail. 

The present calculation is a first step within a larger programme to compute all topolo- 
gies contributing to the massive OMEs, including those which define heavy-hght quark 
transitions in the variable flavor number scheme. Due to the large heavy flavor contribu- 
tions of 0(25 — 30%) to the nucleon structure functions the experimental precision of the 
DIS World Data requires these calculations, both to extract the different parton distri- 
bution functions and to measure the strong coupling constant q;s(M|) at highest possible 
precision. The knowledge of these quantities is instrumental for precise measurements at 
the LHC. q;s(M|) is furthermore one of the central parameters in physics. 

The present calculation built technically on a method to compute the Feynman- 
integrals, containing one massive line, directly, i.e. avoiding traditional methods, as 
integration- by-parts and related techniques. Those usually lead to a large proliferation of 
terms compared to the compact results being finally obtained. Instead we referred to rep- 
resentations in terms of generalized hypergeometric functions. They occur as the analytic 
result of the Feynman parameter integrals. Due to this the expansion in the dimensional 
parameter e can be uniquely performed in an elegant way. At the same time, various 
significant simplifications of intermediary results are possible choosing particularly this 
representation. In the end, a small number of nested finite and infinite sums over hyper- 
geometric terms equipped with Beta-functions and harmonic sums has to be performed. 
In the present calculation up to triple sums occurred. The main technical work to be 
performed consisted in finding these representation by means of computer algebra for the 
contributing graphs in a partly automated way. Intermediary very large expressions had 
to be handled and corresponding codes based on the systems FORM and MAPLE had to 
be designed. 

In the present case the final sums, compared to foregoing massless computations of 
other groups at 3-loops, cf. e.g. [93,94], turn out to be much more involved due to 
the the increased nesting in the massive case. Moreover, it has been unclear whether 
the usual nested harmonic sums form a frame to express all the intermediary results. 
This had to be found out using strict mathematical methods based on the construction 
of sum- and product fields. The corresponding software Sigma [142, 143] could be used 
for this purpose. Indeed, it turned out that in the intermediary results the algebra of 
harmonic sums is to small to describe the corresponding terms and generalized harmonic 
sums had to be invoked. The methods used would have been pointing to other more 
general structures also, if contributing. Thus this method is indispensable in exploring new 
territories in higher order quantum-field theoretic calculations. The present work triggered 
mathematical research to find all relations for this new class of functions, cf. [206], which 
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will help to make the codes being used at present even more efficient. 

The major new results of the present work are the constant parts of the unrenormal- 
ized massive OMEs at 3-loop order for general values of N, a^^\N) in 0{TpnfCF,A). 
Furthermore, the corresponding contributions to the 3-loop anomalous dimensions in the 
vector- and transversity case are obtained. In this way, a first independent recalculation 
of these important quantities given in [93, 94, 203] before, was performed using very dif- 
ferent methods. Here the anomalous dimensions are computed in the massive case. We 
compared the present results mutually with all the available fixed Mellin moments in the 
literature and find agreement. In case of the anomalous dimensions the comparison could 
be performed for the general expressions. The representation we gave, however, is more 
compact, since we applied the algebraic relations between harmonic sums. 

The final results for the matrix elements a\f (N) and the 3-loop anomalous dimensions 

^ij^ (N) are given in harmonic sums only. All generalized harmonic sums cancel already at 
the level of the individual diagrams. The complexity of harmonic sums determining the 
present results is maximally six for the ajf (iV) and two for 7jf (iV) at 0{T^nfCF,A), 

usmg 

also the structural relations, Ref. [132]. This situation is comparable to the case of all 
massive OMEs at 2-loop order expanded in the dimensional parameter to 0{6), cf. [121]. 
However, the present calculation, at 3-loop order, has been by far more complex and had 
to pass much more sophisticated structures intermediary. It is interesting to note that 
the structure w.r.t. harmonic sums for the vector- and tensor fiavor non-singlet operators 
concerning both the matrix elements a'l^\N) and the contributions to the anomalous 

dimensions 7^^^ (N) are the same for the terms given by the harmonic sums in case at the 
present color factors. 

The present calculations unraveled the relevance of hypergeometric structures in com- 
puting Feynman integrals in general. They are given in the present by the generalized 
hypergeometric functions emerging. Likewise, the finite and infinite sums had to be 
performed over hypergeometric summands, equipped with products out of (generalized) 
harmonic sums. 

The present work constitutes a first step to explore the general structure of single 
scale massive observables in QCD at the 3-loop level, which probably bears a rich host 
of yet unexplored structures. They form an interesting topic to be studied in the future 
to understand the final simplicity of seemingly complex problems at the higher loop level 
in Quantum Field Theories such as Quantum Chromodynamics, and hopefully finally the 
yet unknown reason for that. 
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10 Appendix 



A Conventions 

We use natural units 

h=l , c=l, £0 = 1, (A.l) 

where h denotes Planck's constant, c the vacuum speed of light and Eq the permittivity 
of vacuum. The electromagnetic fine-structure constant a is given by 

a^aif, = 0) = = — ^ 137.03599911(46) ' ^ 

In this convention, energies and momenta are given in the same units, electron volt (eV). 

The space-time dimension is taken to be D = 4 + £ and the metric tensor g^j^^ in 
Minkowski-space is defined as 

S'oo = 1 , S'n = -1 ,i = 1 . . .-D - 1 , Qij^O . (A.3) 

Einstein's summation convention is used, i.e. 

^l=o 

Bold-faced symbols represent {D — 1) -dimensional spatial vectors: 

X — {xq, x) . (A. 5) 

If not stated otherwise, Greek indices refer to the D-component space-time vector and 
Latin ones to the D — 1 spatial components only. The dot product of two vectors is defined 

by 

D-l 

p.q = poQo - ^ PiQi ■ (A.6) 

i=l 

The 7-matrices 7^ are taken to be of dimension D and fulfill the anti-commutation 
relation 

{7^,7,} = 2y^, . (A.7) 

It follows that 

l^Y = D (A.8) 
Tr(7^7,) = 4^^, (A.9) 

Tr {-f^-f^-fa^/s) = ^g^ci.gaP + Qnpgi^a- gfiaQup] ■ (A. 10) 

The slash-symbol for a D-momentum p is defined by 

iZ-^l^p"- (A.ll) 
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The conjugate of a bi-spinor u oi a, particle is given by 

u = ii^7o , (A.12) 

where f denotes Hermitian and * complex conjugation, respectively. The bi-spinors u 
and V fulfill the free Dirac-equation, 

(j/ — m)u{p) = 0, u{p)(j/ — m) = (A. 13) 

^+m)v{p) = 0, w (p) + m) = . (A.14) 

Bi-spinors and polarization vectors are normalized to 

'^^u{p, a)u{p, a) = jZ+rn (A. 15) 

^v{p,a)v{p,a) = pj-m (A. 16) 

Y,e\K\)e\K\) = (A.17) 
\ 

where A and a represent the spin. 

The commonly used caret " ' " to signify an operator, e.g. O, is omitted if confusion 
is not to be expected. 

The gauge symmetry group of QCD is the Lie-Group SU (3)c. We consider the general 
case of SU{Nc). The non-commutative generators are denoted by i", where a runs from 
1 to — 1. The generators can be represented by Hermitian, traceless matrices, [158]. 
The structure constants f"''"^ and d"''"^ of SU{Nc) are defined via the commutation and 
anti-commutation relations of its generators, [160] , 

[t^^t^] = if'^bc^'' (A.18) 

{t'^A = rt'^ + ^^^fe . (A.19) 

The indices of the color matrices, in a certain representation, are denoted by i,j, k, I, ... 
The color invariants most commonly encountered are 

SabCA = r'^f'"' (A.20) 

SijCp = ttitf^ (A.21) 

6a,TF = tt,t\i- (A.22) 

These constants evaluate to 

Ca = N,.Cp = ^Tf = \- (A-23) 
At higher loops, more color-invariants emerge. At 3-loop order, one additionally obtains 

ci"'"46c = {Nl - l){Nl - A)jN, . (A.24) 
In case of -SC/(3)c, Ca = 3 , Cf = 4/3 , d^-^^dabc = 40/3 holds. 



65 



B Feynman Rules 



The Feynman rules can be derived, as in any renormalizable Quantum Field theory, from 
the path integral representation [207-210]. For the QCD Feynman rules. Figure 14, we 
follow Ref. [160], cf. also Refs. [211,212]. D-dimensional momenta are denoted by Pi 
and Lorentz-indices by Greek letters. Color indices are a,b,... and i,j are indices of the 
color matrices. Solid directed lines represent fermions, wavy lines gluons and dashed lines 
ghosts. Arrows denote the direction of the momenta. A factor (—1) has to be included 
for each closed fermion- or ghost loop. 




1^, b, P2 




P, C,P3 



11, a 



h/ "^-c^P 



-gsr^P, 



1/, b 




(T, d 



p,c 



+r'n,.gp. - gp.g.p\ 
+r'''n^pg.a-gp.gpa] 



i P j 



m+iO 'J 



a P h 



Figure 14: Feynman rules of QCD. 
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The Feynman rules for the quarkonic composite operators are given in Figure 15. Up to 
0{g^) they can be found in Ref. [81] and also in [213]. Note that the 0{g) term in the 
former reference contains a typographical error. In Ref. [95] these rules were checked and 
agree up to normalization factors, which may be due to different conventions. There also 
the new rule with three external gluons was given. The terms 7-1- refer to the unpolarized 
(+) and polarized (— ) case, respectively. Gluon momenta are taken to be incoming. 



P,J 



N-l 



N > 1 



P2,3 



gt]A'4^1± Ef=o'(A • pi)^(A • P2)^-^-' , N>2 



jj,, a 




N-l-2 



{tH'),,{Ap, 



iV > 3 



(tV),,(Api 



Ap, 




P3,y",a pi.-u.b P5,P,c 



3A^A'^A^47± Ef= 


-4 v^TV- 

=0 2^i=j 


3 ^N- 
+ 1 Z^m= 


:?+i(A.P2)^- 


(A.pi)^-™- 


-2 


{tHH^),,{A.p^ + 


A.P5 + 


A.piY 


-^■-i(A.p5- 


f A.pi)"^-'- 


-1 


+{m')ji(A.p4 + 


A.P5^ 


-A.pi) 


'-^■-HA.p4 


+ A.pi)"-' 


-1 


+{tHH'^)jiiA.ps + 


A.P5^ 


-A.pi) 


'-'-\A.p, 


+ A.pi)"-' 


-1 


+{tHH-),,{A.ps + 


A.p5i 


- A.pi) 


'-'-\A.p^ 


+ A.pi)"-' 


-1 


+(m^),i(A.p3 + 


A.p4^ 


-A.pi) 


'-^-\A.p, 


+ A.pi)'"-' 


-1 


+{t'^tH'')jiiA.ps + 


A.P4 4 


-A.pi) 


'-'-\A.p, 


+ A.pi)"^-' 


-1 



7+ 



7- = 75 • 



iV > 4 

For transversity, one has to replace: 



Figure 15: Feynman rules for quarkonic composite operators. A denotes a light-like 4-vector, 
A^ = 0; is a suitably large positive integer, Ref. [95] 
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C D-dimensional integrals 

In the calculation of the D-dimensional loop integrals [214-216], D — e, we perform 
first a Wick-rotation to Euclidean momenta 



M 

X 



1 

n r_t2 _ 



A(-A;f-mf)«»' 

with Vcj e N. One obtains the following Euchdean integrals, where, cf . [160] , 



kj^ — -\- k , 



{27r)D {k^ - ip)m H ) J (27r)^(A;| + v5)' 



^ r(r + D/2) r(m - r - D/2) 

(47r)^/2 r(D/2) r(m) ^ 



/ 



m, r e N 

gCtia2 ga3a4 _j_ gCtias ga2a4 _j_ ^aia4ya2a3 



(C.l) 

d^'k k'^k'^ fik") = ^-jy- / d^k k^ f{e) (C.2) 



d Jh k-j^ ^3 ^4 ) 



/ 

y d^A; (C.3) 

y d^fc k'^k^k^k2k^k^ f{k') = j^s+iZ^+sD J ^""^ ^^-^^ 



+ 

/ , n , , /I /. 

X 



with 



Pg^^Qjaj j _ gaia2 ^ga3a4 gO^ae _^ gasas ga4ae, ^ gOsaa ^0405-^ 
_j_^ai03 ^ga2a4 gOisae _|_ ^a2a5^a4Q6 _|_ ^«2a6^Q!4a5j 
_j_^ai«4 j^^a2a3^ct5a6 _j_ gOiar, ga2ae _|_ ^a3Q6^a2a5j 
_j_^Q:io:5 |'^03a4^02a6 _|_ ^«2a3^a4a6 ^ ^tt3a6^Q2a4j 
_|_^aia6 ^gCi2.oi4 gCi2ai _j_ gOtzct^ gOi20i4 _|_ ^a2a3^a4a5j 



and 



2M+1 



d^A; A;f /(A;2) = . (C.6) 
For each loop integral a universal factor 



e 



5, = exp (7B-ln(47r)) -J (C.7) 
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emerges, where 7e denotes the Euler-Mascheroni constant 



= hm 



* 1 

1=1 



(C.8) 



The factors arc kept separately and are not expanded in e. In the MS-scheme [97] 
they are set to = 1 at the end of the calculation. 
The F-function obeys the relation 



r(l + £) = exp 



n=2 



(C.9) 



with 



Cn = E^' ^e^' ^^2, 



(C.IO) 



fe=i 



the Riemann ^-function at integer arguments. 

We apply the following Feynman parametrization to combine denominators 



„1 „i / n N 

r(n) / dxi . . . dXjiS I 2_] ~ ^ 

•^0 -^0 \fe=l y 



(xiAi + ... + XnA^y 



(C.ll) 



resp. 



AT---A-r. 



. . . 



/ dXnS Xfe - 1 

•^0 Vfc=i 



11^=1 



;a;iAi + ... + x„A„)(^^=i"'=) ' 



(C.12) 



with Vaj e N. 

The integral over the 5-distribution yields 



„1 / n \ „+oo / n \ / n \ n 

/ dX/ 5 Xfc - 1 = / (iX; 5 Xfc - 1 U' 1 - ^ Xfe J]^ e{Xm) 

•^0 \fc=l / \k=l J \ k=l J m=l 



q 1 - E n ^(^-) ' 



(C.13) 



k=l,k^l / m=l,m^l 

where 9{z) denotes the Heaviside function 



e{z) = 



1, ^ > 
0, z<0 



(C.14) 
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D Results for the Individual Diagrams 



In this appendix we hst the rcsuhs for the individual diagrams contributing to the 



O {TpNfCp^A) terms in the massive operator matrix elements ^g^, 



(3),PS .(3),PS 4(3),NS 



f3l XR, 

and Aqq'q ■ They are all represented in terms of harmonic sums, despite the fact that 
in intermediate results also generalized harmonic sums occur. No sums beyond w — A 
contribute. Furthermore, the individual diagrams do not contain sums, which fully cancel 
in the final result, similar to [120, 121]. The individual contributions were compared to the 
corresponding moments in the calculation Ref. [95,179]. The common pref actor 



has been taken out. 



(D.l) 



D.l A 



Qg 



;(3) 



The results for the individual diagrams contributing toA^^ are: 



64(7V- l)(7V + 2) 
9 



64(2iV4 - 5iV3 + 7N^ + 2N + 3) 
27A^(7V+ 1) 



1 

+- 

s 



on Q 

- — (TV - 1)(N + 2)S2 - -(N - 1)(N + 2)C2 



56, 



32Pi(iV) 



81Ar2(iV + l)2(iV + 2)2 



-16(7V - 1){N + 2)Ss + —{N -1){N + 2)C3 

9 

32(2Ar4 _ 5^3 ^ ^2N + 3) 8(2N'^ - 5N^ + 7N^ + 2N + 3) 



9A^(A^ + 1) 
16P2(A^) 



9A^(A^ + 1) 



243A^3(7v + l)3(Ar + 2)3 



(D.2) 



with 



Pi{N) = 5N^ + 26N^ - 203N^ - 721N^ - 1079N'^ - mN^ + 269N^ 

+ 180 , (D.3) 

P^(N) = 17iV^^ + 107iV^° + 824A^^ - 927iV^ - 8904iV^ - 17883iV^ 
-19301iV^ - 101937V^ + 39737V^ + 15342Ar2 + 110527V 

+2808 . (D.4) 



„ 11 



128 ^ 128 



1 



64 „2 ^ 128 „ 



128(4iV2 + 21Ar + 8) 64(lliV2 + 51iV + 22) 
+ . .w,. . 1- 27{N+l){N + 2) 



27(7V + l)(7V + 2) 



70 



1 

+- 

e 



64 64 736 64 16 



64(4iV^ + 21N + 8) 2 64(A^^ - 15A^ + 2) 16 
^ 27(A^ + l)(A^ + 2) ^ ^ 27(iV+ l)(iV + 2) ' ^ Y^' 

32P3(A^) ^ , 32P4(Ar) 



16o4,32 352 128 112 64 32 



1 -^)-^ 



32.. 2 8. o2^16>. 64 ^, , 64(4iV^ + 21iV + 8) ^ 
" Y^^^^ " 3^'^i + y^^^' " y^^'^*^^ + 9(iV + l)(7V + 2) 

64(4Ar2 ^ 21iV + 8) 16(4iV^ + 21Ar + 8) 

^ 27(iV + l)(iV + 2) ' ' ^ ~9(iVTl)(iVT2)~ 
^ 64(4iV2 ^ 21iV + 8) 3 32(11A^2 ^ 240iV + 22) 

"'-'l oi / Ar , 1 \/ 7\r , o\ '-'3 



81(iV + l)(iV + 2) ' 81(iV+ l)(iV + 2) 
112 8(llAr2 + 5iiV + 22) 8P5(A^) .2 



9 ^ 9(A^ + l)(A^ + 2) ^ 81A^(A^+ l)2(iV + 2)2 



81A^(A^ + l)2(Ar + 2)2 ^ 243A^(A^ + l)3(Ar + 2)3 
16P8(iV) 



243(A^ + l)4(Ar + 2)3 



(D.5) 



P^{N) = 88A^^ + 780A^^ + 1873A^^ + 1938A^^ + 892A^ + 216 , (D.6) 

P4(A^) = 121 A^5 + 1126A^^ + 3280Ar=^ + 3907 A^2 + 1792Ar + 484 , (D.7) 

P^{N) = 176A^^ + 1614A^^ + 4097A^^ + 4686A^^ + 2324A^ + 648 , (D.8) 

Pq{N) = 44A^^ + 282 A^^ + 2057A^^ + 3966 A^^ + 1796 + 648 , (D.9) 
P7(A^) = 824A^^ + 10350A^^ + 46281A^^ + 104256A^^ + 124746A^^ 

+77364Ar2 ^ 24952A^ + 4752 , (D.IO) 
P^{N) = 1187A^^ + 15182Ar6 + 70743A^^ + 161517A^^ + 197367A^^ 

+137022Ar2 + 56836A^ + 9496 . (D.ll) 



(3),Ci 

Qa 



16(13A^2 + 65iV + 6) 



1 

+- 

e 



3(iV + 2)(iV + 3) 

64C 4A(Ar) 

' 3Ar(Ar + l)2(Ar + 2)2(Ar + 3)j 

2(13iV2 + 65iV + 6) 8(13iV3 + 86iV2 + 191iV + 102) 
' ^ (Ar + 2)(iV + 3) (TV + 1) (AT + 2) (AT + 3) ' 

+ Mm \ (D12) 

3A^2(^ + l)3(^ + 2)3(A^ + 3) J ' ^ ■ ' 



71 



with 



Pg{N) = 91 A^^ + 855 A^^ + 2481 A^^ + 3037A^^ + 1436A^^ - 28A^ - 48 , (D.13) 
Pio{N) = 569A^^ + 7152A^^ + 33604A^^ + 79234A^6 + 99271A^^ + 61206A^^ 

+13356A^^ + 584Ar2 ^ ggg^ ^ gg^ p -^4^ 



T^UfCf 



N{N + 2){N + 3) £2 



16(iV2 + 5A^ + 42) 



1 

+- 

e 



16{N + 3){N^ + 4Ar3 + 41iV2 + 32iV + 12) 



3iV(iV + 1)2 



^1- 



3{N +iy{N + 2) 



+2{N^ + 5N + 42)C2 + 



.q2 



8Pi2(iV) 



4Pi3(A^) 



3A^(iV + l)2 ' 3N{N+iy 
Pu{N) 



3iV(iV+ l)3(A^ + 2) 3(iV + l)4(iV + 2) 



(D.15) 



with 



Pii{N) = 7N^ + 70N^ + 748N^ + 2770N^ + 3961N^ + 2320N + 588 , (D.16) 
Pi2(iV) = A^^ + 7A^^ + 53A^^ + 1887V2 + MiAT + 54 , (D.17) 
p^3(iV) = 7Ar^ + 70Ar^ + 730A^^ + 3340Ar^ + 6127Ar=^ + 5074A^2 + 2340Ar 

+504 , (D.18) 
PuiN) = 53A^^ + 665A^^ + 7044A^6 + 36078A^^ + 90789A^^ + 122301A^^ 

+93322 AT^ + 40780A^ + 7752 . (D. 19) 



-A3), El 





64(A^- 1) 




9(A^ + 1) 



1 

+- 



A^(A^ + 1) 

32,,., , , 16(Ar-l)(15Ar3 + 47Ar2 + 44Ar+18) 
--(4Ar + 5)S, + 27Ar(Ar + 1)^ 

8(8A^3 + 29A^2 + _ gg) 



9(Ar + 2)(A^ + 3) 



-Si 



8(A^-i; 



3(Ar + i) 

8{32N^ + 157A^3 + g^oiv^ + 71 AT - 66) 

9(Ar + l)(Ar + 2)(Ar + 3) ^ 

8Pi5(iV) 



C2 



g^isl^vj ^ 4Pi6(Ar) 

"27Ar(A^ + l)(Ar + 2)(Ar + 3) ^ 81Ar2(Ar + l)3(Ar + 2)2(Ar + 3) 

4(16Ar3 + 37Ar2 + 25Ar - 258) 4 
27(iV + 2)(Ar + 3) - 3(4^ + 5)C2-5i 



72 



with 



4(16A^=^ + 37 A^^ + 2bN - 258) 56(A^ - 1) 

9(A^ + 2)(A^ + 3) ' ' ~ 9(A^ + 1) 

8(124Ar4 + 593Ar3 + ii42Ar2 + i99Ar- 258) 16 

27(iV + l)(Ar + 2)(Ar + 3) - y + 5)^2,1 

2Pi7(Ar) 2(Ar - 1) (15Ar3 + 47iV=^ + AAN + 18) 

^ 27N{N + 1) (iV + 2) (TV + 3) ^ ^ 9iV(iV + 1)2 

2Pi8(iV) ^ 2Pi9(iV) 

27iV(iV + l)2(iV + 2)(iV + 3) ^ 81iV(iV + l)2(iV + 2)2(Ar + 3) ^ 

243iV3(iV+ l)4(iV + 2)3(iV + 3) J ' ^ ■ > 

Pi5{N) = 118A^^ + 779A^^ + 1838A^^ + 835A^2 - 114A^ - 432 , (D.21) 
Pi6(A^) = 231A^'^ + 1670A^^ + 4858ArV 1994Ar6 _ 12383A^^ 

-23476A^^ - 23234Ar3 - 9924Ar2 ^ ^944^ ^ i296 , (D.22) 
Pi7(iV) = 236A^^ + 1105A^^ + 2428A^^ - 23A7N^ - 3198A^ - 1296 , (D.23) 
Pi8(A^) = 1052Ar6 + 6813Ar5 + 18221A^^ + 14457A^^ - 3865A^^ 

-8166iV - 2592 , (D.24) 
Pi9(A^) = 2830A^^ + 24935iV^ + 87817iV^ + 134237Ar^ + 95737Ar3 

+28280iV2 - 29532A^ - 19008 , (D.25) 
P2o{N) = 5487A^^2 + 36857A^" + 58060iV^° - 225576iV^ - 995433iV^ 

-1617615A^^ - 1886730A^^ - 2199914iV^ - 1681960Ar^ 

-5448727V2 - 87696iV^ - 518407V - 15552 . (D.26) 



x(3),Fi 



{N + l){N + 2) 



512 512^ 
+ —Si 



+- 



64 2 64(47V2 + 9iv + 8) 
9 ^ 9iV 



9 ■ 9 

128(37V3 -N^ + SON + 12) 
27iV2 



^1 



128(3iV3 + 2iV2 + 17Ar + 6) 



1 

+ - 

e 



27N{N + 1) 



32 3 64(3A^=^ + 20N^ + 3N + 12) 2 



MP2i{N) 



+ y^'^' Y^'^' + 81iV(iV + 1)2(7V + 2) 
^64(4Ar3 + 87iV2 -85Ar 



81A^2(^^1)(^^2) 
64C2 



^1 



27A^2 

128(7V+ l)(Ar + 2) 
37V 



24) 256(3Ar2 + lOiV + 6) 

"^2 



27iV 



'S'2,1 



448 



73 



160^^ 64 4(48A^2 ^ lOlAT + 96) _,2 8{AN^ + 9N + 8) . ^ 

^3^1 + ~5~>J2,1>J1 '^2 ^Tivr ^2'-'2 



27 3 ' 9N ' 3N 

64 64(3iV2 + 7iV + 6) 8(56iV2 + 169iV + 112) 

+ y ^3,1 + ^ ^2,1,1 ^ ^4 

32(6A^3 ^ Qip^2 _ + 24) „ 16(3iV3 - A^^ + SON +12) 

32(6A''^ + 61N- - 21N + 24) 128(A^'^ + 9N- - ION - 6) 

+ 27iV^ 9m 

448 32(9A^^ - 623N^ + 894Ar + 276) 16P23(A^) ^.2 

^"9"^^ 8TiV2 ^ ~ 81Ar2(Ar+l)(7V + 2)^ 

16(3Ar3 + 2N^ + 17iV + 6) 16P24(iV) 

"C2 ~ 01 Ar9/ AT I T\r 1\T I ON '-'2 



9Ar(Ar + l) 81iV2(iV+ l)(iV + 2) 

32P25(iV) ^ 32P26(iV) , 



243iV2(iV + l)2(iV + 2)2 2437V(iV + l)3(iV + 2)2 



P2i{N) = 24A^^ + 196A^^ + 1047iV^ + 1748A^2 + 984A^+120 , (D.28) 

P22{N) = 33iV^ + 364A^^ + 1462iV^ + 2543iV^ + 1904iV + 132 , (D.29) 

P23(A^) = 48A^^ + 746iV^ + 2697iV^ + 2746iV2 + 1104A^ + 240 , (D.30) 

P24(iV) = 124A^^ + 198iV^- 2387A^^-6162A^2- 3632 A^- 480 , (D.31) 
P25(iV) = 264Ar^ + 4046A^^ + 21591Ar5 + 52844iV^ + 74856iV3 

+66812iV2 + 30576iV + 2640 , (D.32) 
P2q{N) = 363iV^ + 6758iV^ + 41285iV^ + 121235iV^ + 190235A^^ 

+150758Ar2 + 46964iV + 2904 . (D.33) 



'Gi _ T^UfCp f_ 1 224 1 
(TV + 1)(7V + 2) I ~£2 + - 



-645, 



16 (A^2 ^ 27A^ + 38) 8 (65 A^^ + 430A^^ + 881iV2 + 608A^ + 44) 
3(A^ + l)(A^ + 2) 3(A^ + l)2(Ar + 2)2 

8 (Ar3 + 30iV^ + 137Ar + 144) 
~ 3{N + l){N + 2){N + 3) ~ ^^^^ 
16(25Ar3 + ISOiV^ ^ 248A^ + 105) 4P27(Ar) 
^ 3(7V + l)(iV + 2)(iV + 3) ^ ^ 3N{N + 1)2(7V + 2)2(7V + 3) ' 



2P28(iV) 



3( AT + 1)3 (AT + 2)3 (AT + 3) 



(D.34) 



P27(Ar) = 13N^ + 373Ar5 ^ 2197Ar4 + 4907Ar3 + 4534Ar2 + 1464Ar + 144 , (D.35) 
P28(Ar) = 455A^^ + 5828Ar6 + 29700Ar5 + 77486Ar^ + 111073Ar3 

+87014Ar2 ^ 35364iv + 6840 . (D.36) 



74 



1 

+- 

e 



with 



16(iV=^ + 41iV-30) 2 
m{N + 2){N + ?,) ^ 

16 {IIN* + 850iV3 + 2707iV2 + 4160iV + 2724) 

27iV(iV + l)(iV + 2)2(iV + 3) ^ 

16P29(iV) 16 (13A^2 ^ _ gg) 

"3(iV + l)2(iV + 2f{N + 3) ^ 9iV(iV + 2)(iV + 3) ^ 

8(7iV^-73iV + 150) 3 40 
27iV(iV + 2)(Ar + 3) ^ 3Ar^2'5^ 

8 (7iV2 - 73A^ + 150) 16 (47iV2 + 667iV - 150) 160 

" m{N + 2){N + 2,) ' ' ^ 27Ar(7V + 2)(iV + 3) ^ ^ 37V 
4 (119A^^ - 962A^=^ - 3137A^2 _ 37847V - 3324) 2 16 

27A^(A^ + 1)(A^ + 2)2(Ar + 3) ^ ~ {N + 2) 

4 (433A^^ + 6746A^=* + 16769Ar2 ^ 131927V + 5052) 

27Ar(iV + 1)(A^ + 2)2(Ar + 3) ^ 

4P3o(iV) ^ 

81Ar(Ar + l)2(Ar + 2)3(Ar + 3) ^ 

9( AT + 1)3 (AT + 2)4 (AT + 3) J ' ^ ^ ^ 



P29(A^) = 7A^^ + IN^ - 203A^^ - 955 A^^ - 1604 - 948 , (D.38) 
P3o(Ar) = 427A^^ - 10201A^^ - 91517A^^ - 319471A^^ - 556466A^^ (D.39) 

-463084A^ - 137784 , 
P^i{N) = 253A^^ + 1884A^^ + 3444Ar5 - 9582A^^ - 50265A^=^ - 82446A^^ 

-58120A^ - 13680 . (D.40) 



{N + l){N + 2) 



+- 



32 2 



32 128(Ar+l) 



3A^ 



Si + 



256 



3 ^ ^ ^ 3 9A^ 9A^ ^ 



64 (20A^2 ^ 377V + 23) _ 2560 

^ 9N(nTT) ^ 9~ 

14 28 14 
+ + ^-^^^i' + ^^2' - 645'2,i,i + 3252,1^1 - 4C2^2 + + 32^3,1 



75 



16(117V-5)^3 16(iV + l) 



+-^>S'3S'i — 285*4 ^5"^ — ^C25'i ^^~;7T7 -S2S1 



16(llA^-5) 



N 



32(197V + 5) 64(iV+l) 

H >J3 77 '-'2,1 + 



16 (179iV^ + 200iV + 75) ^2 
9A^ ~" N ' 27N{N+1) 



3A^ 

Sf + 32C2 



16 (23A^2 + 2iV - 75) 32 {27AN^ + 816N^ + 921iV + 325) 16960 

"'-'2 r.^TA^/A^ , -i \ O '^'l + 



27iV(iV+ 1) 



27N{N+iy 



27 



64 64 (A^ + 4) 
y ^ 9(A^ + 2) 



+ 



16 2^64(7V+1) 16(67V + 5) 



32 (Ar3 + 33A^2 ^ 74^ ^ 33 (jtv^ + 25A^2 ^ gg^v + 80) 

H . - V / . "S^l 



27A^(A^ + l)(A^ + 2) 



27(A^+ l)(A^ + 2)2 



+- 



40^3 , 8^ ^ , 8(18iV+13) ^ ^ , 80(iV + l) ^ 16(54iV + 49) 
"27"^^ + 3^^"^^ + 9 + 3 + 27 



qor/v^lW 16(3iV+^ 32 + 3N - 14) ^ 
-32(7V + l)5_2,i 3 ^2,1 5-2 

16 (24A^^ + 38A^=^ + 15N'^ + SIN + 18) 16 (22Ar3 + qqn"^ + 357V - 1 



27A^(A^ + l)(A^ + 2) 



^2 + 



8(A^ + 4) 



C2- 



16P32(A^) 



3(A^ + 2)^ 81A^(A^ + l)2(Ar + 2)2 
13^4 8(9iV+ll) 



Si + 32{N + l)S_2Si + 



27N{N + 1){N + 2) 



Si 



81(N +1)^(N + 2f 



27*^^ 



52,5.-48(iV+l)5_2.5, + ^il^?^535, 



56 



+40(iV + l)5_35i-— Ca^i 

y 



27 

25^51 + 8(7V+1)C,5_. 



^'J^CA + 8(iV + DS.A + 24(iV + - i5<H^ 



'S'3,1 



+48(iV + + + 5. - 40(iV + 1,S_3,: 



9 



+32(A^ + 1)52,_2 + 



8 (65iV3 + 174iV2 + 31iV - 90) 
81iV(iV + l)(iV + 2) 



8 {72N^ + 34iV3 - 273Ar2 - 1577V + 90) 4 (iV^ + 33iV2 + 74iV + 36) 

-O2O1 H ^ ,r/ AT , \ / AT , ^\ C2'-'l 



27iV(Ar + l)(iV + 2) 

16 (8iV2 + 3iV - 14) 8(iV + l) 

— <b_2<Ji <b_4 



3(Ar + 2) 

56(iV + 4) 2(54A^ + 41) ^ 8 (432Ar^ + 707A^'^ + lllA^^ ^ ^54^ ^ ;^gQ) 
^ 9(A^ + 2) ^ 9 81Ar(A^ + l)(Ar + 2) 

2 2 16 (8A^- + - 14) 8(21A'i + 35A-^ + 33A^- + 106A' + 72) 

C2>-'i + ; ; J — 2 1 + n ^ '-'2 1 

3^ 1 3(Ar + 2) ' '^^'^^^ , , o^ 



9iV(iV + l)(Ar + 2) 
40 (8iV2 + 3N- 14) 

9{N + 2) ' 



9N{N + 1){N + 2) 



4P34(iV) 



81Ar(Ar + l)2(Ar + 2) 



.q2 



4 (7iV3 ^ 25iV2 ^ 807V + 80) 
9(Ar + l)(Ar + 2)2 



C2 



76 



64 (13A^^ + 36A^3 ^ 73^2 ^ ;^53^ ^ ;^3o) 

H ^^/». . . ^-.n '-'-2 + 



4P35(A^) 



27(A^ + l)(A^ + 2)2 

8P36(A^) 



81N{N + iy{N + 2y 



2A3N{N + lf{N + 2y 



2A3{N+lf{N + 2) 



with 



PssiN) = 

P34{N) = 

P35{N) = 

PseiN) = 



173iV^ + 1494iV^ + 4706iV3 + 7191A^2 + 5588A^ + 1872 , 
320iV^ + 2720iV^ + 9797iV3 + 17507A^2 + 15272iV + 5216 , 
700iV^ + 4572iV^ + 10909iV3 + 12447iV2 + 7624A^ + 2448 , 
312iV^ + 962iV^ + 456iV^ - 2947iV3 - 7203iV2 - 6880iV 
-2448 , 

3019A^^ + 33159iV^ + 152637iV^ + 388497iV^ + 593193iV3 
+540399iV=^ + 267080A^ + 52848 , 

5041Ar^ + 58633A^^ + 293478Ar^ + 808539A^^ + 1321572Ar=^ 
+1281873A^^ + 683576A^ + 154784 . 



(D.41) 

(D.42) 
(D.43) 
(D.44) 

(D.45) 

(D.46) 

(D.47) 



x(3),Jla 



Tl-njC'A J 1 



Ar2(Ar+l)2 £3 



32(A^-3)(2^- + A^-4) 



9 



+ - 



16P38(iV) 



27(A^-l)Ar(7V+l)(Ar + 2 
16 



{N - 3)(2A^' + N- ^)S2 + -{N - 3)(2Ar2 ^ ^ _ 4^^^ 



8P39(iV) 



81(iV - iym{N + l)2(iV + 2f 
28 

- — {N - 3){2N^ + N- 4)C3 



9 



+ 8{N -3){2N^ + N -4:)Ss 
9{N -1)N{N +1){N + 2) ^ 



2P38(iV) 



C2 



9{N -1)N{N +1){N + 2) 
4P4o(iV) 

"243(iv - i)m^N + if{N + 2)3 r ' 



(D.48) 



with 



PssiN) 

P39{N) 

P4o{N) 



16A^^ - 17A^^ - 209A^^ + 142A^^ - 65A^^ - 629A^^ + 186A^ 

+216 , (D.49) 

68A^^^ + 75A^^° - 1715A^'^ - 2112 A^^ - 3360A^^ + 24A^^ 

+24544Ar5 + 11469Ar^ - 24037A^^ - 6432Ar2 + 4932iV + 3024 , (D.50) 

460Ar^^ + 1799A^^^ - 11123iV^^ - 40963iV^^ - 38795iVi^ 

+155048A^^° + 356978iV^ - 365210iV^ - 953263iV^ + 458455iV^ 

+924557A^^ - 305717Ar^ - 358086iV3 - 66708Ar2 ^ 870487V 

+38880 . (D.51) 



77 



X(3),J16 



TpUfCA 



A^2(Ar+l)2 £3 



1 

+ - 



1 

+- 

s 



9 

32P4i(iV) 



9 ^ ^ ^ 27Ar(Ar + l)(iV + 2) 

-^(4Ar2 + 4Ar - 5)>5i2 - ^(4A^' + 4A^ - 5),52 



16P4i(iV) 



^1 



16P42(A^) 



81iV2(iV + l)2(7V + 2)2 

[N^ + 47\ 
4P4i(iV) 



- —{4N^ + AN-5)Sl 



44n ^fi 

■^(47V2 + 47V - 5)^3 + y (47V2 + 47V - 5)C3 



^1^ + 



52P4i(iV) 



27A^(A^ + l)(A^ + 2) ' 27A^(Ar + l)(Ar + 2) 



4P4i(iV) 



9N{N + 1)(N + 2y 9 



^2 

C2-—(4Ar2 + 4iV- 5)52^1 



8P42(iV) 



i(47V^ + 4Ar-5)C25,- si^2(;v + l)2(iV + 2)2 

8P43(iV) 



^1 



243A^3(7v + l)3(Ar + 2)^ 



(D.52) 



P4i(iV) = 447V5 + 191iV^ + 99iV3 - 17iV2 + 136iV + 60 , (D.53) 
P42(Ar) = 232Ar8 + 1978A^^ + 4292A^6 + 4447A^^ + 3446A^^ - 1058A^=* 

-5180A^^ - 2712Ar - 720 , (D.54) 
P43(iV) = 1328A^^^ + 16796A^^° + 66200Ar'^ + 134952A^^ + 148833A^^ 
+40020A^^ - 52496A^^ + 66940A^^ + 180160A^^ + 115296A^^ 

+45504A^ + 8640 . (D.55) 



1 ?,2{N^ - 10A^2 + 43A^ + 46) 



(A^ - 1)A^(A^ + l)2(Ar + 2) I £3 

16P44(iV) 



1 



1 

+- 

e 



27(N -1)N{N + 1){N + 2) 
16 



— {N^ - lON^ + A3N + 46)^2 - -{N^ - WN^ + 43A^ + 46)C2 



8P45(iV) 



81(A^ - l)2Ar2(Ar + l)2(iv + 2y 



78 



-8(A^^ - lON^ + 43N + 46)5'3 + y (iV^ - lOA^^ + 43A^ + 46)C3 

9{N - 1)N{N + 1){N + 2) 9(7V - 1)N{N + 1){N + 2)^^ 

4P46(iV) 

243(Ar - lfm{N + lf{N + 2)3 

with 

Pu{N) = llA^^ - lOON^ + 317N^ - 59N^ - 218AN^ - 765N^ + 106AN 

+276 , (D.57) 

P45(Ar) = Q7N^^ - 606iVi° + 221iV^ - 5084A^* - TTOOAT^ + A1990N^ 

+48735iV^ - 45550iV^ - 33703iV^ + 18682iV2 + 7212iV + 1656 , (D.58) 

Pie{N) = 431A^^^ - 3320iV^^ - 7624A^i3 - 6288Ar^2 + 1QU83N^^ 

+372478A^^° - 690814iV9 - 1786393A^^ + 805264iV^ + 2328180Ar6 

-8516727V^ - 1227601iV^ + 257348Ar3 + 1237927V2 + 48240Ar 

+9936 . (D.59) 



(D.56) 



i(3),Xii, 



{N -1)N{N + iy{N + 2) 
32(3iV2 - 23iV - 20) ^ 



64(3A^2 _ 237V - 20) 
9 



32(547V^ - 127 - 289N^ + 322Ar + 376) 
27{N + 1){N + 2) 



1 

+- 

e 



8{3N^ - 23N - 20) 2 104(3A^2 _ 237V - 20) 
9 + 9 ■ 



C2 



16P47(A^) 



^8(3A^2-23A^-20) 

16(54iV^ - 127A^3 _ 289A^2 ^ 322iV + 376) 

27{N + 1){N + 2) ^ ^ 81(7V + l)2(Ar + 2)^ 

4(3A^2 - 23Ar - 20) 3 

+ 27 

440(3iV2 - 23A^ - 20) 56(3A^2 _ 237V - 20) 
+ 27 9 

"9 



27 ^ g 

52(3A^2 _ 237V - 20) 4(3A^2 _ 237V - 20) 
6201 H C2'-'l 



79 



4(54A^^ - 127N^ - 289N^ + 322N + 376) ^ 
27(A^+ l)(A^ + 2) " 

52(54A^4 - 127A^=^ - 289N^ + 322N + 376) 
27{N +1){N + 2) 

4(54Ar4 - 127A^3 _ ggQiV^ + 322Ar + 376) 

9(7V+ l)(7V + 2) ^ 

8P48(iV) c;, 
81(iV+ l)2(iV + 2)2 ' 

8P49(iV) 

243(iV + l)3(iV + 2)3 



(D.60) 



with 



P4y{N) = 468Ar6 + 574A^^ + 199A^^ + 1615A^^ - 556AN^ - 18092A^ 

-10304 , (D.61) 

P^siN) = A68N^ + 57AN^ + 199N^ + 161bN^ - bb6AN^ - 18092 AT 

-10304 , (D.62) 

P49{N) = 3096A^^ + 11084iV^ + 7658A^^ - 72687iV^ - 301881Ar^ 

-4124347V3 ^ 102607V2 + 4332407V + 231808 . (D.63) 



x{3),Lla 



+ ^^47V+5W I 16(27V3 + 57V2 + 6jV + 4) 



1 



8P5o(iV) 



277V2(7\r+ l)2(7\r + 2) 



14 2 

-{m + 5)S! + -{AN + 5)C2Si 



1 4 28 8 

+-{4N + 5)S2Si + -{AN + 5)Ss + -{AN + 5)^2,1 

2(36A^^ + 106A^2 ^ _ 4) 2(2Ar3 + 5N^ + 6N + A) 



9A^(A^ + 1) 
2(84A^=^ + 130A^2 ^ _ 52) 



9A^(A^ + 1) 

2P52(A^) 



^2 + 



3A^(A^ + 1) 
4P5i(A^) 



S1N^{N + 1)^{N + 2y 



S1N^{N + 1)^{N + 2) 



Si 



K in 

+ -{AN + 5)52^? + —{AN + 5)^3^1 + 4(4Ar + 5)S2,iSi 
6 9 



80 



+g(4iV + 5)S', + ^^^(2-^2 + y (4iV + 5)^4 

+4(47V + 5)^3,1 - 4(47V + 5)^2,1,1 - ^ g^^(^ ^ 'Sf 



9N{N + 1) ^ 9 

(404A^=^ + 834A^2 ^ gQjjv - 156) 14(2Ar3 + siv^ + 6iV + 4) 

S2S1 nnT/nT I -i\ Ca 



27A^(A^ + 1) 9A^(A^ + 1) 

2(908A^3 + 1086A^2 _ ^97^ _ ggQ) 4(20A^2 ^ gg^v + 61) 



81A^(A^ + 1) 9(A^ + 1) 



'S'2a 



27N^{N + 1)^{N + 2) ' 9N^{N + iy{N + 2) 

, ^54 (iV) 2P55(iV) 

27Ar2(Ar + i)2(iv + 2) ^ 243Ar3(Ar + 1)3(7V + 2)2 ^ 

P56(iV) 



486iV4(Ar + l)4(iV + 2)3 



(D.64) 



with 



P5o{N) = A3N^ + 301N^ + 587N^ + ASAN^ + 132N^ -SON -A8 , (D.65) 
P5^(Ar) = 190A^^ + 870A^^ + 802Ar^ - 543A^3 - 938Ar2 - 192iV + 144 , ^ gg^j 
P52(A^) = 991iV^ + 9415iV^ + 32215iV^ + 52113iV^ + 43252iV^ 

+17484A^^ + 2928iV'^ + 3104iV2 ^ S648N + 1152 , (D.67) 
P53(iV) = 276A^^ + 1472iV^ + 1976A^^ + 425iV3 - 6747V2 _ 3527V 

+48 , (D.68) 
Pr,i{N) = 948iV^ + 4328iV^ + 6176iV^ + 2945iV3-2iV2 + 608iV 

+624 , (D.69) 
P55(iV) = 1379Ar^ + 9765Ar^ + 20193iV^ + 10203A^^ - 15291A^^ 

-30216A^^ - 27620A^3 - 14160A^^ - 5472A^ - 1728 , (D.70) 
P56{N) = 21055A^^^ + 262474A^^^ + 1324746A^^° + 3564860A^'^ 

+5666307Ar« + 5550126A^^ + 3387956Ar^ + 1293128A^^ 

+273984Ar^ - 120320A^3 - 219648A^^ - 129024A^ - 27648 . (D.71) 



J(3),i^i6 



1 



64 16(2iV3 + 7Ar2 + 6Ar + 3) 
"9" 9A^(iV + l)2 



64^2 



9 ^ 



32(5A^ + 14) 
27(iV + 1) 



^1 



8P57(iV) 



27N^{N +lf{N + 2) 



81 



with 



32o3 , 8 , 64 208 64 16(57V + 14) 



2{2N^ + 7N'^ + 6N + 3) 8(74Ar3 ^ i21N^ + 38N - 27) 

-C2 - 



3A^(A'+1)2 27N{N + iy 

8(A7N^ + 13N^ - 196N - 108) 2P5s{N) 



-Si + 



81Ar(iV + l)2 S1N^{N +iy\N + 2y 

8Sf 56 4 2 32 2 208 64 

+"27 9"^^^ 3^^^ "9" "27" "9" 2,1*^1 

80 2 8 128 32 64 



16(5A^ + 14) 3 4(5iV + 14) 



81(A^+1) ' 9(Ar + l 
32(5iV + 14) 14(2iV3 + 7N^ + 6N + 3) 

-S2O1 AAr/Ar , i\9 C3 



27(A^ + 1) 9A^(A^ + 1) 

4(616A^=^ + 899A^2 + 202A^ - 243) 32(5A^ + 14) 



81A^(A^ + 1)2 27(N+1) 
4(47A^3 + 13A^2 _ igg^ _ 108) PstCA^) 



'S'2,1 



81Ar(A^ + l)2 1 9Ar2(Ar + l)3(7V + 2)^^ 

, 4P59(iV) ^ 2P,o{N) 

81Ar2(iV+i)3(iv + 2) ' 2437V2(^+i)3'^i 

P6i(iV) 



486Ar4(Ar + l)5(Ar + 2)3 



(D.72) 



P57(iV) = 43Ar6 + 302iV^ + 566iV^ + 445iV^ + 129Ar2 - 81iV - 54 , (D.73) 
P^s{N) = 1207A^^ + 11543iV^ + 40819iV^ + 70919iV^ + 65266iV^ 

+27608A^^ + 312iV-'^ + 162iV^ + 3672iV + 1512 , (D.74) 
P^g{N) = 7A2N^ + 3260iV^ + 5157iV^ + 3025iV2 + 251A^^ + 513Ar 

+486 , (D.75) 
Pgo(iV) = 323iV^ - 3972iV^-9291iV3- 4456 Ar2- 1080A^ + 648 , (D.76) 
PeiiN) = 25807A^^2 + 324932iV" + 1659342iV^° + 4520784iV^ 

+7180599iV^ + 6692496iV^ + 3387488iV^ + 763528Ar5 

+119892A^^ + 17388A^=^ - 147960A^2 - 143856A^ - 38880 . (D.77) 



x(3),Mi 
^Q9 



Ar2(iV+ 1)2 



16(4iV2 _ _ 9) 



9 



1 6 

— {N^-2N-2)Si + 
9 



8P62(iV) 



27N{N + 1){N + 2) 



82 



1 

+- 

£ 



Un^ -2N- 2)Sl + ^{5N- 14) (SAT + 4)^2 
9 9 

SlN'iN + iy{N + 2) J + 27^^^ 

4 14 

— (109Ar2 - 2187V - 245)5'3 {4N^ - 8N - 9)C3 

27 9 

26 2 

— (iV^ - 2Ar - 2)>S2-5i + -(iV" - 2Ar - 2)C2-Si 

2^63 (iV) , 2P65(iV) 

27A^(iV + 1) (iV + 2) ^ 277V(iV + 1) (iV + 2) ^ 

, n2(iV) , 4P66(iV) c. 

9A^(A^ + 1) (AT + 2) 81 A^2(^ + 1)2(^ + 2)2 ' 

I 2P67(iV) 1 

243A^3(7v + l)3(Ar + 2)3j ' ^ ^ ^ 



+ 
+ 



Pg2(iV) = 12iV6 - 17iV^ - 155iV^ - 43iV''^ + SaiV^ - 207A^ - 138 , (D.79) 
Pq3{N) = 17N^ + 59N^ + 21N'^ -38N^ + 16N + 2A , (D.80) 
P64{N) = 24A^^ - 91A^^ - 713A^^ - 558Ar5 - 424A^^ - 1202A^^ + 912A^2 

+2055A^ + 918, (D.81) 
Pg5(iV) = 72Ar6 _ ii9Ar5 _ 989Ar4 _ 279A^^ + 536A^^ - 1258A^ - 852 , ^ g2) 
PmiN) = 88A^^ + 553A^^ + 977A^^ + 241A^^ - 436A^^ - 26A^^ - 296A^^ 

-624Ar - 288 , (D.83) 
P67{N) = 26AN^^ + 1381A^^^ - 4166A^^° - 47256A^'' - 152187A^^ 

-251073A^^ - 189276A^^ + 16418A^^ + 68699A^^ - 79947 A^^ 

-148962A^2 - 95364A^ - 23112 . (D.84) 



i(3),Afla 



N{N +1){N + 2) I £3 
32(Ar3 + 2Ar2 + 1 



32 



(2A^^ + 3A^ + 2)Si 



9(A^ + 1) 



+ 



9 



(6iV2 + lliV + 8) Si 



■y (iV2 - 4)S_2 - ^(14A^^ + 33A^ + 12)^2 



16P68(iV) 



-^1 + 



16P69(iV) 



277V(7V + l)(7V + 2) 277V(7V + l)2(Ar + 2) 



83 



1 

+- 

e 



32 4 
_ (Ar2 _ 4)5_25i - -(6Ar2 ^ ^ _ 44)5^5^ 

-y (A" - 4)S_3 + f (A" - 4)5-2,1 + ■^iV(llJV + 16)S2,i 

2T ' ^ 21N{N + 1){N + 2) ^ 3(iV + 1) 

U{1N^ - 17 - 43A^2 + + 6) ^PrsiN) 

^ 27N{N + 1) 27A^(A^ + l)(A^ + 2) ^ 

8Pn(iV) 8P72{N) 



81A^2(7v + l)2(Ar + 2)2 ' Sl^iN + 1)%N + 2y 

+^(307V2 + 59 AT + U)St + -(GTV^ + 117V + 8) (2^1 
54 3 

-8(A^^ - 4),5_25i + -{2N'^ + 63N + 188)52^1 
9 

28 40 
__(27V2 + 37V + 2)C3-Si - y (iV' - 4)^-3-5i 

-^(IMTV^ ^ _ 368)5'35i + 16{N^ - 4)5_2,i5i 

+^(25A^^ + 32A^ - 12),52,i-5i + —{6N^ - 205N - 292)S^ 
9 ' 18 

+ ^(iV2-4)5_4-^(iV2-4)C2^_2 

-^(147V2 ^ 33^ ^ ;^2)C252 - ^(^' - 4)-5_2-S2 

7 4n 

--(54iV2 + 99iV - 4)^4 +^{N^- 4)-S_3,i 

9 3 

39 Ifi 
-^(N' - 4)^2,-2 + y (77V' + SN- 6)^3,1 

-16(7V2 - 4)5_2,i,i - ^(257V2 + 32iV - 12)^2,1,1 

2P74(iV) 3 _ 2P68(iV) ^ ^ 

81iV(iV + l)(iV + 2) ^ 9Ar(Ar + l)(Ar + 2)^' ^ 
32(77V^ - 17iV3 - 43iV2 + 57V + 6) 2P75(iV) 
^ 9iV(iV + l) ' ^ 27iV(ivTT)(iVT2) ' ' 

28(iV3 + 2N^ + 1) 80(7iV^ - 17N^ - A3N^ + 5iV + 6) _ 

' "Cs H o^TAr/Ar I IN '-'-3 



9(iV+l) 27iV(iV + l) 

4P76(iV) 32(7A^^ - 17A^=^ - 43A^2 ^ 5A^ + 6) 

"'-'3 „Ar/Ar , i\ '-'-2,1 



81A^(Ar + l)(A^ + 2) 9N{N + 1) 

:'^2,l + 01 Ar2/Ar , iN2/Ar , r,\7^l 



27N{N + 1){N + 2) ' 81Ar2(Ar + l)2(Ar + 2) 

2P69(iV) ^ 32P79(iV) 

rS.2 - 01 Ar9/ AT , 1\9/ AT , r.\ '-'- 



9Ar(Ar+ l)2(iV + 2)^ 81Ar2(Ar+ l)2(Ar + 2) 



84 



2i^8o(A^) . 4P8i(7V) 

'-'2 — r>joAr-J/Ar , -i\Q/Ar , r.Nci'-'l 



with 



81A^2(7v + l)2(Ar + 2)2 ^ 243A^3(7v + l)3(Ar + 2)3 

4P82(iV) 

'243A^3(7v + l)4(Ar + 2)3 



P68(^^) = 16A^^ + 72A^^ + 75A^=* - ISSiV^ - 232Ar - 48 , (D.86) 

PmiN) = 14A^^ + 67A^^ + 85A^^ + 37A^^ - 18A^^ - 47A^ - 6 , (D.87) 

P7o(iV) = 42A^5 + 217A^^ + 316A^^ - 190A^2 - 532A^ - 144 , (D.88) 
P71 (iV) = 140A^^ + 1131 A^^ + 2545A^^ - 258A^^ - 6077 A^^ - 4596A^^ 

+476A^2 + 816A^ + 288 , (D.89) 
P72{N) = 154A^'^ + 1208A^^ + 3481A^^ + 4866A^^ + 3091A^^ + 612Ar4 

+975A^^ + 1549A^^ + 336A^ + 36 , (D.90) 

P73(Ar) = 82A^^ + 417A^^ + 686A^^ + 210A^^ - 348A^ - 144 , (D.91) 

P74{N) = 94Ar5 ^ 5Q7^4 ^ ygg^s _ 2947V2 _ ii32Ar - 336 , (D.92) 

P^^{N) = 90A^^ - 283A^^ - 1618A^=^ - 410A^2 + 1900A^ + 432 , (D.93) 

P7e{N) = 716A^^ + 2238A^^ + 939A^^ - 2532Ar2 _ 2252Ar - 528 , (D.94) 

P77{N) = 47A^^ + 136A^^ + lOA^^ - 295A^^ - 252Ar - 60 , (D.95) 
P7s{N) = 336Ar^ + 2909A^^ + 7295A^^ + 2096A^^ - 11749A^^ - 9748A^^ 

+1580iV2 + 2448Ar + 864 , (D.96) 

P7g{N) = 50A^^ - A^^ - 70A^^ + 530A^^ + 716A^^ - 97A^2 - 84A^ - 36 , (D.97) 
P8o(Ar) = 344Ar^ + 2775A^^ + 9197 A^^ + 13188A^^ + 6341A^^ + 1560A^^ 

+4404A^2 ^ 2880Ar + 864 , (D.98) 
Psi{N) = 1252A^" + 14118A^^° + 58795A^^ + 122580A^^ + 153807A^^ 
+ 131190A^^ + 43487Ar5 - 64116Ar^ - 66172Ar3 - 14880Ar2 

-7488A^ - 1728 , (D.99) 
P82(Ar) = 1550iV^^ + 16786A^^^ + 74226A^^° + 174527A^^ + 224721 A^® 
+133572A^^ - 9824A^^ - 69685Ar^ - 68916Ar^ - 46499A^2 

-12426Ar2 - 2340A^ - 216 . (D.lOO) 



i(3),iVii, 



N(N + 1)(N + 2) 1 £3 
32Ar(A^2 + 2Ar - 2 



9(Ar + 1) 



16 

— (4Ar2 + 5Ar + 2)5i 
9 



4 1 28 

-(8Ar2 ^ 5^ _ 2)^2 _ -1){N + 2)S^2 



■^{32N^ + 71N +18)S2 
9 



8P83(iV) 



^1 + 



16P84(A^) 



27A^(A^+ l)(Ar + 2) 27(A^+ l)2(Ar + 2) 



85 



1 

+- 

e 



|-(16Ar2 + 5iV - 10)^3 ^ 2 ^4^2 + SAT + 2)C25i 



-^(AT - l)(Ar + 2)5_25i - ^(48Ar2 + 59Ar - 118)^2^1 
160 64 

-^(iV - l)(iV + 2)^_3 + y (iV - l)(iV + 2)5_2,1 

+^(28A^2 + 31 - 26)^2,1 - ^(284Ar2 + 4577V - 158)^3 
2P85(iV) ^2 47V(jV^ + 2Ar-2) ^ 

o/Af , i\ ^2 



27Ar(Ar + l)(Ar + 2) ' 3(A^ + 1) 

64(11A^3 _ i^^2 _ 59^ _ ;^Q) 2P86(A^) 

-0-9. + . , , , ^ . , . rT>-'2 



27(iV + l) 27iV(iV + l)(iV + 2) 

4P87(iV) 8P88(iV) 

+ 0-. AT/AT , -.NO/AT , ^NO '-'I 



81A^(A^ + l)2(Ar + 2)2 81(Ar+ l)3(Ar + 2)2 
+ ^(32Ar2 + - 26) 5f + ^(8Ar2 + 57V - 2)C252 

-(AT - l)(iV + 2)S^2Sl + — (-224Ar2 - 251Ar + Am)S2Sl 



14 320 
__(47V2 + 57V + 2)C3^i - —{N - 1){N + 2)3.^3^ 

2 1 28 

- — (736A^2 + 763A^ - 1510),53,5i + {N - 1){N + 2)5_2,i,5i 

27 3 

+ ^(88iV2 + 79A^ - 134)52,1^1 + ^{12^"^ - 403iV - 1034)^2^ 
9 ' 3d 

+^(7V - l)(Ar + 2)5_4 - y (A^ -1){N + 2)C25_2 

1 64 
+ (_32Ar2 _ 7;Liv - 18)C2-52 - —(AT - l)(Ar + 2)5_2-S2 

D 9 

64 1 

- — (A^ - 1)(A^ + 2),52,_2 + — (-1120A^^ - 1903A^ + 526)^4 

3 ' 18 

+ ^{N - l)(Ar + 2)5_3,i + ^(48Ar2 + 45A^ - 82)^3,1 

-32(Ar - l)(Ar + 2)5_2,i,i - ^(68Ar2 ^ _ 142)^2,1,1 
P89(iV) ^3 PUN) ^ 



81Ar(A^ + l)(A^ + 2) ' 9Ar(A^ + l)(A^ + 2) 
128(11A^=^ - 14A^2 _ 59^ - 10) ^ ^ P9o(A^) 



o_2<Ji + ^_ - - , - ^ - , - rr'-'2'-'l 



27(Ar + l) 27A^(A^ + l)(Ar + 2) 

28Ar(A^2^2Ar-2) 2P9i(Ar) 

9(A^+1) 81Ar(A^ + l)(A^ + 2) ^ 

80(11A^=^ - 14Ar2 - 59A^ - 10) 4P92(Ar) 

^ 27(iVTl) ~ 27iV(iVTl)(iVT2) ^'^ 

32(llAr3_i4^2_59^_10)^ Pg3(iV) 

-2,1 + 01 Ar/Ar , ^\'}JiKT , oN^-^l 



9(Ar + i) 8lAr(Ar + l)2(iV + 2)2 



86 



with 



Ps3iN) 

P84(iV) 

Ps7{N) 
PssiN) 

P89{N) 

P9o{N) 

P9l{N) 
P92{N) 

P9s(N) 

P94{N) 
P95{N) 



:C2- 



128P94(A^) 



9(N + iy{N + 2y" 81(N + iy(N + 2) 

P95{N) ^ 2Pg,{N) 



S1N{N + iy{N + 2)^ 

4P97(iV) 

"243(A^ + l)4(Ar + 2)3 



2A3N{N + lf{N + 2y 



= UA8N^ + 14464iV^ + 45664A^^ + 27660A^^ - 1573 lliV^ 
465972A^^ - 638594iV-'^ - 483184iV2 - 176200A^ - 21120 



P96{N) 

Pgr{N) = 824iV^ + 7474iV^ + 20726A^^ - 2277iV^ - 140187Ar 



-335448Ar^ - 348968iV2 - 155743Ar2 - 258747V - 11616 . 



(D.lOl) 



20N^ + 76N^ + 35N^ - 256N^ - 3WN - 96 , 
8A^^ + 43A^^ + 45A^^ + 29N^ + 22N - 24 , 
40A^^ + 94A^^ - 155A^=^ - 856A^2 - 836A^ - 192 , 
184Ar5 + 622A^^ + 251A^=^ - 1240A^^ - 1596A^ - 384 , 
184A^^ + 1388A^^ + 2586A^^ - 2671A^^ - 12540A^^ 
-13444Ar2 _ 6160A^ - 960 , 

88A^^ + 530 A^^ + 227A^^ - 3946A^^ - 11473A^^ - 14065 A^^ 
-7202A^ - 528 , 

80A^^ + 130A^^ - 535A^^ - 2056A^^ - 1876A^ - 384 , 
624Ar5 ^ 454^4 _ 4951^3 _ yqqON'^ - SAAN + 384 , 

2248A^^ + 4052 A^^ - 9473A^^ - 22736A^^ - 11468A^ - 2208 , 
236A^^ + 374A^^ - 1269Ar=^ - 3032Ar2 _ i628A^ - 288 , 
368A^^ + 2498A^^ + 3042A^^ - 12661A^^ - 37056A^^ 
-35356A^2 - 14144iV - 1920 , 
19A^^ + 237V^ + N^ + 131N^ + 1787V - 28 , 
13767V^ + 7298 A^^ + 11630A^^ - 8683A^^ - 47392A^^ 
-46660A^^ - 15168A^ - 3840 , 



x(3),Oi 



1 



(Ar + l)(Ar + 2)U3 



16(9Ar + 10) 128(2A^ + 3) 

^ ~ 9(Ar+ l)(Ar + 2) 



9Ar 



1 



4(23Ar + 26) 2 4(31Ar + 34) 
9iV 9iV 



8(96A^3 + 383Ar2 ^ 552^ + 268) 
^ 27Ar(Ar + l)(Ar + 2) ^ 

16(99A^^ + 562A^3 ^ io35Ar2 + 656A^ + 48) 
^ 27(Ar+ l)2(Ar + 2)2 



2(51A^ + 58) 
277V 



87 



with 



2{9N + 10) . ^ 2{67N + 74) 4(87Ar + 94) 



3iV ^ 9Ar " ' 27N 

8(23iV + 26) 16(2iV + 3) 

"'-'2,1 ~ O/AT I IN/AT I 



9N ' 3{N + l){N + 2) 

2(644iV3 + 2135iV2 + 2692iV + 1516) 
^ 27iV(iV + l)(iV + 2) ^ 

2(244Ar3 + i003iV2 + 14367V + 668) 2 4Pio2(iV) 
^ 27Ar(A^ + l)(A^ + 2) ^ ^ 81A^(Ar + l)2(Ar + 2)2 ^ 



4Pio3(iV) 



(107iV + 122) 4 
108A^ ^ 



81(A^+ l)3(iV + 2) 

(540iV-^ + 2243A^2 ^ 3204iV + 1468) 3 _ 4(23A^ + 26) 

^ 81A^(A^ + l)(A^ + 2) ^ 37V ^'^ 

, 4(61 TV + 70) ^ (1397V + 154) ^ (23iV + 26) 2 

+ ^iv 18iV 67V 

(7557V + 842) ^2 (317V + 34) 

367V 67V 
2(1797V+194) 4(617V + 70) 14(97V+n)) 

-D3O1 — — 'J2,l'-'l -\ TTTZ S3'-'l 



277V 9N ' 9N 

{511N + 562) 2(23 lOA"'^ + 7337A'2 ^ 8856A' + 5284) 

-'-'4 H 7TTTT 



18A^ 817V(Ar + l)(Ar + 2) 

836A^=^ + 3283A^2 + 4636 + 2300 112(2Ar + 3) 

^ 277V(7VTT)(7VT2) ' ' 9(7VTl)(7V + 2)^^ 

4(2147V3 + 8897V2 + 1298A^ + 596) 
^ 27Ar(Ar + l)(A^ + 2) ^'^ 

(967V3 + 3837V^ + 552Ar + 268) P98(A^) ^.2 

^ 97V(Ar+ l)(7V + 2) ' 81Ar(iV + l)2(7V + 2)2 ^ 

81A^(A^ + l)2(Ar + 2)2 ^ 
2(99A^^ + 562 A^=^ + 1035A^2 ^ ggg^ ^ 45) 
^ 9(A^ + l)2(Ar + 2)2 

2A3N{N + l)%N + 2)^^ 243{N + l)\N + 2y ^ ' ^' ' 



Pgs{N) = 560A^^ + 4964A^^ + 24041A^3 ^ ^^jqqn^ + 58136A^ + 22256 ,(D.118) 
P99(A^) = 5646A^^ + 27934Ar6 - 28488A^^ - 426969A^^ - 1064262Ar3 

-12028807V2 - 648456A^ - 133672 , (D.119) 
Pwo{N) = 66465A^^ + 706420A^^ + 3202098Ar6 + 8134152A^^ 

+12726405Ar^ + 12559980A^^ + 7477056Ar2 + 2214080A^ 

+120576 , (D.120) 
Pioi {N) = 49967V^ + 40516Ar^ + 140065A^^ + 238954Ar2 + l91872Ar 

+54736 , (D.121) 
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Pl02(iV) 



= i2N^ - 388N^ - 6171 A^^ - 18898A^^ - 21564Ar - 8416 , (D.122) 
= 2745A^^ + 20977A^^ + 59601A^^ + 76519A^^ + 40266A^^ 

+3700A^-480. (D.123) 



A, 



^(3),Pla 

Qg 



(Ar + l)(A^ + 2) I £3 



64(Ar + 4) 64(A^ - 4) 



9(A^ + 2) 



9A^ 



^1 



+- 



16(3A^ - 8) 



,2 , 16(57V + 16) 

97V 97V ' 



32 (SiV^ - 98iV2 - 251iV - 112) 32 (23A^2 



149iV2 + 277iV + 160) 



27Ar(A^ + l)(Ar + 2) 



27(iV+ l)(iV + 2)2 



+- 



8(7iV - 16) 3 8(7V-4) 8(3iV-32) 

"'-'l Win- S,2'-'l ;777 020i 



27N 



3N 



9N 



16(29Ar + 52) 128(iV-2) 

H ^ZTTT '-'3 TTTT '-'2,1 



27A^ 
8 {30N^ - 184A^2 



9N 

525N - 224) 2 8(iV + 4) 



27A^(A^ + l)(A^ + 2) 

(9 lA^'^ + 500A- + 797A^ + 1 18) 



^1 - ^)^^^C2 



27Ar(Ar + l)(A^ + 2) 



So 



81(Ar + l)2(Ar + 2)3 

- 64) 



- 321 



27Ar 



St 



3{N + 2) 

81N{N + 1)^{N + 2)^ 

2(3N - 8) 



.Si 



3N 



<2St 



9N 
128(A^-2) 



9A^ 



27A^ 



9A^ 
64(Ar - 2) 



S2AS1 + 



{9N + 320) _2 2(5A^ + 16) 



.53,1 + 



9N 
128(A^-2) 



-S' + 



3N 



C2S2 



3N 9N 
4 (74A^3 _ 355^2 _ io73Ar - 448 

81Ar(Ar + l)(Ar + 2) 
4 (8Ar3 _ 9gjY2 _ 25ijv - 112) 

9Ar(Ar + l)(A^ + 2) 
4 {6N^ + 8UN^ + 1893A^ + 896) 



2(113iV + 256) 

'-'2,1,1 -\ TTTl '-'4 



9A^ 



■s! 



C2'5l 



27Ar(Ar + l)(Ar + 2) 
8 (484Ar3 + 1778Ar2 + 2297A^ + 1456) 

81Ar(Ar + l)(Ar + 2) 
64 (11A^=^ - 43Ar2 - 137 A^ - 56) 



56(iV + 4) 

'-'2'-'l + ^/ AT , 0\ ^3 



9(Ar + 2) 



S 



2,1 



27Ar(Ar + l)(A^ + 2) 

IPumiN) ^2 , 4 (23A^^^ + 1 lOA-^ + 277A^ + 160) 



81Ar(Ar + l)2(Ar + 2)2 

4Pio7(A^) 
"81Ar(Ar + l)2(Ar + 2)2 



St 
S2 



9{N + 1){N + 2)^ 

8PMN) ^ 
^1 



C2 



243Ar(Ar + l)3(Ar + 2)3 
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with 



8Pio9(iV) 



243(A^+ l)3(iV + 2)4 



(D.124) 



P^oiiN) = 38A^^ + 1660A^^ + 6647 A^^ + 10051A^^ + 7013A^ + 2464 , (D.125) 

PiobiN) = 307A^^ + 2794Ar4 + 9631A^^ + 16150A^2 + 13219A^ + 4240 , (D.126) 

Pioe{N) = 24A^^ - 3140A^^ - 13401A^^ - 20075A^2 - 13587 - 4928 , (D.127) 
P^07{N) = 1052A^^ + 8260A^^ + 25625A^=^ + 40447A^2 + 32003A^ 

+9856 , (D.128) 
Pios{N) = 1210A^^ + 22502A^^ + 119904Ar5 + 319074A^^ + 502572Ar3 

+486540A^2 + 259511A^ + 49600 , (D.129) 
Pioq{N) = 3395A^^ + 40451A^^ + 201441Ar5 + 547869A^^ + 879225A^^ 

+832767A^2 + 431665A^ + 94720 . (D.130) 



j(3),Pl6 



N{N + l){N + 2)\ £3 



9 ^ ' 9(Ar + 2) 



-{UN + 30)^2 - - 30)-s2 + 



16(2A^3 _ 2577V2 - 542A^ - 172) 
27(A^ + l)(A^ + 2) 



32Pno(iV) 



27N{N + 1){N + 2f 
4 



-^(7V- 70)^1^ -^(iV-10)C25i 
16 



g(7V - 70)^2-51 + -(537V + 70)^3 - -j{N - 10)^2,1 

4(8A''^ - 663A'- - 115 lA' - 516) 2 8(2A'2 + TA" + 2)C2 

H ^^,^r . -.X/,. . '-'1 



27(Ar + l)(Ar + 2) 
4(124Ar3 + 663Ar2 ^ ^q^ON + 660),S2 



3(Ar + 2) 

SPiuiN) 



27{N + 1){N + 2) 



S1{N +1)^{N + 2) 



16Pn2(A^) 



+ ^(150 - N)St + ^{227N + 510)^4 



81N^{N + iy{N + 2y 
-^{N - 150)S,S, -liN- 30)52,1^1 + ^(30 - N)C2Sl 
+^(150 - N)S2S! + ^{IIQN + S70)Si + hllN + 30)C252 



Si 



2{20N^ - 1475A^2 _ 3278A^ - 1204) 
81(Ar + l)(Ar + 2) 



--(A^ - 30)^3,1 + 30)^2,1,1 + 

2(2Ar3 - 257Ar^ - 542Ar - 172) . , 28 
+ 9(Ar + l)(Ar + 2) + " 



Sf 
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with 



2(20A^=^ - 1475A^2 _ 3278Ar - 1204) 56(2A^2 ^ 7^ ^ 2) 

^ 27{N + T){N + 2) ' ' ^ 9(7VT2) 

4(574Ar3 + 1475A^2 + i^^qn + 1852) 8(2Ar3 _ 2577V2 - 542A^ - 172) 

"'-'3 H r./ Ar , -1 \ / AT , r.\ '-'2,1 



81(A^ + l)(A^ + 2) " 9(A^+ l)(A^ + 2) 

2Pn3(iV) ^2 ^ 4Pno(iV) 2Pn4(iV) 



2 



81(A^ + l)2(iV + 2)2 ^ 9Ar(Ar + i)(iV + 2)2^ 81N{N + iy{N + 2) 

^Pu.iN) ^5^^ 8^,5^^, (0.131) 



243(Ar+ l)3(Ar + 2)3 243N^{N + 1)^{N + 2y 



Puo{N) = 22A^5 + 134A^^ + 204Ar3 + 25A^2 - 88A^ - 12 , (D.132) 
Pni{N) = 140A^^ + 3958Ar^ + ISIOIA^^ + 22156Ar^ + 15566A^ + 6088 , (D.133) 
Pii2{N) = 458iV^ + 4306iV^ + 15574iV^ + 28027A^^ + 26233A^^ 

+ 12193A^2 + 2798A^2 + 636iV + 72 , (D.134) 
Pii3{N) = 290A^^ + 9918Ar^ + 39839 A^^ + 62088Ar2 + 46622Ar 

+18264 , (D.135) 
Piu{N) = 1094A^6 + 9414iV^ + 313917V^ + 522967V2 + 41630Ar2 

+11064iV - 864 , (D.136) 
Pii5{N) = 4276A^^ + 68264A^^ + 351618iV^ + 915849iV^ + 1394436Ar^ 

+ 1270848A^2 + 610766Ar + 89584 , (D.137) 
PiwiN) = 4426A^" + 52400iV^° + 256348A^^ + 673350iV^ + 1012896iV^ 

+836430iV6 + 281102iV^ - 71327iV^ - 82696iV3 - 22728Ar2 

-4464iV-432. (D.138) 



i(3),Qi 

Aq^ =0. (D.139) 



x(3),Rl 

Aq, =0. (D.140) 



i(3),Ria 

An, =0. (D.141) 



X(3),fl'ii, 

Aq, =0. (D.142) 
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x(3),5i„ 



{4 


16(iV + 2) 




9 



1 

+- 



1 

+- 

e 



--(TV + 2)^1 + 



8(327V4 + 73iV3 + lOiV^ - UN - 24) 
277V(7V+ l)(7V + 2) 



--(TV + 2)5? - ^(7V + 2)S, --{N + 2)C2 



4(32Ar4 + 73Ar3 ^ ^giV^ - 64iV - 24) 



27N{N + 1){N + 2) 
4Pu7{N) 



81Ar2(Ar+ l)2(Ar + 2)2 



-^{N + 2)S!-^{N + 2)Ss 



+^{N + 2)C3 - y + 2)^2^1 - ^(A^ + 2)C25i 



l(32Ar4 + 73iv3 + 10iV2 - 6AN - 24) ^2 
27A^(A^ + l)(A^ + 2) 



^1^ 



13(32A^^ + 73A^=^ + lOA^^ _ _ 34) 



27A^(A^ + l)(A^ + 2) 
(32A^^ + 73A^=^ + lON^ - 6AN - 24) 



9A^(A^ + l)(A^ + 2) 



C2 

2Pii8(A^) 



81N^{N + iy{N + 2y ' 2A3N%N + iy{N + 2y 



(D.143) 



with 



Pu7{N) = 352N^ + UGGN^ + 775N^ - ddON"^ + 4:36N^ + m6N^ 

+1200A^ + 288, (D.144) 

PusiN) = 2432Ar^° + 12004Ar'' + 24210A^^ + 36105A^^ + 56184A^^ 
+49188A^^ - 15320A^^ - 62080A^^ - 47040A^^ - 19584Ar 

-3456 . (D.145) 



T^UfCA \ 1 


16(iV-2) 


Ar2(Ar+ 1)2 


9 



1 



8 8(4Ar4 + 43iV3 + 26Ar2 - 64iV - 24) 

9^ ~ '^^ 27Ar(iV + l)(iV + 2) 

2(jv - 2)Sl + ^(7V - 2)S, + \{N - 2)C2 
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Si 



4Pii9(iV) 



A{AN^ + A3N^ + 26N^ - 6AN - 24) 
^ 27N{N + 1){N + 2) 

1 n n 14 13 

+-{N - 2)S! + —{N - 2)Ss - -{N - 2)C3 + -{N - 2)S2Si 



81A^2(7v + l)2(iV + 2)2 



+ -(N-2)C2Si + 



{AN^ + A3N^ + 26N'^ - 6AN - 24) 



27N{N + 1){N + 2) 
13(4Ar4 + 43iV3 + 2QN^ - QAN - 24) 

27N{N + 1){N + 2) ^ 
(4Ar4 + 43iV3 + 26N^ - QAN - 24) 



Si 



9N{N + 1){N + 2) 
2^119 (iV) . ^ 



C2 

2Pi2o(iV) 



81Ar2(iV + l)2(iV + 2)2 243iV3(iV + l)3(7V + 2)3 



(D.146) 



with 



Pii9(iV) = 152A^^ + 722A^^ + 461A^^- 1186A^^- 364Ar3 + i976Ar2 

+1200A^ + 288, (D.147) 

P^2o{N) = 1168A^^° + 6028A^'^ + 8238A^^ + 5523 A^^ + 27480 A^^ 

+51132A^5 - 1240A^^ - 62080A^^ - 47040A^^ - 19584A^ 

-3456 . (D.148) 



j(3),Tla 



1 64 1 



(A^ - 1)A^(A^ + 2) I 9 e 



6A{7N^ + 6N- 10) 
" 27(A^+ l)(Ar + 2) 

32(7iV2 + 6A^ - 10) 



32 



1 

+ - 

e 



8^2 104^ 8^ 



^1 



27(Ar+ l)(Ar + 2) 
64(25iV^ + 45Ar3 _ I4iv2 + 36iV + 124) 
81(Ar+ l)2(iV + 2)2 



4 e3 440^ 

+ 27"^^ + ~2rf^' 



^52 4 8(7iV^+6iV-10) 2 

+ 'Q^^^' + 3^^"^^ " 27(iV+l)(Ar + 2)'^^ 
104(7A^2 + 6A^ - 10) 8(7A^2 + 6A^ - 10) 

-'-'2 ^/^T- . -N/^^ . C2 



27(A^+ l)(A^ + 2) 9(A^ + l)(Ar + 2) 

32(25A^^ + 45A^3 _ ^4^2 ^ 35^ ^ ^24) 



81(A^ + l)2(Ar + 2)2 

61-f\2i(A0 



-^1 



243(Ar + l)3(Ar + 2)3 



(D.149) 



with 



P121 {N) = 79 iV^ + 336Ar5 ^ 725^4 ^ 1488^3 ^ i392Ar2 _ 744Ar 

-1480 . (D.150) 
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x(3),rii, 



TlufC, 



{N -l)N{N + IY{N + 2)\ e 
1 

+- 





64(iV + 3) 




9 



lOiV - 36) 



1 

+- 

£ 



32 32(lliV3 + 39iV^ 

+ 27(7V+l)(iV + 2) 

|(iV + 3)5,^ + :^(7V + 3)^2 + + 3)C2 



16(11A^3 + 39A^2 + lOA^ - 36) 
27(A^ + l)(A^ + 2) 

?3 



^1 + 



16Pl22(A^) 



81(A^ + l)2(Ar + 2)2 



4 44n 56 

+ ^(^^ + 3)^? + ^{N + 3)^3 - -{N + 3)C3 

,52 +4 3V5 4(lliV3 + 39iV^ + lOiV - 36) 
+-(7V + 3)^2^1 + -(iV + 3)C25i 27(iV+l)(iV + 2) 

52(llAr3 ^ 39JY2 ^ ;lq^ _ 4(lliV3 + 39A^2 + ^Q^y _ gg) 

27(A^ + l)(A^ + 2) ^ 9(Ar + l)(Ar + 2) 

81(Ar+l)2(iV + 2)2 ' 243(Ar+l)3(iV + 2)3 j ' V -^^ 



■C2 



with 



P^22{N) = 58A^^ + 225A^^ + 181A^3 ^ 324A^^ + 1372A^ + 1296 , (D.152) 
Pi23(iV) = 332Ar7 + 2418A^6 + 9531A^^ + 27393A^^ + 45426A^3 

+24108A^^ - 24536A^ - 25920 . (D.153) 



x(3),Tic 



1 



{N -l)N{N +lf{N + 2) I £3 
1 



32(A^2^3iV + 4) 
9 



+ - 



+ - 



^^(N^^^N^A\R 16(14iV^ + 51iV3 + 57iV^ - 187V - 56) 
-(TV + 37V + 4)^1 27(iV+l)(7V + 2) 

4 c;2 4 

-(iV2 + 3iV + A)Sl + — (iV2 + 3iV + 4)5'2 + -(iV^ + 3iV + 4)C2 

y y o 



8(14Ar4 ^ 5^^3 ^ 57^2 _ ;^g^ _ 
27(A^ + l)(A^ + 2) 



^1 + 



S1{N + 1Y{N + 2Y 



+ ^{N^ + 3iV + 4)^3 + ^(iv2 + 3iV + 4)^3 - ^(iV^ + 37V + 4)C3 
27 27 9 
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26 2 

+ (iV2 + 3iV + 4)S25i + -{N^ + 3N + 4)C2^i 
y o 

2(14Ar4 + 51N^ + 57iV2 - 18N - 56) 2 
27(iV+ l)(iV + 2) ^ 
26(14iV^ + 51A^3 + 57Ar2 _ _ gg^ 

27(A^+ l)(A^ + 2) ^ 
_2(14A^4 + 51A^3 + 57A^2_i8^_56) 4Pi24(iV) 

9(A^ + l)(Ar + 2) 81(A^ + l)2(Ar + 2)2 ^ 

4Pl25(iV) 

243(A^ + l)3(Ar + 2)3 

with 

lOOA^*^ + 438A^^ + 629A^^ + 393Ar^ + 1052A^^ + 2508A^ 

+1792 , (D.155) 
632A^^ + 4284A^^ + 14546A^^ + 36363A^^ + 68457A^^ 

+73650A^^ + 11500A^2 - 54168A^ - 37760 . (D.156) 



(D.154) 



Pl24{N) = 
Pl25(iV) = 



with 



(iV- 1)N^{N + 1){N + 



1 512 



+- 



+- 



128(3iV4 - 14iV3 - 21iV2 + 16iV + 4)" 
3(A^- 1)A^(A^+ l)(A^ + 2) 

256 9 64^^ + 64A26(A^) 
-256^2 - 64C2 + 3(^_ 1)2^2(^ + 1)2(^ + 2)2 

448 64(3A^4 - 14A^3 _ 2iN'^ + iqn + 4) 

— H 7TT TTTTTTT TTTTT So 



QC4. c I " r -L • 

' 3^' {N -l)N{N + l){N + 2) 

16{3N^ - UN^ - 21Ar2 + 16Ar + 4) 

^ {N -1)N{N + 1){N + 2) 

32Pi27(iV) 

3(A^ - l)m^N + lf{N + 2)3 



(D.157) 



Pi26{N) = A^s + 34Ar7 _ 457V6 _ 239A^^ + 38A^^ + 269A^^ - 86A^2 

-36Ar - 8 , (D.158) 
Pi27{N) = 4iV^2 + 34JY11 + 286iVi° - 20N^ - 18A5N^ - 826iV^ 

+3276A^^ + 632N^ - 2657N^ + 396N^ + 192N^ + SON 

+16 . (D.159) 



95 



X(3),B2 



N^{N + 1)2 I £3 



64(A^- l)(iV + 2) 



1 

+- 

£ 



32Pi28(iV) 



3{N -1)N{N + 1){N + 2) 
-32{N -1){N + 2)S2 - 8{N - 1){N + 2)C2 

16Pl29(A^) 



3{N - iyN^{N + iy{N + 2y 

56 



-48(A^- l)(A^ + 2)53 

16Pl28(A^) 



+ — (N -1)(N + 2)C3 + 7 ^^ov- y 



+ 



4Pl28(iV) 



C2 



{N -1)N{N + 1){N + 2) 

8Pi3o(iV) 

"3(Ar-l)3iV3(Ar+l)3(Ar + 2)2 



(D.160) 



with 



Pl28{N) 
Pl29{N) 

Pl30{N) 



2N^ + AN^ - 13N^ + UN^ + AIN^ - UN - 12 , (D.161) 

jyw + - ISN^ + lOON^ + A6N^ - AMN^ - 212A^^ 

+A73N^ + nON^ - WON - 56 , (D.162) 

4A^^^ + 26A^^=^ + 67N^^ + 13AN^^ - 719A^^° - 1462Ar9 

+3010A^^ + 5074iV^ - 4187 A^^ - 5412Ar5 ^ 2609A^^ + 

2240A^=^ + 296Ar2 - 576A^ - 240 . (D.163) 



j(3),C2 



1 





2048 




3 



+- 



256P 



131 



+ - 



3(7V- l)Ar(Ar+ l)(Ar + 2) 

128Pi32(iV) 



+15365. 



1024^2 + 256C2 + 
1792C3 



3(Ar-l)2Ar2(Ar+l)2(Ar + 2)2 
128Pi3i(iV) 



3 {N -1)N{N + 1){N + 2) 
32Pi3i(iV) 



^2 



C2 



{N -1)N{N + 1){N + 2) 

64Pi33(iV) 

'3{N-lfN^{N+lY{N + 2f 



(D.164) 
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with 



Pl3l{N) 
Pl32{N) 



+ 9iV^ - 25N^ - 53iV2 - 4 , 



-UN' 



?4iV- 



4iV^^ + 10A^^2 + 38iV" + 568iV^° + 201iV^ - 3687iV^ 
-3872iV^ + 3330A^^ + 2677iV^ - 15737V^ + 6487V^ 
_7447V2 _ 1767V - 16 . 



(D.165) 
(D.166) 

(D.167) 



A, 



with 



Qg 



1 



Ca\_ 

2 J N^(N + lf ] £3 

128Pi34(A^) 
3N{N + 1){N + 2y{N -1 

32Pi35(iV) 



1 256(A^ - 1) 



+32 (AT - 1)^2 

224 



:(7V-i)C3 



3(iV- l)2iV2(iV+ l)2(iV 

64Pi34(iV)52 



128(Ar - 1)^2 

+ 192(Ar - 1)^3 



A^(A^ + l)(A^ + 2)2(iV 

16Pl34(A^) 



C2 



N{N + 1){N + 2y{N -1) 

16Pi36(iV) 



(D.168) 



(D.169) 



Pi34{N) = 2N^ + 6N^ - 9A^^ + 16 A^^ + 53N^ + SN-4:, 
Pi35{N) = A^ii + 4A^^° + 17A^'^ + 58A^^ - 237A^^ - 278A^^ + 1057A^^ 

+1156A^^ - 358A^^ - 420A^2 - 152A^ + 16 , (D.170) 
Pi36{N) = 4Ari5 + 22A^^^ + 79A^^^ + 95N^^ - 575A^^^ + 335A^^° 
+3743A^'' - 3193A^^ - 13937 A^^ - 1437A^^ + 10542Ar5 
+3490A^^ - 832A^^ - 2752 A^^ - SOOAT + 32 . (D.171) 



j(3),B2 



T^UfC, 



(Ar + l)2(Ar + 2) £3 



+- 



32(2A^2 _ 13JY _ 
9(Ar-l)Ar 

_32 16Pi37(iV) ' 

9 ^ 27(A^-l)2Ar2(Ar + l)(Ar + 2) 



+- 



8(13A^2 _ 7gjY _ 193) 
9(Ar- l)Ar 



^2- 



4{2N^ - 13A^-33) _8 2 
3(Ar- l)Ar ^'~9 ' 



16(26Ar4 + 103Ar3 + 253A^2 + 338Ar + 144) 
^ 27(Ar- l)Ar(Ar + l)(Ar + 2) ^ 
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with 



8Pl38(iV) 



88{5N^ - 32N - 81) 
27(A^- 1)A^ ' 



81(A^ - 1)^N^{N + l)2(iV + 2)2 
, 28(27V2_i37V-33) ^ 52 4 4 



S2 + 



27{N -iyN^{N + 1){N + 2) " 9{N - iym{N + 1){N + 2) 



:C2 



4(26A^^ + 103A^3 + 253N^ + 338N + 144) 
27{N - 1)N{N + 1){N + 2) 
8Pi4i(iV) 



^1^ 



Si 



81{N - 1)N{N + l)2(Ar + 2)2 

4Pl42(iV) 

' 243(Ar - l)4Ar4(7V + 1)3(7V + 2)^ 



(D.172) 



Pis^{N) = 52N^ + 24:N^ + 28N^ + 765N^ + 532N'^ - 555N - 54 , (D.173) 
Pi38{N) = 572A^^° + 1696 A^'' + 2460A^^ + 207A^^ - 18024A^^ - 24537A^^ 

+16705A^^ + 20618A^=^ - 12729A^2 - 900A^ - 324 , (D.174) 
Pi39(Ar) = 338A^^ + 221Ar^ + 318A^^ + 4675A^^ + 2998A^^ - 3474A^ 

-324 , (D.175) 
Pi4o{N) = 52A^^ + 24Ar5 + 28A^^ + 765A^^ + 532A^2 - 555A^ - 54 , (D.176) 
Pui (N) = 286A^^ + 1856A^^ + 5575A^^ + 8375A^^ + 4624Ar2 _ 1276A^ 

-1296 , (D.177) 
Pu2{N) = 4276Ar^^ + 16804Ari3 + i0992iVi2 _ gujgNii _ i81276Ario 

+209316A^^ + 783074iV^ - 83677iV^ - 1043730A^^ + 235534Ar5 



+660610A^* - 261267Ar'^ - 4554Ar2 _ 6372A^ - 1944 . 



(D.178) 




+- 



+ ■ 



16 



16(5A^3 + 15iV2 + 5N- 17) 
9iV(iV + 1)2 

SPu3{N) " 

9 27{N -1)N^{N +1)^{N + 2) 

A{31N^ + 92N^ + 31N - 102) 
^ 9N{N + iy ' 



^1- 



2(5A^=^ + 15A^2 + - 17) 4 2 
^ 37V(7VTlj2 ^ 9^ 



27A^(A^ + l)(A^ + 2) 



Si + 



81{N - iym{N + 1)^{N + 2)2 
4(136A^3 + 407A^2 ^ 135^ _ 459) 14(5A^=^ + 15A^2 + 57V - 17) 

H ^/ ^ s„ : r-^ 



27A^(A^ + 1)2 
26 2 2 n 



9A^(A^ + 1)2 
2{5N^ + 3QN^ + 3AN + 60) 



27A^(A^ + l)(Ar + 2) 



-St 
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2Pi45(iV) ^ PMN) , 



27(N -1)N^{N + 1)^(N + 2) " 9(N - 1)N'^(N + 1)^(N + 2) 



81N^{N + iy{N + 2)^ 

2Pl48(iV) 

'2A3{N - l)^m{N + 1)^{N + 2)3 



^1 

(D.179) 



with 



Pu3{N) = A6N^ + 251N^ + 255 A^^ - 250A^^ + 234 A^^ + 737 

-24:7N - 234 , (D.180) 

Pi44{N) = 146Ar^^ + 1928iV^° + 4491iV'' + 869iV^-3135iV^-10143iV6 

-26309A^^ - 3731iV^ + 25311iV'^ + 4129iV2 - 4824A^ - 2988 (D.181) 

Pu5{N) = 281 A^^ + 1541iV^ + 1589 AT^ - 1441iV^ + MOOAT^ + 4328Ar2 

-1542A^ - 1404 , (D.182) 

PueiN) = AQN'^ + 251N^ + 255N^ -250N* + 23AN^ + 737 -2A7N 

-234 , (D.183) 

Pl^^{N) = 35N^ - 243N^ - 1303Ar^ - 2082Ar3 - 136A^^ + 2136Ar 

+720 , (D.184) 

Pi48(Ar) = 454iV^^ + 15564A^^^ + 57172A^^^ + 34999A^^^ - 143229A^^^ 
-317074A^^° - 15876A^^ + 8501 16A^^ + 667430Ar^ 
-871408A^^ - 712308A^^ + 468109A^^ + 287657A^^ 
+51030Ar2 - 79092Ar - 36936 . (D. 185) 



A. 



Qg 



N{N + 1){N + 2) £3 



32(Ar-4) 32 (Ar2 + 2Ar + 4) 
-Si + ■ 



9 



9(A^ + 2) 



1 8(7V-4) 8(7V-4) 
[ 3 + 3 

16 (26A^^ - 16N^ - S7N^ - S3N - 128) 
27(A^- l)(A^ + l)(A^ + 2) 

16Pi5o(A^) 
27(Ar - 1)A^(A^ + 1)(A^ + 2)2 



Si 



1 

+ - 

e 



28(Ar-4) ^3 
27 



9 



16(7A^ + 12) 8(17A^ + 28) 

^^2,1 7ZZ '-'3 



27 



4 {28N^ - 12N^ - 35N^ + 19N - 96) ^2 4 (A^^ + 2Ar + 4) 



9(Ar- l)(Ar+ l)(Ar + 2) 
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st + 



3{N + 2) 



<2 



with 



4 (12Ar4 + i72Ar3 + 215N^ - 159N - 144) 
^ 9(A^-l)(A^ + l)(A^ + 2) ' 

8Pi5l(7V) „ 8Pi52(A^) 

T'-'l — 



81(A^ - 1)(A^ + iy{N + 2)2 ' 81(A^ - 1)N'^{N + iy{N + 2f 
20(57V - 4) 8(17iV + 36) 29(7V - 4) 

+ ^ -^3-^1 ^ ^2,1^1 + ^ -^2 

(181iV + 204)^ 8(5iV + 36) 



9 

2(200A^^ - 76A^=^ - 141A^2 ^ 337^ _ 598) 

'-'1 



81(A^- 1)(A^+ l)(A^ + 2) 

2(26A^^ - 16Ar3 - 87iV2 _ §3^ _ ^28) 8(17A^ + 36) 

-C2>-'i H p; '-'2,1,1 



9(iV- l)(iV+ l)(iV + 2) ' 9 

2(296iV4 - 940 _ 2l33iV2 - 359Ar - 608) 28 {N^ + 2iV + 4) 

-S2S1 , , ^3 



27(iV- l)(Ar+ l)(iV + 2) 9(Ar + 2) 

4(676iV^ + 2332A^3 ^ i947jv2 - 2035Ar - 2632) 

^ 81(A^- l)(iV+ l)(iV + 2) ^ 

8(2A^3 - 150A^2 _ 491^ _ 384) 

+ 27(7V + l)(7V + 2) -'^'^ + - 

, 2PMN) 2 , 2Pi5o(A^) 

27(Ar - l)(Ar + l)2(Ar + 2)2 1 ^ 9(A^ - 1)A^(A^ + 1)(A^ + 2)2^^ 

2Pi53(iV) ^ 4Pi54(iV) ^ 

27{N -l)N{N+lf{N + 2f ^ 243{N-l){N + iy{N + 2f ^ 

243(Ar-l)Ar3(Ar+l)3(Ar + 2)4 ' ^ ^ ^ 



Pi49(A^) = 242 A^^ + 246 A^^ - 495 A^^ - 958A^^ - 2384A^2 - 1979 

+1296 , (D.187) 

Pi5o{N) = 13Ar6 + 82A^^ + 16Ar^ - 353Ar=^ - 446A^^ - 152Ar - 24 , (D.188) 
Pi5i{N) = 268A^^ + 412A^^-1147A^^-4532A^^- 6520A^^-2513Ar 

+1936 , (D.189) 

Pi52{N) = 53A^'^ + 341A^^ - 693A^^ - 6717A^^ - 13305A^^ - 8433Ar^ (D.190) 

+2422 A^^ + 4324 A^2 + 1128Ar + 144 , 
Pi53{N) = 342 A^^ + 3314Ar6 ^ 8767 A^^ + 8398 A^^ + 2732 A^=^ + 879 A^^ 

+48Ar + 288, (D.191) 
Pi54{N) = 1940Ar« + 8156Ar^ + 9782A^'^ - 20076Ar5 - 94500Ar^ 

-153048A^^ - 142667A^^ - 123635Ar - 76928 , (D.192) 
PibbiN) = 2833A^^2 + 32644A^^i + 146294A^i° + 316196A^^ 

+266628A^^ - 211206A^^ - 715870A^^ - 709393A^^ 

-362414Ar^ - 119936Ar3 - 36528Ar2 - 8064Ar - 864 , (D.193) 
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with 



X(3),i?2 



TlrifCA 



(Ar + l)2(iV + 2) £3 



1 



9(iV- l)iV 

32 16Pi56(iV) ' 

9 ^ 27(iV-l)2iV2(iV + l)(iV + 2) 



1 

+ - 



8 2 8(13iV2 - 79iV - 198) 



4(2iV2 - 13iV - 33) 



9'' %N-l)N 3{N-1)N 

16(26iV^ + 103N^ + 253N^ + 338iV + 144) 



C2 



27{N - 1)N{N + 1){N + 2) 

8Pl57(iV) 

■81(Ar - 1YN%N + iy{N + 2y 

28(2A^2 _ _ 



-^1 

88(5iV2 - 32A^ - 81) 
27{N-1)N ' 



52 



9(A^- l)iV 



Cs ;:-'S'2S'i — -C2'S'i — —S 



9 



27 



A{26N^ + 103N^ + 253N^ + 338N + 144) 
27{N - 1)N{N + 1){N + 2) 

4Pl58(iV) 

27(A^- l)2Ar2(Ar + i)(Ar + 2) ^ 

2Pl59(A^) 



C2 



^1 



9(N -iyN^(N + 1)(N + 2) 

8Pi6o(A^) 

'81(A^ - 1)N{N + 1)^(N + 2)2 

4Pi6i(A^) 

'243(A^ - l)4Ar4(7V + 1)^{N + 2)= 



(D.194) 



(D.195) 



Pi5g(iV) = 52Ar^ + 6A^^ - 26A^^ + 783A^^ + 658A^^ - 555A^ - 126 , 
Pi57(iV) = 572A^^° + 1480A^'^ + 1920A^^ + 1233A^^ - 15432A^6 

-26697A^^ + 12385A^^ + 23048A^^ - 10029A^^ - 1980A^ 
-756 , (D.196) 
Pi58{N) = 338A^^ + 113A^^ - 6A^^ + 4783A^=^ + 3754A^2 - 3474A^ 

-756 , (D.197) 
Pi59{N) = 52Ar^ + 6A^^ - 26A^^ + 783A^^ + 658A^2 - 555A^ - 126 , (D.198) 
Pieo{N) = 286 A^^ + 1856 A^^ + 5575 A^^ + 8375 A^^ + 4624A^2 - 1276 

-1296 , (D.199) 
Pi6i{N) = 4276Ar^^ + 15508Ar^=^ + 8400Ar^2 _ 729i6Arii - 173500A^^° 
+169140Ar'^ + 778700A^^ + 23729A^^ - 1060254Ar6 
+98482Ar5 + 694468Ar^ - 190635Ar^ - 20106Ar2 - 14148Ar 
-4536 . (D.200) 
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D.2 A 



PS 
Qq 



The individual contributions to Aq,^ are: 



j(3),PS,a 
An„ 



T^UfCp j 1 



N^{N + iy 



^{2 + N){-l + N) 



1 



1 

+ - 

e 



32 32 
■¥P + ^)(-l + '^)^' + 27]V(l + iV)(2 + iV) 



104 



— (2 + N){-1 + N)S2 - g (2 + N){-1 + N)Sl 



8 16 
-3^2 + ^)(-l + ^)C. + 27iV(l + 7V)(2 + 7V) 



-Si 



16 



Pl63{N) 



81 Ar2(i + Ar)2(2 + Ar)2 
52 



440 



{2 + N){-l + N)Ss 



-- (2 + N){-1 + N)S2S, - 27 (2 + + N)S! 

56 4 
+- (2 + N)(-l + N)C3 - 3 (2 + N){-1 + N)SiC2 



+ 
+ 



52 Pi62(iV) 



27 A^(l + A^)(2 + A^) 

4 /\o2(A0 



^2 + 



Pl62(iV) 



C2- 



27A^(1 + A^)(2 + A^) 



9 N{l + N){2 + Ny 81 N^{l + Ny{2 + N)^ 

8 Pi(ii(AO 



-Si 



243 N^{1 + N)^{2 + N)^ 



(D.201) 



with 



Pw2{N) 
PwsiN) 

Pl64{N) 



8N^ + 29 -9 -8 + QAN + 2A , (D.202) 

25 A^^ + 151 A^^ + 116 + 352 + 1052 A^^ (D.203) 

-428 A^^ - 2264 A^^ _ 1200 A^ - 288 , 

158 A^^^ + 1505 N^^ + 5261 A^^ + 12912 A^^ + 13860 A^^ 

-16140 A^^ - 33344 A^^ + 27880 + 77344 A^^ 

+50496 + 19584 TV + 3456 . (D.204) 



~ (A^ - l)N{N + l)(Ar + 2) j~i3-9- 

256 1024 (2Ar2-5) 
~9~ ^ ^ ~W {2 + N){l + N) 
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1 

+- 

£ 



_^c2_ 64 512 (27V2-5) 
9 ' 9 ^ 3 27 (2 + A^)(l + A^) ' 



128 (25 - 42 iV3 - 107 iV^ + 348 TV + 532) 
"sT (2 + Ar)2(l + Ar)2 

3520 



416 32 3 32 448 



1664 (2A^2_5) ;^28 
H — :rz— t::^- — tttt:: — ■ — ttt '-'2 + 



(2A^2_5) 



27 (2 + A^)(l + Ar) 
128 (2A^2_5) 



27 (2 + A^)i(l + AT) 



1 -^1 



C2 



9 {2 + N){l + N) 
64 (25 - 42 iV^ - 107 iV^ + 348 + 532) 
81 (2 + iV)2(l + Ar)2 1 

128 Pi65(iV) 



243 (2 + 7V)3(l + 7V): 



(D.205) 



with 



P\Qb{N) = 79 iV^ + 411 iV^ + 2085 Ar^ + 5289 Ar3 + 3252 

-5484 iV - 6064 . (D.206) 



-PS 



-PS 



The contributions toA^^ Q read: 



^(3),PS,a 



(iV- l)iV(iV + l)(iV + 2) 





512 


1 




9 





512 512(2iV+ l)(4iV + 7) 



9 



27(iV + l)(iV + 2) 



+ - 



256 



- -C2 - —Si + ^^^^^^ + ^^^^^ + 



64 



256 
~9" 



27(Ar + l)(Ar + 2) 



^1 



128(181A^^ + 894A^3 ^ 15977V2 + 1248A^ + 400) 



256 ^3 



27 



81(A^+ l)2(Ar + 2)2 

512 448 (2Ar + l)(4Ar + 7) 
+ + (]v + l)(iV + 2) 



H — ^'S'2S'i + "^C2'S'i 



128(181Ar4 + 894iV3 + 1597iV2 + 1248A^ + 400) 
81(iV+ l)2(7V + 2)2 

128Pi66(iV) 



256 



Si 



256 



So 



64 



243(Ar+ l)3(Ar + 2)2 



(D.207) 



with 



Pim{N) = 1037iV^ + 8247iV^ + 26940A^^ + 46191Ar3 + 43809Ar2 

+21648Ar + 4192 . (D.208) 
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j{3),PS,b 
A 

^qq,Q 



T^UfCp 1 64 ,,,, ^ 



1 



^(2 + 7V)(-l + 7V)5, 



32 16 + 26 _ 25 - 11 AT + 6 
'27 



1 

+ - 



iV(l + N) 



(2 + N)i-l + N){-fs,-fsl-lc2 



32 {16N'^ + 26N^ -25N'^ -UN + 6) 16 Pi67(A^) 
^27 ]V(TTaO ^ ~ 81 N^{l + Ny 

+{N + 2) (TV - 1) I +g ^3 + ^ 5,5, + + l 



56 \ 16 (16A^^ + 26iV3-25iV2-lliV + 6) 
^"9" ~ 27 iV(lTiV) ' 

16 (16 + 26 iV3 - 25 iV2 - 11 TV + 6) , 
~27 Ar(l + iV) ^ 

4 (16iV4 + 26iV3-25iV2-117V + 6) 16 Pi67(iV) 
~9 N{l + N) ^'^81 Ar2(i + iV)2 ' 



8 Pl68(A^) 



243 iV3(l + iV)3 



(D.209) 



with 



Pl67(iV) 



181 + 447 iV^ - 32 AT^ - 297 A^^ - 92 iV^ + 15 iV - 18 , (D.210) 
2074 A^* + 7210 A^^ + 4927 A^^ - 3503 A^^ - 5309 A^^ - 929 A^^ 
+231 AT^ + 9Ar + 54 . (D.211) 



D.4 A 



qq,Q 



;NS 



The contributions toyl^^ ^ are given by : 



j(3),NS,a 



N{N + 1) 1 £3 



^(Ar + 2)(Ar-i) 



1 

+- 

e 



6A {N^ + 2N^ -lON^ -5N + 3) 
81 iV(lTAO 



9 ^ 81 (1 + Ar)2Ar2 
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with 



+- (N + 2)(-l + N)Cs - 2^ ^ ^^^j^, ^C. 

_26__^7o(A^l 
729 (1 + NfN^ J ' 



(D.212) 



Pl69{N) 

Pi7o{N) 



+14 + 42 A^^ - A^^ - 50 A^^ - 19 A^^ + 2 - 3 , (D.213) 

+308 + 1232 A^^ + 395 A^^ - 2353 A^^ - 2413 A^^ 

-391 A^^ + 153 AT^ + 9 AT + 27 . (D.214) 



A 



;(3),NS,b 
qq,Q 



TpnfCF< 



128 128 



+ 



128 64 320 

"sT +27 '^^""sT'^^ 



+- 



27 9 ^ 81 9 ^ 27 



Si 



m 

81 



+ ^S +^^2^9 +^r 9 ^^^^ 9 

16 160 13888 9856 
+ 27 ^^ + ^'^^"^29"'^^ + 729" 



(D.215) 



x(3),NS,c I 1 16 120 2 ^ 337 



qq,Q 



= TinfCfl — C2 - 

^ ' £2 9 ^27 3 ^ 



(D.216) 



D 5 4NS,TR 



;NS,TR 



The contributions io A^^ q read: 



-(3),TR,a 



1 64 

~2 fiY 



1 r32(14A^2 ^ 14^ _ 3) 
~e 81iV(ivTT) 

16(308Ar4 + 616A^3 + 263A^2^9Ar + 27) 8 

729iV2(iV + 1)2 27^^ 



56 



(D.217) 
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x(3),TR,b 
A 

^qq,Q 



128 128 



1 



64 



320 128 



1 

+- 

£ 



32 16 160 16 320 896 

27 9 ^ 81 9 ^ 27 81 



S 8Q 40 160 

+ 27"^' + + 9^"^' " Sl"^' " 27^"^' " + ^"^^ 

16 13888 98561 
+ 27^^"^29~'^^ + T29"/ ■ 



(D.218) 



x(3),TR,c 



99>< 



1 16 1 20 2 337 



(D.219) 
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E Variable Transformations 



In many cases the structure of the emerging Feynman parameter integrals can be simphfied 
by applying transformations to the integration variables, which were given in Ref. [217] . 

• To define the product x' — xy as the new integration variable, one maps: 

/ / x{l-y) 

X -.^ xy , y := , 

l-xy 

x' 

X = x' + y' - x'y' , y = — — — , 

y' + x' — x'y' 

d{x, y) 1 — x' 



d{x', y') x' + y' — x'y' 

One obtains 

/' /' d-dy fix, y){xyr = /' /' dx'dy' ^^^^ 
Jo Jo Jo Jo X + y — X y 



(E.l) 



f(y' + x'- x'y', ■ f ) . (E.2) 
^ x' + y' — x'y' ' 



Terms of the form [x — y) can be combined by 



X > y : X < y : 

x :— X — y , x :— y — X , 

y / 1 -y 



1 — X + y 1 + X — y 

x^x' + y' -x'y' , X ^ (1 - x'){l - y') , 

y = (1 - x')y' , y = 1 - (1 - x')y' , 

d{x',y') ~ d{x',y') ~ ^ ' 



Thus, one obtains 

/ / dxdy f{x,y){x-yY ^ dx'dy' 
Jo Jo Jo Jo 



x (1 — x') [/(y' + — x'y , (1 — x')y'^ 
+(-l)^/((l - y')(l - 1 - (1 - Ay') ] ■ (E.4) 

If one applies Eq. (E.3) to factors [x — y) in the denominator, special care is needed 
to avoid possible divergences. The above transformations allow to simphfy the Feynman 
parameter integrals analytically. In this way higher transcendental functions like gener- 
alized hyper geometric functions, which obey single-sum representations are obtained, cf. 
Eq. (F.ll). 
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F Special Functions 



In the following we summarize for convenience some relations for special functions which 
occur repeatedly in quantum field theory calculations and are mutually used within the 
present work. 

F.l The Euler Integrals 

The F-function, cf. [218,219], is analytic in the whole complex plane except at the non- 
positive integers, where it possesses single poles. Euler's infinite product defines 



1 



F(^) 



; exp(7£z) Y[ 



i=l 



1 + - ) exp ( -- 



(F.l) 



The residues of the F-function at its poles are given by 



Res[F(^)],=_^ = 



m 



N eNUO . 



In case of Re{z) > 0, the F-function can be expressed by Euler's integral 

poo 

r{z) = / dtexp{-t) t^'^ , 
Jo 

from which one infers the well known functional equation of the F-function 

r(z + 1) = zr(z) , 



(F.2) 



(F.3) 



(F.4) 



which may be used for its analytic continuation. Around z — 1, the following series 
expansion is obtained 



F(l - £) = exp(£7E) expi XI Ci^ > , kl < 1 ■ 



i=2 



(F.5) 



Here and in (F.l), 7^ denotes the Euler-Mascheroni constant, cf. (C.8), and Riemann's 
^-function for integer arguments k, cf. (C.IO). A shorthand notation for rational functions 
of F-functions is 



ai, .... a-i 
bi,...,bj 



F(ai)...F(a^) 

m)...r{bj) 



(F.6) 



Functions closely related to the F-function are the Euler Beta-function B{A, C) function, 

the '4){x) — , and the /3(a;)-function. 

The Beta-function can be defined by Eq. (F.6) 



A + C 



B{A,C) = F 

If Re(74), Re(C) > 0, the following integral representation is valid 

B{A, C) = f dx x^-\l - xf-^ . 
Jo 



(F.7) 



(F.8) 
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For arbitrary values of A and C, (F.8) can be continued analytically outside of the re- 
spective singularities using Eqs. (F.l, F.7). Its expansion around singularities can be 
performed via Eqs. (F.2, F.5). The ip-iunction and /3{x) are defined as logarithmic 
derivatives of the F-function via 



1 d 



V{x) dx 
1 
2 



V{x) , 



x + 1 



(F.9) 
(F.IO) 



F.2 The Generalized Hypergeometric Functions 

The generahzed hypergeometric function pFg is defined by, cf. [136,137,220], 



fll, Up _ 



E 



i=0 



{ai)i...{ap)i 



(F.ll) 



Here (c)„ denotes Pochhammer's symbol 

(c)n 

for which the following relation holds 

(7V + 1)_, = 



F(c + n) 



(F.12) 



N e n . 



(F.13) 



In (F.ll), there are P numerator parameters ai...ap, Q denominator parameters bi-.-bg, 
and one variable z, all of which may be real or complex. Additionally, the denominator 
parameters must not be negative integers, since in that case (F.ll) is not defined. The 
generalized hypergeometric series pFq are evaluated at a certain value of z, which is 
always z — 1 ior the final expressions in this thesis. 

Gaufi was the first to study this kind of functions, introducing the function 2-^1, and 
proving the theorem, cf. [136], 



2Fi[a,6;c; 1] = F 



c, c • 



a — b 



Re(c - a - 6) > , 



(F.14) 



which is called GauB' theorem. An integral representation for the hypergeometric GauB- 
function is given by 



2F1 



a,b + 1 
c + 6 + 2' 



T 



c + b + 2 
c+ 1,6+ 1 



[ dx x\l - xY{l - zxy , (F.15) 
Jo 



cf. [136], provided that the conditions 

1^1 < 1 , Re(c+ 1), Re{b+1) > , 



(F.16) 
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hold. Applying Eq. (F.15) recursively, one obtains the following integral representation 
for the general hypergeometric function p+iFp: 



ao, fli, . . . , ap_ 
bi,... ,bp 



hi,... ,hp 
Qi, . . . , Op, h\ — (ii, • • • , hp — dp 



X 



/ dxi... [ dxp - xif'-'''-^ . . . x^/'-\l - xp)^^-"^-^(l -zxi... xp)- 

Jo Jo 



-ao 



(F.17) 

under similar conditions as in Eq. (F.16). As generalized hypergeometric functions ap- 
peared frequently during this computation, taking advantage of their properties is of 
importance. 

If one considers the fraction 



(F.18) 



with integers m > n > 1 one may transform the arguments in the following way by using 
the definition of Pochhammer's symbol : 



Mi 

in). 



r{n) r{m + i-l) m + i-l 
r(m) r(n + i-l) n + i-1 
r{n) r(m + i-l)/ m-n 
r(m) T{n + i~ 1) V ^ n-l + i 
r(n)T(m-l) . . r 1 



rim) 



[m 



(n-l)T(n-l) 



+ 



(m - n)r(A^) 



m — 1 {m — !)■ ni — n (ni — !)■ 



(F.19) 



n — 1 {m — !)■ m — 1 (n)^ 
This relation can be applied to an arbitrary function pFq and proves to be especially useful 
if one considers the generalized hypergeometric function of the form 3F2 
with m, n integers, n > m > 0. In this case repeated application of 



ai ,02 ,m . -j^ 
bi,n ' 



ai, 02, m 
bi,n 



n-1 
m — 1 

m — n 



+ 



ai,a2,m - 1_ ^ 
bi,n — 1 ' 

fli, 02, m — 1 



m — 1 

leads to a linear combination of terms of the form 



3^2 



hi,n 



(F.20) 



3^2 



hi,n 



(F.21) 



and terms that do not contain any sum any longer. These relations belong to the class 
of contiguous relations [136,220-223]. If ai, 02,^1 are non-integers or 01,02,61 > n — 1 
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(F.21) can be simplified by considering the following relations: 



Tin) 



Tin) 



{A) 



T{n + i) (n + i - 1)! ' 
T{A + i) 
T{A) 

r(A + i-(n-l) + (n-l)) 

T{A) 
T{A-{n-l)) 



(F.22) 



r(^) 

After applying (F.22) and (F.23) one obtains 



{A + l-n)^ 



i+n—l 



(F.23) 



3-^2 



iai)iia2)i 1 



r(ai + l-n)r(a2 + l-n)r(6i) 



I i 



i=0 



r(ai)r(a2)r(6i + l - n) 
(fli + 1 - n)i(a2 + 1 - n) 



(6i + 1 - n). 



[n 



z-1) 



r(ai + l-n)r(a2 + l-n)r(6i) 
r(ai)r(a2)r(6i + l-n) 

ai + 1 — n, 02 + 1 — n ^ 
bi + l-n 

{ai + l-n)i{a2 + l-n)i 1 



X 



2-^1 
n-2 



E 

i=0 



(6i + 1 - n)i 



(i)! 



(F.24) 



The infinite sum contained inside the hypergeometric function in (F.24) can now be 
evaluated by applying GauB' theorem, (F.14), such that no infinite sum remains. The 
relations (F.20) and (F.24) have been implemented into FORM-algorithms and in many 
cases made it possible to perform the infinite sums before expanding in the dimensional 
regularization parameter e. A similar useful relation is the following, which has been 
given in Ref. [224] and holds under the condition, that Re(e — a — b — 1 — k)>0 : 



3-f^2 



a,b,d + k ^ 


= r 


e 








d, e ' 




e — a,e — b 





r(e — a — b — 



AA(aM% 
ij {d)i 



(F.25) 



The following relations are restricted to special cases, but prove often to be useful. One 
defines the parametric excess of the series bys :— d + e — a — b — c. Saalschiitz's theorem, 
cf. [136], states that 



3-t^2 



a, b, c _^ 


= T 




d, e 





d, 1 + a — e, 1 + 6 — e, 1 + c — e 
1 — e,d — a,d — b,d — c 



(F.26) 



provided that s = 1, i.e. the series is Saalschiitzian, and one of the numerator parameters 
is equal to a negative integer. Another theorem is a generalization of Dixon's theorem. 
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[136], 



3-^2 



a,b,c ^ 


= r 




d, e ' 





d, e, s 
a,b + s,c + s 



3-1^2 



a,e — a,s ^ 



s + 6, s + c ' 

F.3 Harmonic Sums and Nielsen— Integrals 



(F.27) 



Expanding the F-function in £, its logarithmic derivatives, the ■^^'^^ -functions, emerge. In 
many apphcations of perturbative QCD and QED, harmonic sums occur, cf. [124,125], 
which can be considered as generahzation of the ■0-function and the ^-function. These 
are defined by 



N m 



Sai,... ,am{-^) 



;sign(ai))"i (sign(a2))"' (sign(a„))"'^ 



ni=l n2=l 



n™=l 



n 



ail 



a2\ 



N e N, V/a; eZ\o, 

1 . 



We adopt the convention 



-"ai,... ,am 



■'ai, 



(F.28) 
(F.29) 



(F.30) 



i.e. harmonic sums are taken at argument N, if no argument is indicated. Related 
quantities are the Z-sums defined by 



■'mi,... ,mk 



iV>n>i2...>iA:>0 



n ti[sign(m;)]^ 
Ami 



(F.31) 



The depth d and the weight w oi a harmonic sum are defined by 

d :— m , 

m 



W 



(F.32) 
(F.33) 



i=l 



Harmonic sums of depth d — 1 are referred to as single harmonic sums. The complete 
set of algebraic relations connecting harmonic sums to other harmonic sums of the same 
or lower weight is known [131]. Thus the number of independent harmonic sums can be 
reduced significantly, e.g., up to w = 4 the 80 possible harmonic sums can be expressed 
algebraically in terms of 31 basic harmonic sums only [132]. One introduces a product 
for the harmonic sums, the shuffle product LU , cf. [131]. For the product of a single and 
a general finite harmonic sum it is given by 

,b2,... ,bm,ai 

(F.34) 

For sums Sai^...,an{^) 'S'5j^...^5^(A^) of arbitrary depth, the shuffle product is then the 
sum of all harmonic sums of depth m + n in the index set of which Oj occurs left of aj for 
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i < j, likewise for bk and bi ior k < I. Note that the shuffle product is symmetric. 
One shows that the following relation holds, cf. [131], 



-SaiAbl,b2,-,bmi^) - ■■■ - Sbi^b2,...,aiAbm{^) > (F.35) 

where the A symbol is defined as 

aA& = sign(a)sign(6) (|a| + . (F.36) 

Due to the additional terms containing wedges (A) between indices, harmonic sums form 
a quasi-shuffle algebra, [225,226]. By summing (F.35) over permutations, one obtains the 
symmetric algebraic relations between harmonic sums. At depth 2 and 3 these read, [124], 

Sm,n ~l~ 'S'n,rn — ^m^n ~l~ ^mAn ; (F.37) 

^ ] Sl,m,n — SiSmSn + ^ ] SiSmAn + 2 Si/\mAn i (F.38) 

perm{i,m,n} inv perm{Z,m,n} 

which we used extensively to simplify our expressions. In (F.37, F.38), "perm" denotes 
all permutations and "inv perm" invariant ones. 

The limit — )> oo of finite harmonic sums exists only if Oi 7^ 1 in (F.28). Additionally, 
one defines all cr-values symbolically as 

lim Sa,,...,aXN) . (F.39) 

iV— >-oo 

The finite a- values are related to multiple (-values, [124,125,134,135,202,227,228]. Fur- 
ther we define the symbols 

00 

(70 (F.40) 
^ — ' I 

1=1 

00 

a, := (F.41) 
1=1 

It is useful to include these cr-values into the algebra, since they allow to treat parts of 
sums individually, accounting for the respective divergences, cf. also [124,125,202]. These 
divergent pieces cancel in the end if the overall sum is finite. 

The relation of single harmonic sums with positive or negative indices to the -0^^^- 
functions is then given by 

S^{N) = V(^ + 1)+7B, (F.42) 
Sa{N) = tl)ll^(a-i)(iv+l) + C„,A;>2, (F.43) 

r(a) 

S_^{N) = (-l)^/5(Ar + l) -ln(2) , (F.44) 
S-a{N) = -i_A_^(«-i)(7V + l)-(l-2i-")Ca,^>2. (F.45) 

Single harmonic sums can be analytically continued to complex values of A^ by these rela- 
tions. At higher depths, harmonic sums can be expressed in terms of Mellin-transforms of 
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polylogarithms and the more general Nielsen-integrals, [229-231]. The latter are defined 

by 



S.,,(.) = I ^ \oi^-\x) log^(l - zx) (F.46) 

and fulfill the relation 

If p = 1, one obtains the polylogarithms 

Li„(x) = S„_i,i(x) , (F.48) 

where 

Lio(x) = . (F.49) 

\ — X 

These functions do not suffice for arbitrary harmonic sums, in which case the harmonic 
polylogarithms have to be considered, [232]. The latter functions obey a direct shuffle 
algebra, cf. [131-134,202]. The representation in terms of Mellin-transforms then al- 
lows an analytic continuation of arbitrary harmonic sums to complex A^, cf. [233-236]. 
Equivalently, one may express harmonic sums by factorial series, [219,237,238], up to 
polynomials of S\{N^ and harmonic sums of lower degree, and use this representation for 
the analytic continuation to e C, cf. [132, 133,202]. 
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G Examples for Sums 



In the present calculation numerous single- to triple finite and infinite sums of an ex- 
tension of the hypergeometric type had to be calculated. For these sums, depending on 
various summation parameters, n^, the ratio of the summands, except the part containing 
harmonic sums. 



a{...,ni + 1, ■■■) 
a(...,ni,...) 



(G.l) 



is a rational function in all variables nj. Sums of this type can be represented by basic sums 
of a certain type, which are transcendental to each other and form sum- and product- 
fields, cf. [140-143] and references therein. The general form of these sums is 



Ni(N)N2{ki,N)N3{k2,ki,N) 4 

£ E E R{h,h,h,N)l[Sa;,{s{k,N))T 

ki=l fc2 = l fc3=l 1=1 



si{ki,N)...Sp{ki,N) 

Si{ki,N)...Sq{ki,N) 



, (G.2) 



with R a rational function, s{ki, N) a linear combination of the arguments with weight 
±1, CLi an index set, p,q G N, iVj e N U oo. The generahzed F-function usually includes 
both Beta-functions and binomials. 

In the present calculation one faces more complicated sums than occurring in earlier 
two-loop calculations up to 0(e), [120, 121]. We present a few examples. 



■1) 



N 



6 - 23N + 9N^ + 2N^ 
2{N -iyN^{l + N){2 + N) 



+ 



27(-l) 



N 



N + 2 (N -1)N(N + 1)(N + 2) 



Si 



N{N + 2) 
1 



•SI 



+ 



48f-l 



N + 2 {N -1)N{N +1){N + 2) 



S3S2 



2S^ 



N{N + 2) 



^ ^ ' 20{-l + Nfm{l + N){2 + N) 

7 



■(-1 



21 



10(A^ + 2) 



^1 + 



10A^(A^ + 2) 



+ 



+ 



3(-l 



+ 



(N -1)N(N + 1)(N + 2) N + 2 

^ 18(-1)^ 

N + 2 {N -1)N{N + 1){N + 2) 



Si 



N 



+ N)N{1 + N){2 + N) 
6 - 23A^ + 9A^2 ^ 2N^ 



2(-l + Nym(l + N)(2 + N) 
3 



N(N + 2) 

2 



' Sa 



s,+ 



2Sl^ ^ (-l)^(3iV- 1) 



N + 2 



{N - l)3Ar3 



Si 



2 + N 



S^2S-s +{-{-!) 



N 



3(12 - 6iV - UN^ + 7N^ + 12N* + N^) 
(-1 + N)^N^{1 + N){2 + N) 
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+(-1) 



N 



9Sf 



+ 



+ 



{-1 + N)N{1 + N){2 + N) 
3(-l)^ 1 



{N -1)N{N + 1){N + 2) N + 2 



^ 2(l2-37iV + 9iV^ + 4iV3) ^ 



^^""^^ (-1 + Ar)27V2(i + Ar)(2 + Ar) 

'S'2|'<S'2,l 

24 



5i 



(-l + iV)2iV2(l + iV)(2 + iV) 



(-l + A^)iV(l + A^)(2 + iV) 



5i 



'S'3,1 



2S- 



3,2 



iV + 2 



12(-1) 



N 



2 + AT 
-(-1)- 

+ (-1)^ 
-(-1)- 



(AT- l)Ar(iV + l)(iV + 2) N + 2 
-(-1) 



'S'4 1 



Ar + 2 



2(6-23A^ + 9Ar2 + 2Ar3) 



12 



-1 + Ar)2Ar2(l + Ar)(2 + AT) 

25'_2,l,-2 



(-1 + Ar)Ar(l + Ar)(2 + Ar) 
1 

(-1 + A^)A^(1 + A^)(2 + Ar) 
30 



'S'2,1,1 - 



N + 2 



42^2,2,1 - 24^3,1,1 + 54^2,1,1,1 



-1 + N)N{1 + N){2 + N) 



8i8i - ~8i{2) 



+(-1) 



N 



30 



(-l + iV)iV(l + iV)(2 + iV) 



8i~8i - ^3(2) + 



+ 



-(-1) 



AT 



2(6- 12N + 7N^ + N^) 
{~l + NfN^ 



+ (-1) 



TV 



{-l + NfN 

&8l 



{-l + N)N{l + N){2 + N) 



2(6-23iV + 9iV2 + 2iV3) 



(-1) 
+ (-1)^ 



(-1 + NYN'^{1 + N){2 + N) 
^ 2{Q-2'3N + m'^ + 2N^) 



81 



(-1 + A^)2Ar2(i + N)(2 + N) 
6 



+ (-1) 



TV 



12 



+ < 



+ 



(-1) 
-(-1) 

N 



(-1 + N)N{1 + N)(2 + N) 

AT -L^'Jl 



'^3(2) }8,,i 



(-1 + N)N(1 + N)(2 + N) 
1 



Si 



^2(2) 



2' 



(-1 + A^)A^(1 + A^)(2 + iV) 



2(6-23iV + 9iV2 + 2Ar3) 
^^""^^ (-1 + Ar)2iV2(l + iV)(2 + iV) 



5i 



^1 



^ 2(6-23iV + 9iV2 + 2Ar3^ 
-l + iV)2Ar2(l + Ar)(2 + AT) 

12 



-1 + Ar)A^(l + A^)(2 + AT) 



^1 



8i,i {-,!]} 8i,i{2,l) 



+ 



2(6-l2Ar + 7Ar2 + Ar3) 
^ (_l + Ar)3Ar3 ^ 



6-52 



(-1 + A^)A^(l + A^)(2 + A^) 



2(6-23Ar + 9Ar2 + 2Ar3) 
^ ^ (-l + Ar)2Ar2(l + A^)(2 + A^) ' 
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+(-1)" 
-(-1)^ 



30 



(-1 + N)N{1 + N){2 + N) 
6 



S, 



1 + N)N{1 + N){2 + Nf''' [r'^l 



+ 



X 



^ 2(6-23iV + 9iV^ + 2iV3) , ^ 

' (-l + NyN^(l + N)(2 + N) 



12 



(-l + iV)A^(l + iV)(2 + iV) 



^1 



(-1) 



N 



{-1 + N)N{1 + N){2 + N) 



30S,{-]S,,s{2,l) + 36S,,,[-,2 



+ 



+ < 



-30S,S, ( ^"j - 30^2 + 30^1,1 Q, 1 



^2,l(l,2) 



(-1) 



^ 2(6-23iV + 9iV2 + 2iV3) 
(-1 + 7V)27V2(1 + Ar)(2 + AT) 



+ (-1) 



AT 



24 



{-1 + N)N{1 + N){2 + N) 

i-ir 



Si 



Si{-]}S,,i{2,l) 



(-1 + N)N(1 + N)(2 + N) 



-12^1,1 (^,lj^2,i(2,l) + 30^2,3 (^,2 



-12^2,3 ( 2, ^) - 18^1 ^3,i(2, 1) + 30^3,2 2\ - 30^4,1 2 



+ 



(-1) 



N 



2{6-12N + 7N^ + N^) 

(-l + iV)3iV3 

^ 2(6-23A^ + 9iV2 + 2iV3^ 



-(-1) 



AT 



6^2 



-1 + N)N{1 + N){2 + N) 



-(-1) 



AT 



1 + Ar)2iV2(l + iV)(2 + AT) 
30 



Si 



So 



+ 



-(-1) 



iV 



1 + Ar)Ar(l + Ar)(2 + AT) 

2(6-23Ar + 9A^2^2A^=^) 
(-1 + A^)2A^2(i ^ ^)(2 + A^) 



-^wl ^,1,2) +^1,1,1 Q,2,l 



-(-1)^ 
+(-1)^ 



12 



(-l + A^)Ar(l + A^)(2 + A^) 
12 



^1 



Si{^]S,,iAW) 



+ 



X 



■1) 



^ 2(6-23A^ + 9A^2^2A^3) 
(-1 + A^)2A^2(i ^ ^)(2 + A^) 



■1) 



N 



12 



(-1 + A^)A^(1 + A^)(2 + A^) 



^1 



^1,1.2 (l, 1, 2] - ^1,1,2 (2, 1, 1) - 2^1,1,2 (2, 1, ^ 
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N 



(-1 + N)N{1 + N){2 + N) 



665i,i,3U'l'2 +36^1,1,3 -'2,1 



+30>5i,i,3 l,-,2 +30>Si,i,3 1,2,- 



+ 



X 



4(6-23iV + 9iV^ + 2iV3) _ 
^ ^ (-1 + Ar)2Ar2(i + Ar)(2 + AT) ^ ^ I 



24 



-l + 7V)iV(l + 7V)(2 + 7V) 



5i 



-51,2,1 (^,2,1) -5i,2,i(2,-,l)- 2,1,- 



-(-1) 



iV 



12 



(-1 + A^)A^(1 + A^)(2 + N) \2 



-(-1) 



TV 



5i - 5i,2,l(l,2,l) 
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(-1 + A^)A^(1 + A^)(2 + Ar) 



30^1,2,2 U,l,2 +365i,2,2 U,2,l 
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1 
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1 



1 



+30^2,1,2 1, 7^, 2 - 12^2,1,2 2, -, 1 - 24^2,1,2 2, 1, - + 30^2,2,1 1, 2, - 
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1 



-24^2,2,1 2, -, 1 - 12^2,2,1 2, 1, - + 30^3,1,1 1' 2 + 30^3,1,1 -, 2, 1 
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+12,Si,2,i,i -,2,1,1+ 30^1,2,1,1 1, 1, -, 2 + 30,Si,2,i,i 1, 1, 2, 



30,52,1,1,1 1, -, 1, 2 + 30^2,1,1,1 1, 2, 1 + 30^2,1,1,1 1, 1, ^, 2 



+30,^2,1,1,1 1,1,2,- + 30^2,1,1,1 1, 2, -, 1 + 30^2,1,1,1 1, 2, 1, 
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-24^2,1,1,1 ( 2, i 1, 1 j - 24^2,1,1,1 (2, 1,^' 1) - 24^2,1,1,1 (2, 1, 1, \ 



-12^1,1,1,1,1 ( \, 1, 2, 1, 1] - 365i,i,i,i,i Q, 2, 1, 1, 1 
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(G.6) 



The nested sums emerging in this work, which were not given before in Refs. [120, 
121] and those being closer related to the structure of harmonic sums [125], are of the 
type illustrated above. The latter have been calculated using C. Schneider's packages 
Sigma [140-143], EvaluatelnfiniteSums [239] and J. Ablinger's package HarmonicSums [202, 
240]. The involved structure of Eq. (G.3) requires simplifications due to explicit algebraic 
and structural relations for the generalized harmonic sums. These will be given in an 
upcoming publication Ref. [206]. 
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H Fixed Moments 



H.l 3-loop Moments to the anomalous dimensions 

In the following we list the O {TpHfCp^A) contributions to the fixed moments of the anoma- 
lous dimensions related to the present calculation. We have used them for comparison 
in recalculating the corresponding contributions to the anomalous dimensions at general 
values of AT, cf. Eqs. (194,210,218). 

■,(2) 



.... ^(2), PS 
(^^) Tqq 



'~(2),.^ rr.2 /4481539^ 9613841 \ ^, 
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7^(8) = T,n,^ 1377810000 + 810152280000 ^""^^ 
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- ^^"Al9413375^^-+ 161800390434375 ^-) ^^'^^ 



^(2)-PS(4) ^ _t2^^ 123734 

/qq V y F / i- 151875 V ) 

in = -T|„,C, (H^8) 

iqq y > F } t 72930375 ^ ' 

^ ^ 6751269000 ^ ^ 

^g),PS(io) ^ -r^„,C, 531694610144 

^ ^ ^ ^ 420260754375 ^ ' 

^g),P5(i2) = 2566080055386457 

'11 ^ ' F / i- 2851637832143100 ^ ' 



^(2),NS,+ (2) = -T>,C.^ (H.13) 
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^(2),NS,+ /«^ ^2 ^ 321390284 

7^' '^(6) = -^i^-/^- 1041862^ 

,(2),NS,+ .o^ _ ^2„ ^ 38920977797 
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,(2),NS,+ ...^ _ ^2^ ^ 27995901056887 

(10) - -T,n,C, 748828253250 
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7< 
7. 
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^ ^ 562605750 
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7,, (Idj - IpUfL^F 1645175672390250 
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H.2 The 0{£^) Contributions to A. . 



We list the contributions O {TpnfCF,A) to the fixed moments of the constant part of the 
unrenormalized massive OMEs, a-^'^^ from Ref. [95]. We used these values for comparisons 
to the general A^-results computed in the present paper, Eqs. (176,200,206,215). 
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H.3 3— loop Moments for Transversity 



We list below the contributions O {TpnfCF,A) to the fixed moments of the transversity 
anomalous dimension, cf. [179], to which we compared the result for general values of N 
calculated in the present paper, Eqs. (222), cf. also [203]. 
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Here we list the O {TpnfCF,A) moments of the constant part of the unrenormalized massive 

f3) TR. 

OMEs Ci)jqQ from Ref. [179]. Wc used these values for comparisons to the general N- 
result computed in the present paper, Eq. (221). 
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